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CHAPTER  I 


INTRODUCTION 


The  need  for  widespread  wave  data  of  improved  quality  has 
been  recognized  by  the  ocean  science  and  engineering  community  for 
some  time.  In  an  editorial  in  Shore  and  Beach  (October , 197*0  ,  M.  P. 
O'Brien  noted  that: 


Observations  and  measurements  of  ocean  waves 
have  been  made  at  points  along  the  coasts  of  the 
United  States  —  at  some  locations  over  a 
considerable  period  of  time  —  but  the  accumulated 
data  fall  short  of  the  need  in  geographical 
distribution,  duration  of  the  records,  and  detail  of 
wave  characteristics....  [Data] ...  suffers  from  a 
number  of  deficiencies,  notably  the  lack  of  wave 
direction  measurements.... 


Significant  strides  have  been  made  in  the  measurement  and 
recording  of  wave  amplitude  and  period  since  Dean  O'Brien's 
commentary,  but  the  success  of  these  systems  has  not  been  matched  by 
the  development  of  numerous  reliable  wave  sensors  with  directional 
capability.  Despite  the  importance  of  directional  wave  information 
to  offshore  and  coastal  engineers  as  well  as  ocean  scientists,  the 
development  of  directional  wave  sensors  has  lagged  behind  other  ocean 
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instrumentation  advancements.  Indeed,  a  recent  conference  of 
scientists  and  engineers  skilled  in  wave  measurement  cited  the 
improvement  of  the  measurement,  analysis,  and  reporting  technology  of 
directional  wave  spectra  as  the  first  priority  of  future  ocean  wave 
research  and  development  (Dean, 1982). 

Some  directional  wave  data  are  available  today  using  both 
remote  sensing  and  in-situ  techniques.  Remote  sensing,  a  promising 
alternative  for  the  future  because  of  its  capability  to  monitor  large 
ocean  surface  areas,  presently  suffers  from  the  expense  and  lack  of 
availability  of  air  and  space  craft  platforms  as  well  as  sometimes 
lengthy  and  expensive  data  analysis.  Further,  the  development  of 
remote  sensing  still  requires  dependable  ground  verification. 
Although  in-situ  systems  are  single  point  sensors,  they  offer  some 
advantages  over  remote  sensing  techniques.  In-situ  instruments  can 
monitor  the  sea  surface  continuously,  are  composed  of  less  expensive 
and  less  complex  hardware,  and  can  generally  offer  relatively 
inexpensive  data  reduction  that  approaches  real-time.  Although  more 
difficult  to  install  and  service,  submerged  in-situ  systems  offer 
advantages  over  surface  riding  or  piercing  instruments  in  that 
submerged  hardware  is  less  prone  to  loss,  does  not  interfere  with  the 
propagating  wave  field,  and  is  more  likely  to  operate  without 
malfunction  during  extreme  wave  events.  Such  underwater  systems, 
however,  are  generally  limited  to  use  in  near-shore  applications  or 
at  a  shallow  depth  from  a  fixed  platform  in  deep  water. 
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Several  types  of  underwater  in-s i tu  systems  have  been 
suggested;  some  of  which  are  operational  presently.  Such 
instruments  include  tri-axial  current  meters  (van  Heteren  and 
Keyser , 1981) ,  bi-axial  current  meters  with  one  pressure  transducer  or 
wave  staff  (Aubrey, 1981 ;  Nagata,1964;  Bowden  and  White, 1966; 
Simpson, 1969) ,  arrays  of  pressure  transducers  (Peacock, 1974; 
Seymour , 1978;  Mobarek , 1965 ;  Chakrabart i , 197 1 ;  Chakrabarti  and 
Snider, 1973;  Panicker  and  Borgman, 1970;  Ploeg,1972),  and 
measurement  of  the  bi-axial  components  of  the  force  on  a  sphere 
together  with  bottom  pressure  fluctuations  (Suzuki , 1969) .  A  review 
of  many  procedures  that  have  been  used  to  make  directional  wave 
measurements  is  presented  by  N.  N,  Panicker  (1974),  K.  Rikiishi 
(1977),  and  more  recently  within  the  Proceedings  of  the  National 
Reasearch  Council  Workshop  on  Wave  Measurement  Technology  (1982) . 

Data  from  some  current  meter  systems  is  often  suspect  after 
deployment  because  of  the  suseptibility  of  the  instrumentation  to 
bias  from  bio-fouling  and  corrosion.  Bottom-mounted  pressure 
transducers,  on  the  other  hand,  have  been  used  with  considerable 
success  for  wave  sensing  for  several  years  (Forristal 1 , 1981) . 

The  concept  of  measuring  the  directional  characteristics  of 
waves  is  fairly  straightforward.  If  one  considers  a  long  wave  crest 
approaching  parallel  to  the  alignment  of  two  pressure  sensors,  as  in 
Figure  I  —  1 ,  it  is  easy  to  imagine  both  sensors  recording  the  crest 
passage  simultaneously.  If  the  wave  crest  approaches  obliquely  to 


the  pressure  sensors,  as  in  Figure  I-2,  one  could  imagine  that  one 
sensor  would  report  the  passage  of  the  crest  before  the  other.  The 
difference  in  pressure  over  time  between  the  sensors'  records  can  be 
analyzed  by  well-known  techniques  and  a  first  estimate  of  the  wave 
direction  calculated,  (Longuet-Higgins , Cartwright , and  Smith, 1963). 
There  is,  of  course,  a  directional  ambiguity  with  only  two  sensors 

since  a  wave  approaching  as  in 
Figure  I~3  would  create  the  same 
pressure  record  over  time  as  in 
Figure  I~2 .  Such  ambiguities  may 
be  resolved  by  additional 
sensors . 

Presently,  the  analysis 
of  data  from  bottom-mounted 
arrays  of  pressure  sensors 
usually  develops  pressure 
differences  through  the  sub¬ 
traction  of  individually  meas¬ 
ured  absolute  pressures. 
However,  instantaneous  differ¬ 
ences  in  the  transducers'  records 
due  to  wave  action  are  typically 
very  small  when  compared  to  the 
large  pressure  values  that  they 


are  required  to  record.  Hence,  the  pressure  difference  across  an 
array  is  calculated  by  subtracting  two  large  numbers  in  order  to 
generate  a  very  small  number.  Such  a  system  lends  itself  to  inherent 
inaccuracies.  For  example,  a  bottom-mounted  transducer  in  twenty 
feet  of  water  records  the  passage  of  a  three  foot  wave  crest  of  nine 
second  period  with  a  value  of  24.78  psia  compared  to  the  23.64  psia 
reading  under  still  water.  A  similar  transducer  placed  twenty  feet 
away  at  a  45  degree  angle  to  the  wave  crest  would  report  24.69  psia 
at  the  moment  the  crest  passes  its  neighboring  gauge.  (Linear  wave 
theory  is  utilized  for  this  example.)  The  arithmetic  difference 
calculated  between  the  transducers,  0.09  psia,  is  less  than  one  half 
of  one  per-cent  of  the  still  water  value  each  transducer  would  record 
or  less  than  nine  percent  of  the  dynamic  pressure  caused  by  the 
passing  wave.  The  difference  can  be  improved  by  increasing  the 
distance  between  transducers  —  at  the  expense  of  a  more  physically 
unmanageable  array,  greater  error  in  the  assumption  of  linear  water 
surface  slope,  and  introduction  of  directional  ambiguities  for  the 
higher  frequencies  present. 

It  is  the  contention  of  this  thesis  that  the  accuracy  of  the 
directional  wave  spectrum  calculated  from  such  small  differences 
between  large  pressure  values  is  questionable.  The  accuracy  and 
physical  size  of  the  pressure  sensor  array  can  be  improved  by 
utilizing  differential  pressure  transducers  designed  to  record  small 
pressure  differences  between  two  points.  This  thesis  will  document 
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the  design,  installation,  and  performance  of  such  an  instrument  as 
carried  out  at  the  University  of  Delaware  and  the  Coastal  Engineering 
Research  Center's  test  and  field  research  facilities  in  Fort  Belvoir, 
Virginia,  and  Duck,  North  Carolina. 

Since  the  majority  of  industries  and  machinists  associated 
with  the  project's  instrumentation  work  in  the  English  system  of 
units  in  the  United  States,  the  English  system  was  used  in  the 
development  and  reporting  of  this  project. 


CHAPTER  II 


CONCEPTS  OF  A  DIFFERENTIAL  PRESSURE  GAUGE  SYSTEM 


A.  DPG  Response  to  Surface  Waves 

As  previously  mentioned,  it  is  the  difference  in  pressure 
between  points  that  must  be  considered  in  the  determination  of  wave 
direction.  Differential  pressure  transducers  generate  an  electrical 
voltage  proportional  to  the  fluid  pressure  difference  on  opposite 
sides  of  a  mechanical  diaphragm.  Absolute  (or  gage)  pressure 
transducers  report  only  the  fluid  pressure  at  one  point.  Using 
absolute  (or  gage)  transducers,  the  difference  in  pressure  between 
two  points  must  be  arithmetically  created  by  subtracting  the 
pressures  reported  at  each  of  the  two  points  of  interest.  Since  the 
differential  transducer  measures  pressure  differences  directly,  the 
differential  gauge  can  be  considered  to  be  an  inherently  more 
effective  instrument  for  determining  wave  direction.  Furthermore, 
the  response  with  depth  of  the  differential  pressure  is  a  maximum 
over  a  range  of  frequencies  that  is  more  representative  of  typical 
ocean  gravity  waves.  For  a  pressure  transducer  located  at  height  s 
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above  bottom,  the  dynamic  pressure  at  the  sensor  can  be  expressed 
over  time  t  as: 


PdynU'Y's;t)  =  Y  I  Kp  COstkxx  +  *yy 


at  +  e  ) 


(2.1) 


where  k 

P 


cosh  ks 
cosh  kh 


and  H  m  wave  height 

kx“  k  cosQ  *  wavenumber  in  the  x-direction 
ky=  k  sinO  *=  wavenumber  in  the  y-direction 
Q  =  wave  direction 
O  =  wave  frequency 
h  ■  water  depth 
Y  -  specific  weight  of  water 
e  “  phase  angle 

The  components  of  differential  pressure  differs  in  phase  and  by  a 
factor  of  wavenumber: 


dP  (x,  t, s ;  t)  =  -  Y  2  k 


-  at  +  e  ) 


cos  0 
sin  0 
3x 
3yJ 


K  sin  (k  x  +  k  y 
P  x  y 


(2.2) 
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A  current  meter  similarly  located  at  a  depth  s  will  record  a 
signal  under  waves  such  as: 


U(x,y,s;t)  =  —  0  K  cos  (k  x  +  k  y 
2  u  x  y 


at  +  e  ) 


where  : 


K 

U 


cosh  ks 
sinh  kh 


(2.3) 


Although  the  magnitude  of  both  the  pressure  and  current 
decreases  steadily  over  increasing  frequency,  the  differential 
pressure  signal  reaches  a  maximum  at  an  intermediate  frequency.  This 
is  illustrated  graphically  in  Figure  II-l,  which  depicts  the  dynamic 
response  of  pressure  and  differential  pressure  gauges  located  near 
the  bottom  with  a  wave  crest  of  2.5  foot  amplitude.  In  the 
upper-most  curves,  the  differential  gauge  is  assumed  to  sample  two 
points  separated  by  three  feet  coincident  with  the  wave  ray.  The 
lower  curve  represents  the  difference  in  pressure  monitored  by  two 
absolute  pressure  gauges  separated  by  twenty  feet  along  the  wave  ray. 
Each  curve  is  normalized  by  typical  values  of  the  rated  capacity  of 
the  instruments,  (35  psia  for  the  pressure  transducer,  and  +0.5  psid 
for  the  differential  transducers).  Figure  I I — 1  represents  a  water 
depth  of  fifteen  feet,  while  Figures  1 1  —  2  and  1 1—3  represent  twenty 
and  twenty-five  feet  water  depths  respectively.  In  each  case,  the 
transducers  are  located  three  feet  off  the  seafloor.  It  is  clear 
that  the  differential  pressure  signal  utilizes  a  greater  portion  of 
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FIGURES  11-1,2:  Instrument  response  to  a  five- 
foot  surface  wave  in  varying  water  depths. 
Sensors  are  three  feet  from  the  seafloor. 
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FIGURE  II-4:  Differential  preaaure  transducer 
response  to  a  five-foot  surface  wave  in 
twenty-five  feet  of  water  with  varying 
deployment  depths. 
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the  instrument's  dynamic  range  and  reaches  its  greatest  response 
level  for  waves  between  four  and  six  second  periods.  Sensitivity  to 
longer  period  waves  increases  with  increasing  water  depth. 
Differential  pressure  response  increases  as  the  transducers  are 
raised  off  the  seafloor,  as  expected,  and  the  greatest  sensitivity 
shifts  to  waves  of  higher  frequency  (Figure  II~4) . 

In  Figures  I I— 1  through  II-4,  the  response  of  the 
differential  transducer  was  developed  over  a  pressure  sample  spacing 
of  only  three  feet  —  less  than  one-sixth  of  the  gage  length  of 
conventional  pressure  sensor  array  systems.  Typical  arrays  generally 
develop  two  orthogonal  measurements  of  water  surface  slope  by 
sampling  points  separated  by  twenty  to  thirty  feet.  The  differential 
pressure  gauge  concept,  then,  suggests  that  more  efficiently  measured 
directional  information  can  be  obtained  with  a  considerably  smaller 
instrument  than  is  presently  used. 

B.  The  Differential  Pressure  Gauge  Directional  Wave  Monitor 

To  test  the  effectiveness  of  an  array  using  differential 
pressure  transducers,  the  differential  pressure  gauge  directional 
wave  monitor  (DPG)  was  developed.  Using  small  distances  between 
differential  pressure  sampling  points,  it  was  possible  to  design  a 
system  which  directly  measures  the  water  surface  slope  yet  is 
physically  smaller,  more  manageable,  and  potentially  more  accurate 


FIGURE  II-5:  Schematic  representation  of  the  DPG. 

The  DPG  samples  the  pressure  about  its  center  using  an 
absolute  pressure  transducer  and  simultaneously  samples  the 
differential  pressures  on  the  bow,  port,  aft,  and  starboard  sides  of 


the  instrument  center  using  four  differential  pressure  gauges, 
(Figure  I I — 5 ) .  The  differential  gauges  directly  infer  water  surface 
slope,  and  since  there  are  two  such  measurements  along  two  axes 
through  the  instrument  center,  it  is  possible  to  arithmetically 
develop  the  water  surface  curvature  along  each  axis.  It  will  be 
demonstrated  in  Chapter  IV  that  these  measurements  permit  the 
estimation  of  the  first  seven  Fourier  coefficients  of  the  directional 
wave  spectrum.  Conventional  analysis  from  submerged  pressure  sensor 
arrays  arithmetically  develops  one  measurement  of  the  water  surface 
slope  along  each  of  two  axes.  Typical  surface  riding  buoys  similarly 
measure  the  slope  along  each  of  the  two  axes.  This  data,  along  with 
the  water  surface  displacement  record,  permits  the  generation  from 
these  systems  of  the  first  five  directional  Fourier  coefficients. 
The  cloverleaf  buoy,  a  cluster  of  three  surface  riding  buoys,  is 
capable  of  estimating  the  first  nine  Fourier  coefficients, 
(Mi tsuyasu, 1973 ;  Cartwright  and  Smith, 1964). 


It  can  be  shown  that  it  is  possible  to  generate  up  to  the 
first  M(M-1)+1  directional  Fourier  coefficients  for  an  array  with  M 
number  of  gauges  (Borgman, 1969) •  Using  this  technique,  pressure 
sensor  arrays  such  as  the  Scripps  Sxy  gauge  (Seymour , 1978)  could 
generate  the  first  nine  directional  Fourier  coefficients.  The  DPG, 
in  the  configuration  described,  could  develop  the  first  eleven 
coeff icients .  The  DPG  could  potentially  develop  the  first  twentyone 
coefficients  with  the  same  number  of  transducers  in  a  different 
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configuration.  This  analysis  technique  is  discussed  in  Appendix  E. 


C.  Selection  of  Transducers  and  Gage  Length 


The  selection  of  the  gage  length  between  the  two  points  to  be 
compared  by  the  differential  pressure  transducers  is  dept  dent  upon 
three  criteria:  (1)  the  characteristics  of  the  transducer,  (2)  the 
error  in  approximating  the  water  surface  slope  as  a  linear  function 
between  two  sampled  points,  and  (3)  reasonable  size  limitations  of 
the  instrument.  The  third  criterion  limits  the  gage  length  per 
transducer  to  a  maximum  of  about  five  feet  if  one  imposes  a  design 
constraint  of  easy  instrument  manageability.  At  such  small  spacings, 
the  maximum  error  in  linearly  approximating  the  water  surface  slope 
between  sample  points  is  less  than  one  and  one-half  percent  for  the 
shortest  waves  of  interest,  (say,  3.2  seconds). 


Having  satisfied  the  second  and  third  criteria 
criterion  remains.  In  the  selection  of  the  differential 
one  must  consider  the  physical  size,  ruggedness, 
availability  of  a  transducer  that  can  measure 
bidirectional  pressure  fluctuations.  The  desired  rated 
the  transducer  is  a  function  of  the  maximum  wave  height 
the  ability  of  the  transducer  to  detect  and  report 
difference  of  a  minimum  wave  condition  beyond  the 
electrical  noise  level.  One  might  determine  the  rated 
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first  selecting  an  instrument  deployment  depth  according  to  the  site 
and  the  wave  frequencies  that  are  of  most  interest  or  are  most  likely 
to  be  encountered  —  remembering  that  deeper  installations  require 
more  sensitive  transducers.  A  maximum  wave  height  at  the  most 
sensitive  frequency  for  the  chosen  depth  is  then  considered  with  the 
manageable  gage  length  restriction  in  mind.  The  maximum  differential 
pressure  is  then  calculated  from  Equation  2.2.  The  result  is  an 
estimate  of  the  rated  capacity  of  the  transducer  needed.  This 
estimate  is  then  used  to  select  a  transducer  available  from  industry. 
One  next  considers  the  noise  level  of  the  system  to  determine  whether 
the  selected  instrument  is  capable  of  reporting  the  minimum 
differential  pressure  of  interest,  (i.e.,  a  small  amplitude  wave  of 
frequency  higher  or  lower  than  that  of  the  instrument's  most 
sensitive  range  at  the  chosen  depth).  Satisfied  with  the  results, 
the  gage  length,  ax,  is  fine-tuned  based  on  the  rated  capacity  of  the 
transducer,  R,  and  the  maximum  wave  of  interest.  From  Equation  2.2, 


Ax 


_ R _ 

Y  H  k  cos  0  K 
P 


(2.  A) 


The  instrument's  response  to  a  minimum  wave  condition  is  checked 
again  using  the  newly  calculated  gage  length  to  ensure  that  a 
reasonable  signal  to  noise  ratio  is  maintained  for  small  amplitude, 
long  period  waves. 
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The  installation  site  of  the  DPG  was  proposed  as  the  United 
States  Army  Corps  of  Engineers  Field  Research  Facility  (FRF)  at  Duck, 
North  Carolina,  operated  by  the  Coastal  Engineering  Research  Center 
(CERC) .  The  instrument  was  to  be  placed  near  and  hard-wired  to  the 
Facility  research  pier.  After  examining  the  bathymetry  near  the  pier 
and  considering  the  high  cost  of  underwater  cable,  a  nominal 
deployment  depth  of  20  feet  was  selected.  This  indicated  that  the 
instrument  would  be  most  sensitive  to  waves  of  about  5  second  period, 
(Figure  I I— 2) .  This  was  considered  acceptable  at  the  time  since  the 
mean  monthly  nearshore  wave  periods  for  this  site  were  reported 
between  A. 8  and  6.5  seconds  (SPM,  Figure  4-11,  1977).  (It  was 
discovered  much  later  that  this  estimate  was  poor.  A  better  estimate 
is  given  as  about  7.5  seconds  from  Birkemeier,  et.al.,  (1981).)  A 
design  wave  height  of  16  feet  was  chosen  and  considered  with  a  gage 
length  of  five  feet.  This  suggested  the  use  of  a  +0.6  psid 
transducer  which  was  unavailable.  A  +0.5  psid  transducer  was 
selected  instead,  mandating  a  gage  length  of  4.15  feet  (1.27  meters) 
using  Equation  2.3. 

The  overall  dimensions  of  the  final  instrument,  then,  became 
9  feet  -  9  inches  (2.97  m)  along  each  axis  (including  protective  caps 
on  the  ends  of  each  arm)  and  40  inches  (1  ro)  in  height  to  accomodate 
the  electronics  package. 
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The  response  of  an  instrument  of  this  configuration  in 
varying  depths  of  seawater  is  illustrated  in  Figure  II-6.  The 
contours  represent  lines  of  equal  transducer  response  per  ^ feet  of 
wave  height  normalized  by  the  rated  capacity  of  the  transducer.  The 
fields  of  highest  value  indicate  the  environmental  conditions 
conducive  to  greatest  instrument  performance. 

The  absolute  transducer  was  selected  with  a  rated  capacity  of 
50  psia  (because  of  cost  and  availability)  which  provides  the  DPG  the 
capability  to  measure  a  pressure  head  of  over  75  feet. 


DEPTH  f  feet  ] 


FREQUE  NCY  [rad/sec] 
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FIGURE  I 1—6 :  Differential  pressure  response  with  water  depth  and 

frequency.  Contours  are  lines  of  equal  dP  response  normalized  by 
the  rated  capacity  of  the  instrument,  (1.0  psid).  (Sensors  are 
3.875  feet  from  the  seafloor;  differenial  gage  length  *  4  feet 

along  wave  ray;  wave  height  “  5  feet.) 
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CHAPTER  III 
DPG  HARDWARE 

A.  INSTRUMENT,  SERVICE  BOX,  AND  CRADLE 

1.  Introduction 

The  in-situ  instrumentation  is  contained  in  a  poly-vinyl 
chloride  (PVC)  structure  that  is  secured  to  a  steel  cradle  moored  to 
the  seafloor.  The  PVC  structure,  or  "instrument",  a.',  it  is  hereafter 
called,  consists  of  a  central  tube,  or  "fuselage",  and  four  "arms" 
that  extend  from  near  the  top  of  the  fuselage.  A  plexiglass  service 
box,  also  attached  to  the  cradle,  provides  storage  for  eighty  feet  of 
unarmored  cable  at  the  installation  site.  Figure  1 1  I  —  1  illustrates 
the  instrument  and  Figure  1 1 1  —  2  depicts  the  steel  cradle  with  the 
instrument  in  place  suspended  under  the  Field  Research  Facility's 
CRAB  vehicle. 

This  chapter  describes  the  DPG  system  as  installed  at  the 
FRF,  Duck,  N.C.  Several  components  of  the  system  were  altered  from 
the  original  design  during  DPG  assembly  and  laboratory  testing.  The 
discussion  of  these  modifications  is  presented  in  Chapter  V. 
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FIGURE  III-2:  The  DPG  instrument  secured  to  the  steel  cradle.  The  assembly  is 
suspended  beneath  the  Coastal  Research  Amphibious  Buggy  (CRAB) 
operated  by  the  CERC  Field  Research  Facility. 


FIGURE  I I 1-3 :  Rendering  of  the  DPG  fuselage, 
water-tight  instrumentat i on  cylinder. 
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2.  The  Instrument  Fuselage  and  Arms 

The  fuselage  section  of  the  instrument  is  ten  inches  in 
diameter  and  forty  inches  in  length,  (Figure  III-3).  It  secures  a 
circular  plate  with  five  pressure-sensing  "isolation"  diaphragms  and 
a  water-tight  cylinder  that  contains  the  pressure  transducers.  Both 
ends  of  the  fuselage  are  covered  by  plexiglass  plates  to  prevent 
large  marine  animals  from  entering  the  inside  of  the  instrument.  The 
top  plate  contains  a  number  of  5/8  inch  holes  to  allow  the  fuselage 
to  flood  and  ensure  that  the  pressure-sensing  diaphragms  therein  are 
communicating  with  the  outside.  The  bottom  plate  has  a  one  inch  wide 
slit  cut  across  its  radius  to  allow  the  cable  entry  to  the  fuselage 
and  water-tight  cylinder.  Both  plates  are  secured  to  the  cylinder  by 
titanium  bolts  attached  to  PVC  blocks  on  the  inside  of  the  fuselage. 
The  blocks  are  made  fast  to  the  fuselage  by  adhesives  and  nylon 
screws . 


Each  of  the  four  arms  penetrate  the  fuselage  and  are  held 
therein  by  a  PVC  "spider."  A  spider  is  a  common  plumbing  junction 
which  connects  four  pieces  of  tubing  at  ninety  degree  angles  to  one 
another.  Four  bolts  are  tapped  through  the  spider  to  secure  the 
arms.  In  this  way,  any  or  all  of  the  arms  can  be  disconnected  for 
greater  ease  in  transportation,  or  potentially,  a  change  in  arm 


orientation. 
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The  arms  are  three-inch  diameter  PVC  pipe  and  fifty-two 
inches  in  length.  Each  contains  a  pressure-sensing  diaphragm  near 
its  end.  The  extreme  end  of  each  arm  is  covered  by  a  threaded 
pipe-cap  which  can  be  removed  to  inspect  the  diaphragm  and  replace 
copper  scouring  pads  placed  near  the  diaphragm  to  reduce  biological 
fouling.  The  arms  are  punctured  with  11/16  inch  diameter  holes  about 
their  ends  to  allow  diaphragm  communication  with  the  environment. 
Each  arm  is  labelled  with  strips  of  white  cloth  tape  to  identify  the 
differential  transducer  channel  corresponding  to  that  arm. 

3.  The  Transducers 

A  typical  differential  pressure  gauge  senses  the  pressure 
difference  between  the  diaphragm  at  the  end  of  an  arm  and  one  of  the 
five  diaphragms  mounted  inside  the  fuselage.  The  differential 
pressure  transducers  are  located  inside  the  water-tight  cylinder. 
The  differential  transducers  are  Setra  Systems  Model  228.  They  are  a 
high  line,  low-di f ferent ial  wet-wet  capacitance  type  sensor  with 
infinite  resolution  and  a  ±0.5  psid  range.  Transducer  output  noise 
is  rated  as  below  100  microvolts  RMS.  Burst  pressure  is  2500  psid 
(either  side)  such  that  there  is  less  danger  of  catastrophic  failure 
during  construction  or  installation.  The  differential  transducers 
require  28  volts  DC  nominal  excitation  and  feature  0  to  5  volt 
positive  output.  The  transducer  circuits  have  internal  protection 
against  reversed  excitation  voltage  for  at  least  5  minutes, 
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short-circuit  of  signal  output  leads,  and  short  duration  power  line 
transients  up  to  150  volts.  Each  is  supplied  with  independent  remote 
zero  adjustment,  although  this  feature  is  not  used  in  the  present  DPG 
model . 


The  fifth  pressure-sensing  diaphragm  in  the  fuselage  supplies 
a  signal  measured  by  an  absolute  pressure  transducer  also  located 
inside  the  cylinder.  The  absolute  gauge  is  a  Setra  Systems  Model 
205-2  of  50  psia  range.  It  also  operates  on  28  volt  DC  nominal 
excitation  for  a  0  to  5  volt  positive  output. 

A.  The  Isolation  Sensor  Diaphragms 

The  stainless  steel  diaphragms  of  the  transducers  are 
protected  from  the  harsh  marine  environment  by  the  nine  "isolation" 
diaphragms  located  within  the  fuselage  and  arms,  (Figure  III-A).  The 
isolation  diaphragms  are  made  of  13  mil  DuPont  Fairprenee  elastomer 
mounted  on  an  acrylic  housing.  Fairprene®  is  a  durable  nylon 
material,  coated  with  neoprene,  that  is  flexible  perpendicular  to  the 
plane  of  the  fabric.  The  elastomer  is  sealed  to  its  acrylic  housing 
by  a  90-10  copper-nickel  alloy  ring.  The  90-10  alloy  was  picked  for 
its  anti-fouling  properties.  Bio-fouling  across  the  diaphragm  ring 
could  puncture  the  elastomer  or  restrict  its  ability  to  deflect  and 
thereby  introduce  a  bias.  The  rings  are  designed  to  last  at  least 
five  years  assuming  a  corrosion  rate  of  1.0  mpy,  (Dexter , 1979) .  The 
ring  is  secured  to  the  diaphragm  housing  with  six  Monel  1/A  inch 
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bolts.  Diaphragm  housings  are  2  —  7/S  and  3  inches  in  outside  diameter 
within  the  arms  and  fuselage  respectively,  and  2-5/16  inches  long, 
including  the  Cu-Ni  ring.  The  diameter  of  the  elastomer  exposed  by 
the  ring  is  l~5/8  inches  on  all  diaphragms.  The  acrylic  housings 
behind  the  elastomers  are  hollowed  to  a  conical  shape  that  funnels  to 


FIGURE  I I I— A :  Assembly  drawing  of  an  isolation  sensor  diaphragm. 


FIGURE  I I 1—5 :  Isolation  sensor  diaphragms  shown  at¬ 
tached  by  flexible  and  nylon  tubing  through  the  top 
of  the  water-tight  instrumentation  cylinder  and  into 
the  transducers. 


a  small  opening  in  the  back.  This  opening  mates  with  1/4  inch  (1/8 
inch  I.D.)  flexible  teflon  tubing  which  leads  to  the  sensing  ports  on 
the  pressure  transducers  inside  the  water-tight  cylinder.  The 
flexible  tubing,  manufactured  by  Cajon®,  is  armored  with  stainless 
steel  braid.  Specifically,  it  connects  to  the  back  of  the  isolation 
diaphragm  chamber  and  is  fastened  to  the  top  of  the  water-tight 
cylinder  with  Swage-Lock®  connectors,  (Figure  III-5).  The  connectors 
penetrate  the  top  of  the  cylinder  and  connect  to  nylon  tubing  inside. 
Nylon  tubing  is  used  inside  the  cylinder  because  of  its  partial 
flexibility  and  transparency.  (Nylon  is  easily  impregnated  and 
hardened  by  seawater  and  is  used  outside  the  cylinder  only  as  a 
permanent  fastener,  (Dexter , 1979) . )  The  nylon  tubing  is  heated  and 
bent  to  shape  to  mate  with  the  transducer  pressure  ports.  The 
housing  and  tubing  assembly  is  filled  with  fluid  which  transmits  the 
deflections  of  the  elastomer  to  the  stainless  steel  transducer 
diaphragms.  The  fluid  is  dyed  with  food  coloring  for  ease  of 
identification  and  inspection  for  air  bubbles  that  might  be  seen  in 
the  nylon  tubing  or  acrylic  housings. 


During  assembly,  the  isolation  diaphragm  housings  are  drawn 
through  the  arms  and  secured  inside  them  by  brass  screws  which 
penetrate  the  arm  walls  from  the  outside.  The  five  isolation 
diaphragms  inside  the  fuselage  fit  snugly  through  three-inch  holes  in 
a  circular  plate  positioned  just  above  the  spider.  A  two-inch  hole 
is  cut  in  the  center  of  the  plate  for  ease  in  assembly. 


5.  The  Water-Tight  Instrumentation  Cylinder 


The  water-tight  cylinder,  located  underneath  the  spider,  is 
an  anodized  aluminum  pressure-tested  housing  developed  by  Sea-Data 
Corporation.  Four  sacrificial  zinc  anodes  are  positioned  about  the 
top  of  the  cylinder.  All  of  the  pressure  tubing  penetrates  the  top 
(removable)  section  of  the  cylinder,  while  the  power  and  data  cable 
connects  to  the  bottom  via  Brantner  underwater  XSL-20  connectors, 
(Figure  III-6).  The  cylinder  is  held  in  place  within  the  fuselage  by 
PVC  blocks  bolted  to  the  fuselage  walls  with  Monel  fasteners. 

The  transducers  are  held  inside  the  cylinder  by  a  series  of 
five  wafers  separated  by  acrylic  rods.  The  bottom  two  wafers  support 
two  of  the  differential  transducers,  the  middle  two  support  the 
remaining  two  differential  and  one  absolute  transducer,  and  the  upper 
attaches  to  the  underside  of  the  cylinder  top  and  suspends  the  wafer 
assembly.  The  wafer  and  transducer  stack  is  shown  in  Figure  I I I— 7 • 

6.  The  Cable  and  Service  Box 

The  instrument  is  cabled  to  shore  for  power  requirements  and 
data  delivery.  Eleven  hundred  feet  of  Blake  BC  4960-1  double  armored 
cable  is  laid  between  the  instrument  and  the  land-based  recording 
facility.  The  cable  is  an  eighteen  conductor  (18  gage  -  tin/copper) 
strand  wrapped  around  a  polypropylene  filler.  Th$  conductors  are 
water-blocked,  wrapped  in  mylar  tape,  and  sheathed  under  a 


polyurethane  jacket  with  double  steel  armoring  and  an  external 
polyurethane  jacket.  The  cable  weighs  0.98  pounds  per  foot  in  air. 

Eighty  feet  of  the  seaward  end,  stripped  of  its  armoring,  is 
stored  in  an  18  inch  square  by  9  inch  plexiglass  service  box  just 
before  mating  with  the  instrument.  The  top  plate  of  the  box  is 
removable  by  unthreading  four  titanium  nuts.  This  accesses  the  spare 
cable  if  the  instrument  is  to  be  moved  or  taken  to  the  surface.  If 
the  instrument  requires  service  in  the  boat  or  on  shore,  the  cable 
may  be  disconnected  above  water,  sealed  with  a  dummy  connector,  and 
released.  A  dummy  connector  is  then  attached  to  the  instrument  for 
safety.  For  re-installation  at  some  later  time,  the  cable  is 
recovered,  re-connected  in  the  boat,  and  the  instrument  is  brought 
underwater.  It  is  hoped  that  operating  the  underwater  connectors 
only  above  water  will  help  alleviate  some  of  the  flooding  and 
corrosion  problems  experienced  in  the  past  with  underwater-pluggable 
connectors . 

7.  The  Electrical  System 


Fourteen  of  the  available  twenty  conductors  are  utilized. 


Ten  carry  the  five 

transducers'  signal 

outputs , 

two 

carry 

the 

positive  excitation  voltage  (one 

for  back- 

up) ,  and 

two 

carry 

the 

negative  excitation 

ground, 

(one  for 

back-up) . 

Each  of 

the 

transducers  share  the 

excitat ion 

pos i t i ve 

and  negative 

leads . 

A 

simple  RC  network  can  be  installed  at  the  landward  end  of  each  signal 
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output  as  a  passive  low-pass  filter.  The  network  is  easily  designed 
for  3  dB  attenuation  at  50  cycles  per  second  with  the  aim  of 
eliminating  60  hertz  noise  from  the  power  supply  and  local  electrical 
systems.  Such  high  frequency  noise  would  alias  into  the  calculated 
water  wave  spectra  and  bias  the  analysis  results.  An  active  filter 
is,  however,  more  strongly  recommended.  To  guard  against  the 
possibility  of  unwanted  output  oscillation  caused  by  the  capacitance 


introduced 

by 

the  long 

cable 

to  shore,  100  ohm  resistors 

are 

installed 

in 

series  in 

each 

of  the  output  leads  inside 

the 

water-tight 

cylinder. 

Since 

the  electrical  circuitry 

of 

the 

transducers 

is 

equivalent 

to  a 

A-terminal  network  which 

can 

be 

grounded  at  only  one  point,  it  is  essential  that  the  negative  output 
leads  of  the  transducers  not  be  commoned  since  the  negative 
excitation  leads  are  already  commoned  inside  the  instrument. 

A  complete  wiring  diagram  is  presented  in  Figure  III~8. 

8.  The  Instrument  Cradle 

Both  the  instrument  and  the  service  box  are  mounted  onto  a 
steel  cradle  secured  to  the  seabed.  The  cable  is  constructed  of 
C4x5.A  steel  channel  and  3/A  inch  thick  steel  plate,  two  inch  steel 
pipe  legs  and  Llxlxl/8  angle  iron  stringers.  The  channel  is 
connected  primarily  by  the  steel  plate.  A  break  is  made  in  one  plate 
to  allow  the  cable  to  slip  through  the  structure  when  the  instrument 
is  pulled  out  of  the  cradle.  This  is  closed  by  a  plate  which  bolts 
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over  this  area  to  retain  structural  integrity  after  the  instrument  is 
installed.  The  legs  are  welded  onto  the  underside  of  the  channel 
such  that  they  balance  the  cradle  about  its  center  when  no 
instrumentation  is  installed  upon  it.  Onto  the  bottom  of  the  forty 
inch  legs  are  welded  10  by  12  inch  thin  steel  plates  to  reduce  burial 
into  the  seafloor.  The  cradle  is  secured  to  the  bottom  using  5/16 
inch  galvanized  chain  that  attaches  underneath  each  channel  and 
cinches  tight  with  a  chain  binder  to  one  of  four  screw  anchors  driven 
into  the  seafloor. 

The  instrument  arms  lay  into  the  steel  channel  and  the 
fuselage  slips  into  the  center  of  the  cradle.  The  arms  are  made  fast 
to  the  cradle  with  heavy  electrical  cable  ties.  The  sevice  box  is 
mounted  onto  the  cradle  by  four  titanium  bolts  secured  through  the 
bottom  of  the  box.  Attaching  plates  and  rings  are  welded  throughout 
the  cradle  for  buoy  markers,  flotation  devices,  pingers,  diver-assist 
lines,  and  weights  to  balance  the  asymmetrically  placed  service  box 
if  necessary. 

9.  Biofouling  and  Corrosion  Protection 

The  cradle  was  sandblasted,  primed,  undercoated,  and  then 
covered  with  anti-fouling  paint.  All  surfaces  of  the  plexiglass 
service  box  and  the  PVC  instrument  structure  (with  the  exception  of 
the  fuselage  end-cover  plates)  were  sanded  and  also  coated  with 
anti-fouling  paint.  The  paint  used  was  red  PETIT  BIOGUARD  1665  with 
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19- 1 A 71  active  ingredient  bis  (tributyltin)adipate. 

B.  Data  Retrieval 

At  the  FRF  installation  site,  the  cable  interfaces  with 
Facility  equipment  at  the  seaward  end  of  the  pier.  The  power  supply 
board  is  located  inside  a  trailer  at  the  seaward  end  and  the 
transducer  output  leads  are  hard-wired  down  the  length  of  the  pier 
and  recorded  on  a  NOVA  computer  system  located  inside  the  facility's 
main  building  on  shore.  The  computer  digitally  records  the 
transducer  signals  at  1/4  second  sampling  for  a  17  minute  period 
every  six  hours. 

C.  Some  Considerations  of  the  Hardware  Design 
1.  Fuselage  Isolation  Diaphragms 

The  absolute  pressure  measurement  and  one  side  of  each  of  the 
differential  pressures  are  all  to  be  taken  at  the  instrument  center. 
Five  individual  diaphragms  in  the  fuselage  are  employed  for  this 
task.  Since  not  all  of  the  diaphragms  can  occupy  the  center  of  the 
instrument  in  one  plane,  there  is  a  small  error  inherent  when  saying 
that  the  differential  measurements  are  exactly  adjacent  and  share  an 
endpoint  with  the  absolute  measurement.  The  midpoint  of  the  fuselage 
isolation  diaphragms  are  each  about  two  inches  (7.1  cm)  from  the 
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actual  fuselage  center.  Further,  the  arrangement  of  the  dP3,  dP2  and 
dP4  diaphragms  within  the  fuselage  is  such  that  they  are  not 
precisely  aligned  with  the  x-  or  y-  axes  respectively.  This  design 
was  the  result  of  space  allocation  problems  within  the  top  of  the 
fuselage.  The  errors  are  such  that  the  sensor  alignment  of  dP2  is 
91°,  178°  for  dP3,  and  269°  for  dP4  with  respect  to  dPl  “  0°.  It  was 
considered  that  these  one  and  two  degree  errors  were  negligible  for 
the  present  investigation. 

2.  Height  of  Sensors 

Once  again  due  to  practical  space  allocation,  the  instrument 
was  designed  such  that  the  center  isolation  diaphragms  are  nine 
inches  (22.9  cm)  above  the  arm  diaphragms.  The  pressure  response 
function  is  developed  during  data  analysis  for  the  mean  of  these  two 
elevations  above  the  seafloor.  In  20  feet  of  water,  the  error 
between  the  pressure  response  function  for  the  actual  sensor  height 
and  the  mean  sensor  height  is  less  than  two-and-one-half  percent  for 
the  highest  wave  frequencies  of  interest  (periods  of  3.14  seconds) 
and  less  than  one  percent  for  wave  periods  greater  than  four  seconds. 
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D.  Cost 

The  hardware  used  to  construct  and  install  the  final  DPG 
prototype  at  the  Field  Research  Facility  cost  approximately  $13,000. 
Machine  shop  fees  for  the  entire  project  were  $3090.  The  total 
project  cost,  including  design,  development,  laboratory  testing, 
hardware,  travel,  computer  time,  salaries,  and  indirect  costs  was 
approximately  $31,000. 


« 


« 
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CHAPTER  IV 
DPG  SOFTWARE 


A.  Introduction 

The  theory  of  data  analysis  and  computer  software  for  the 
original  design  of  the  DPG  were  developed  during  the  drafting  and 
machining  stages  of  the  project.  The  complete  theory  and  software  to 
directly  develop  the  first  seven  directional  Fourier  coefficients  for 
a  fully  functioning  4~arm  DPG  are  presented  in  this  chapter.  After 
initial  evaluation  of  the  instrument's  first  sets  of  field  data,  it 
was  considered  that  only  one  differential  pressure  gauge  on  each  axis 
was  reliably  accurate.  Hence  the  software  package,  as  used  at  the 
FRF  installation  site,  analyzes  only  two  orthogonal  differential 
pressure  channels  along  with  the  absolute  signal.  The  FRF  package, 
then,  is  an  abridged  version  of  the  program  described  herein.  It 
develops  only  water  surface  slope  in  the  x  and  y  directions  and 
vertical  displacement,  and  the  first  five  directional  Fourier 


coefficients . 


B.  Theory  of  Analysis 


Selecting  the  x-y  coordinate  system  in  the  horizontal  with 
the  origin  at  the  instrument's  center  as  in  Figure  IV-1,  the  DPG  is 
assumed  to  faithfully  record  the  water  surface  displacement  ~n(x,y,t) 
at  the  center,  the  slopes  r;x  at  A  and  C  and  r)y  at  B  and  D.  The 
curvatures  ^xx  and  T]yy  across  the  center  can  be  ari thmetical ly 
approximated  by  (ilx(A)  ~  1]x(C))/ax  and  (Tly(B)  -  r\y(D))/Ay 
respectively.  With  the  present  DPG  configuration,  it  is  not  possible 
to  arithmetically  create  the  cross  curvature  term  r^Xy.  An  additional 
measurement  would  be  necessary  to  define  such  a  term. 


y 


FIGURE  IV-1:  Coordinate  system  adopted  for  DPG  analysis. 
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The  analysis  procedure  follows  from  Longuet-Higgins , 
Cartwright,  and  Smith  (1963).  The  water  surface  displacement  is 
expressed  in  terms  of  the  complex  directional  amplitude  spectrum, 
F(6,0),  as: 


r)  (x,y ,  t)  = 


r2lT 


J  _-0O 


,  i (at-k  x-k  y ) . .  , 
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frequency 

direction  (assuming  G  measured  counter-clockwise 
from  the  positive  x-axis) 
kcosG,  wave-number  in  the  x-direction 
ksinQ,  wave-number  in  the  y-direction, 

-1. 


The  water  surface  slope  and  curvature  follow  from  (4.1)  as: 
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Information  providing  a  partial  description  of  the 
directional  energy  spectrum  S(6,0),  (equal  to  |F(6,8)|  ),  is 
contained  in  the  auto-  and  cross-  spectra  of  the  water  surface 
displacement,  slope,  and  curvature  terms.  These  spectra  are  obtained 
through  Fast  Fourier  Transform  procedures  using: 
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The  directional  spectrum  at  a  frequency  6  is  then  expressed 
as  a  partial  Fourier  series  in  terms  of  the  wave  direction  0  as: 


N  N 

S(a,0)  =  a  +  £  a  cos  n  0  +  £  sin  n  0 

n=l  n=l 


(4.14) 


The  Fourier  coefficients  are  determined  in  the  usual  manner 
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b  =  [sn  n  (0)] 

2  .2  rj  rj 
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lirk 

The  truncated  Fourier  series  directional  spectrum  is  now 
easily  developed  from  Equation  4.14. 

If  two  measurements  each  are  accurately  made  of  T)x  and  T\y,  it 
would  be  possible  to  generate  all  of  the  first  and  third  coefficients 
four  different  ways  and  the  second  coefficients  two  different  ways 
using  various  combinations  of  the  differential  gauges  in  auto-  and 
cross-spectra.  The  instrument  software,  (re-printed  in  Appendix  F)  , 


generates 

each  of 
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possible 

coefficients . 

There 

are 

two 

advantages 

here.  The 

program  will 

use  the  mean 

value 
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each 

calculated 

angular  Fourier 

coefficient  in  generating 
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spectrum  if  it  is  decided  that  each  transducer  is  working  correctly. 
If  it  is  concluded  that  a  transducer  is  malfunctioning  during  a 
sampling  interval,  the  program  will  disallow  that  transducer's  signal 
during  the  final  coefficient  calculation.  The  redundancy  of  the 
transducers'  measurements,  then,  would  offer  a  potentially  more 
accurate  estimate  of  the  directional  Fourier  coefficients  as  well  as 
a  first-stage  failsafe  mechanism  if  one  of  the  differential 
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transducers  fails.  This  was  indeed  the  case  in  the  field  evaluation. 

C.  Correction  for  Centering  of  Measurements 

The  analysis  theory  used  herein  assumes  that  all  of  the  water 
surface  measurements  are  taken  about  one  point  in  the  generation  of 
the  first  five  directional  Fourier  coefficients.  However,  the 
physical  orientation  of  the  instrument  is  such  that  the  water  surface 
slopes  are  each  taken  about  a  point  midway  along  each  arm 
(approximately  two  feet  from  the  instrument  center  in  the  present 
design).  This  analysis  difficulty  can  be  overcome  by  averaging  the 
two  slope  measurements  along  each  axis  and  using  the  resulting  mean 
slope  values  in  the  analysis.  This  is  done  in  the  present  program  if 
both  transducers  along  an  axis  are  deemed  operational. 

The  measured  slope  values  could  be  transferred  to  the  center 
of  the  instrument  using  a  complex  response  function.  If  one 
considers,  for  example,  the  absolute  pressure  signal  to  be  monitored 
at  the  instrument  center  (x0,y0),  and  an  x~axis  differential  pressure 
signal  measured  at  (xj.y^),  the  water  surface  and  displacement  terms 
could  be  expressed  (from  Equations  4.1  and  4.2)  as: 


r00  r2v 

n(xo'yo'fc)  =  F(o)eX 

J  -co  Jq 


.  (a t-k  x  -k  y  ) 

x  o  y  o  30  3a 


(4.25) 


and  for  the  cross-correlation  of  slope: 


S  (a)  =  S*  (o) 
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(A. 31J 


where  the  starred  spectra  are  those  measured  by  the  instrument. 

The  transfer  function  alters  only  the  phase,  not  the 
amplitude,  and  is  a  function  of  wave  frequency,  direction,  and  gage 
length  of  the  instrument  arms.  The  corrections  require  one  to  know 
the  wave  direction  of  each  frequency  component  a  priori .  These  can 
be  calculated  (if  two  transducers  along  an  axis  are  functioning 
properly)  using  the  simultaneous  solution  of  a  set  of  equations  such 
as : 
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where  (xj.yj),  (xj.yj),  and  (x2,y2),  (x4,y4)  are  the  known  positions 
of  the  x-  and  y-axis  slope  measurements,  respectively. 
Alternatively,  the  analysis  could  be  carried  out  with  the 
non-corrected  spectra,  the  directions  at  each  frequency  estimated, 
and  the  process  repeated  using  spectra  corrected  by  the  previously 
estimated  direction  to  yield  continually  more  accurate  directional 
information.  This  analysis,  however,  can  correct  only  for  the 
principal  direction. 

Such  corrections  are  assumed  small  enough  to  be  neglected  in 
the  present  analysis.  For  the  worst  case  of  small  period  waves, 
(say,  3.2  seconds),  the  phase  shift  is  less  than  1  degree. 

D.  Program  Operation  Notes 

1.  Introduction 

The  FORTRAN  analysis  program,  presented  in  Appendix  F, 
accesses  and  analyzes  the  five  raw  digitized  voltages  recieved  from 
the  DPG  data  storage  tape.  A  crude  flow  chart  of  the  program  is 


shown  in  Figure  IV-2. 
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FIGURE  IV-2:  Flow  chart  of  DPG  analyaii  program 
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2.  Input  and  Screening  of  Data 

A  preliminary  scan  of  the  voltage  records  is  made  tc 
calculate  the  mean  and  standard  deviation  of  each.  In  each  record, 
sampled  points  that  are  greater  or  less  than  3  standard  deviations 
from  the  mean  are  truncated  to  the  average  value  of  their  two 

neighboring  points.  The  number  and  magnitude  of  high  and  low 

truncations  made  in  each  record  are  reported  at  the  terminal.  If  an 
unusually  large  number  of  alterations  are  necessary  for  a  record,  the 
operator  may  choose  to  disallow  the  particular  record  from  further 
spectral  analysis.  The  mean  of  each  record  is  recalculated  and  is 
then  subtracted  from  each  of  the  transducer  signals. 

3.  Conversion  from  Voltage  to  Pressure  Signals 

The  differential  signals  are  converted  from  volts  to  pounds 
per  square  inch  differential  (PSID)  which  corresponds  to 

P (arm) -P (center) .  These  values  are  then  divided  by  the  distance 
between  the  sensor  positions  for  each  arm  to  yield  the  pressure 

differential  per  inch  of  length  along  an  arm.  To  maintain  a 
consistent  sign  convention  among  the  pressure  differences  across  the 
instrument,  the  signs  of  the  dPl  and  dP2  signsals  are  reversed, 
four  differential  channels'  signals  are  now  in  the  form: 

dPl  =  (P(center)  -  P(arml))  /  axI 
dP2  =  (P(center)  -  P(arm2))  /  ax2 
dP3  =  (P(arm3)  -  P(center))  /  ax3 
dPA  =  (P(arm4)  -  P(center))  /  ax4 


The 
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ft 


With  this  convention,  the  program  as  written  will  generate  the 
directional  distribution  with  the  greatest  energy  concentrated  about 
the  compass  heading  to  which  the  waves  are  moving  towards  with  DPG 
arm  1  =0  degrees.  (The  heading  is  later  adjusted  to  represent  the 
compass  direction  from  which  the  waves  propagate  with  respect  to  true 
north. ) 


The  absolute  signal  is  then  reduced  to  pounds  per  square  inch 
absolute  (PSIA) . 

4.  Instrument  Tilt  Approximation 

The  mean  of  each  differential  pressure  record  is  compared  to 
the  zero-wave  condition  tare  as  determined  in  the  laboratory.  If  the 
mean  of  a  differential  record  is  assumed  to  represent  a  zero-wave 
condition  in  the  field,  then  a  discrepancy  between  the  mean  and  the 
laboratory  tare  could  indicate  instrument  tilt  above  the  seabed.  An 
approximation  of  the  tilt  of  each  arm  is  calculated  using  this 
difference  between  the  mean  and  tare  for  each  transducer  record.  If 
the  arm  tilts  about  the  horizontal  an  angle  oc ,  as  shown  in  Figure 
IV-3,  then  a  pressure  difference  between  the  sensors  is  detected. 
This  pressure  differential  is  converted  to  a  pressure  head,  Q,  and 
the  tilt  a  approximated  from  the  arctangent  of  Q/ax.  Specifically, 
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where  aP  *  dP (record  mean)  -  dP (horizontal)  and  ax  =  gage 
length  of  arm.  The  specific  gravity  of  the  fluid  back-filling  the 
sensor  diaphragms, 

>*  j  must  be  dif¬ 
ferent  than  that  of 
the  surrounding  sea¬ 
water,  /.  If  both 
differential  trans¬ 
ducers  along  an  axis 
are  operating  cor¬ 
rectly,  the  angle  of 
tilt  calculated  for 
col  inear  arms  should 
agree  with  opposite 
signs.  The  direction 
of  tilt  is  defined 
after  considering 
the  manner  in  which  the  differential  pressures  are  reported  by  the 
instrument,  (i.e.,  dP  =  P (arm) -P (center) ) .  An  increase  in  the  mean 
value  of  a  differential  transducer  indicates  that  the  end  of  the 

slanting  upwards  (neglecting  any  changes  in 


5L 


FIGURE  IV-3:  Instrument  tilt  above  the  sea¬ 
floor. 


corresponding  arm  is 
temperature) . 
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If  it  is  decided  that  the  instrument  is  stable  over  time, 
(i.e.,  there  is  no  further  settlement),  then  changes  in  the  mean  of 
the  differential  gauge  signals  could  be  used  to  calculate  the 
ambient  seawater  temperature.  This  estimate  could  then  be  used  to 
develop  the  temperature-corrected  specific  gravity  of  the  seawater 
and  back-fill  fluid.  (This  calculation  is  not  included  in  the 
present  software  package.)  Temperature  estimates  could  not  be  made 
if  the  back-fill  fluid  volumes  on  both  sides  of  the  differential 
pressure  transducers  were  equal.  (See  Section  V.C.6.) 


5.  Tide  Calculation 


The  mean  voltage  of  the  absolute  transducer  record  is 
converted  to  pressure  and  is  used  to  calculate  the  tide  level.  The 
total  average  water  column  above  the  sensor,  z5,  can  be  found  from 
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mean  absolute  pressure 
atmospheric  pressure 
specific  weight  of  seawater 
specific  weight  of  back-filling  fluid 
height  of  sensor  above  the  transducer 
back-fill  pressure. 


Both  the  specific  fluid  weights  should  be  temperature  corrected.  The 
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difference  between  zs  and  the  measured  water  column  abcve  the  sensor 
at  mean  water  represents  the  average  tide  level  during  the  record. 

6.  Development  of  Curvature  Terms 

As  the  next  step,  the  curvature  terms  are  created  by 
subtracting  the  colinear  slope  terms  and  dividing  by  the  distance 
between  the  center  of  the  slope  measurements. 

7.  Significance  of  the  Signals  in  Linear  Theory 

The  absolute  pressure  transducer  signal 
considered  to  be  in  the  form: 

P(x,y,t)  =  y  ^  K  cos  (k  x  +  k  y  -  ot) 

and  the  differential  signals  in  the  form: 

P  (x,y,t)  -  -  y  ^  k  K  sin  (k  x  +  k  y  -  at)  (psid)  (A. 37) 

A  A  X  p  A  y 

P  (x,y,t)  =  -  y  T  k  K  sin  (kvx  +  k„y  "  at)  (paid)  (A. 38) 

y  ^  y  P  x  y  ; 


can  then  be 


(psia)  (A. 36) 
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The  wavenumber,  k,  and  the  pressure  response  function  are 
generated  for  the  frequencies 

2tt 

°(n)  =  NAt  (n-1)  ;  n  =  2'3'---N/2  (A. 39) 

where  N  is  the  number  of  sample  points  analyzed  and  At  is  the 
sampling  interval.  The  water  depths  used  in  the  calculations  include 
the  tidal  variation  computed  earlier. 

8.  FFT  and  Conversion  to  Water  Surface  Terms 

The  absolute  record,  the  four  differential  records,  and  the 
two  curvature  records  are  processed  by  the  Fast  Fourier  Transform 
(FFT)  and  the  resulting  Fourier  coefficients  divided  by  the  pressure 
response  function  and  specific  weight  of  the  seawater.  Frequency 
bands  greater  than  2  rad/sec  (wave  periods  less  than  3.1A  seconds) 
are  disregarded  because  of  the  error  inherent  in  dividing  these  high 
frequency  bands  by  their  respectively  small  K  values. 

The  absolute  record  is  then  comparable  to: 

n (x,y ,t)  =  j  cos  (kxx  +  kyy  -  at)  (in.)  (4. AO) 


and  the  differential  records  to: 
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n  (x,y,t)  =  -  ij-  k  cos  0  sin  {k  x  +  k  y  -  at)  (in. /in.)  (A. 41) 
x  2  X  y 


and 

n  (x,y ,t)  =  -  k  sin  0  sin  (k  x  -  k  y  -  at)  (in. /in.)  (4.42) 

y  ^  *  y 

while  the  curvature,  or  acceleration,  terns  are  as: 

n  (x,y,t)  =  -  ^  k2  cos^0  cos(k  x  +  k  y  -  at)  (inf1)  (4.43) 

xx  J  2  x  y 


and 


nyy(x,y,t)  = 


H  2  2 

-  —  k  sin  cos(kxx  -  k^y  -  at) 


(in.-1)  (4.44) 


Twice  the  sun  of  the  squares  of  the  modified  Fourier 
coefficients  represents  the  power  spectrum  of  the  displacement,  slope 
components,  and  curvature  components.  The  significant  wave  height  is 
determined  from  four  times  the  square  root  of  the  sum  of  the  absolute 
channel's  power  spectrum: 


(4.45) 
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9-  Failure  of  the  Absolute  Channel 

If  it  is  decided  that  the  absolute  transducer  has  failed,  the 
water  surface  displacement  and  power  spectrum  can  theoretically  be 
estimated  from  the  x~  and  y-axis  differential  transducers  that  are 
deemed  as  operational.  This  procedure  necessarily  assumes  that  there 
is  only  one  wave  direction  per  frequency.  The  record  of  the  absolute 
transducer  has  been  developed  in  the  time  domain  to  that  of  Equation 
(4.  AO): 

n(x,y,t)  =  —  cosik^x  +  k^y  -  ot) 

so  that  the  power  spectrum  of  the  absolute  transducer  is  as: 

2  -  2 

Snr,(o)  =  T  cos2<kxx  +  V  "  at)  "  T”  (A. 46) 


and  the  power  spectrum  of  the  x  and  y  differential  gauges  are 
correspondingly: 
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The  power  spectrum  of  the  absolute  record,  then,  mifeht  be  developed 
from  the  sum  of  the  x  and  y  differential  power  spectra  divided  by  the 
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10.  Failure  of  3  Differential  Channels 


If  only  the  absolute  and  one  differential  pressure  channels 
are  considered  operational,  another  estimate  of  the  wave  direction 
might  be  made  with  a  180  degree  ambiguity.  With  the  power  spectra  of 
the  absolute  and  differential  channels  in  the  form  of  Equations  4.46 
through  4.48; 
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where  SXx  (6)  is  the  auto-spectrum  of  a  transducer  signal.  This  is 
an  estimate  of  the  "effective  wave  period"  during  the  record  as 
determined  by  each  transducer. 


12.  The  Directional  Fourier  Coefficients 

The  directional  (or  angular)  Fourier  coefficients  are 
calculated  for  each  frequency  band  using  Equations  4.18  through  4.24. 
The  coefficients  are  taken  as  the  mean  of  all  those  calculated  with 
different  combinations  (or  cross-spectra)  of  operating  channels. 
These  angular  coefficients,  as  well  as  the  water  surface  displacement 
power  spectrum,  are  block-averaged  over  a  specified  number  of 
adjacent  frequency  bands.  Since  an  ocean  wave  time  series  is 
composed  of  random  functions,  the  Fourier  coefficients  of  ea<'K 
channel  must  not  be  block-averaged  before  generation  of  the  angular 
coefficients.  The  average  of  a  random  function  exhibits  zero 
correlation  with  the  average  of  another,  so  that  the  angular 
coefficients  must  be  calculated  from  the  raw  cross-spectra  terms  and 


then  averaged.  Similarly,  the  directional  estimates  obtained  from 
Equations  A. 52,  A. 53,  and  A.5A  are  block-averaged. 

13.  Block-Averaging 

The  stability  of  the  generated  spectra  is  enhanced  by 
block-averaging  over  adjacent  frequency  bands.  If  the  water  surface 
displacements  are  assumed  to  be  Gaussian,  the  Fourier  coefficient 
spectral  estimates  are  distributed  according  to  the  chi-square 
distribution  with  two  degrees  of  freedom.  The  block-averaged,  or 
smoothed,  Fourier  coefficient  spectral  estimates  are  represented  by 
the  chi-square  distribution  with  the  number  of  degrees  of  freedom 
equal  to  two  times  the  number  of  adjacent  frequencies  averaged 
(Borgman, 1972) .  When  the  analysis  program  operates  upon  records  from 
computer  simulation  or  a  controlled  laboratory  environment,  each 
spectral  estimate  is  considered  separately;  that  is,  there  is  no 
block-averaging.  This  is  equivalent  to  two  degrees  of  freedom  and  is 
done  to  retain  the  highest  resolution  possible  for  the  narrow  spectra 
and  shorter  time  records  generated  in  simulation  and  the  laboratory. 
When  the  program  operates  upon  a  record  from  the  ocean  environment, 
each  signal  is  block-averaged  over  four  to  eight  adjacent  spectral 
frequencies  to  yield  eight  to  sixteen  degrees  of  freedom 
respectively.  It  is  assumed  that  the  frequency  resolution  thereby 
lost  in  the  block-averaged  analysis  of  the  ocean's  broader  spectra  is 
balanced  by  the  increased  stability  of  the  record. 
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14.  Frequency  Bands  of  Greatest  Energy 

The  block-averaged  frequency  bands  of  maximum  energy  are  then 
identified  from  one  of  the  operative  channels  (generally  the 
absolute)  and  the  amount  by  which  each  exceeds  the  average  energy  in 
the  spectrum  is  calculated.  The  truncated  Fourier  series  directional 
spectrum  (Equation  4.14)  is  then  calculated  for  each  of  these  high 
energy  bands.  The  direction  of  peak  energy  for  each  spectrum  is 
reported  as  that  band's  direction  of  wave  propagation  —  adjusted 
accordingly  to  instrument  orientation  relative  to  true  north. 

The  "mean  direction"  of  each  band  can  also  be  calculated  as: 


n  0 

n 


tan 


-1 


b 

n 

a 

n 


(4.56) 


where  n»l,2,3.  (If  no  water  surface  curvature  terms  are  developed, 
then  n*l,2.)  There  exists  a  directional  ambiguity  for  n>l  since  there 
exists  more  than  one  peak  in  sin  n9  for  n>l  (from  which  the  bn  terms 
are  developed). 


Since 

the 

ultimate  intent  of  the 

DPG  as  a  tool  is 

to  report 

predominant 

wave 

characteristics  at 

the  installation 

site. 

the 

presentation 

of  wave  information  at  the 

dominant  energy 

levels 

is 

considered  to  be  the  preferred  data  display. 
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E.  DPG  Simulation  and  Error 

A  simple  simulation  program  was  utilized  to  test  the 
capabilities  of  the  software.  A  directional  wave-pressure  field  over 
water  depth  h,  at  transducer  depth  s,  was  considered  using 


P(x,y,s;t) 


J  H  . 
Yh  +  l  Y 

j=l 


K  cos  (k  x+k  y-ot) . 
p.  x  y  j 

D 


(A. 57) 


where  J  *  number  of  wave  trains  present.  Multiple  wave  trains  of 
different  frequencies  and/or  directions  were  simulated  by 
superimposing  all  of  the  dynamic  water  surface  elevations  at  (x,y,t) 
onto  the  static  water  column,  h.  The  pressure  at  those  (x,y)  points 
corresponding  to  the  five  DPG  diaphragm  locations  were  generated  over 
time.  (It  was  assumed  that  the  five  diaphragms  in  the  center  of  the 
instrument  occupy  the  same  point.)  From  these  values,  four 
differential  signals  and  one  absolute  signal  were  developed  at  At 
intervals  for  a  record  length  of  Nat  to  simulate  DPG  in-situ 
performance.  Non-linear  wave-wave  interactions  were  not  considered 
in  this  simple  model. 

The  integrity  of  the  FFT  analysis  was  first  tested  by 
simulating  a  single  wave  train  at  a  principal  frequency  of  the  FFT. 
Results  of  the  analysis  of  the  simulated  record  indicated  acceptable 
FFT  performance  since  all  of  the  energy  was  indeed  found  to  be 
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concentrated  at  the  simulated  wave  frequency. 

The  ability  to  discriminate  a  number  of  different  wave  trains 
was  next  tested.  Two  wave  trains  of  frequencies  within  one  spectral 
bandwidth  of  one  another  and  of  different  energy  and  direction  were 
superimposed  upon  a  train  of  lower  frequency.  Analysis  was  carried 


TIME  (seconds) 


FIGURE  IV-5:  Simulated  time  record  of  water  surface 

elevation.  Case  1.  H  *  4  ft.,  0  *  46  deg.,  5-98 
sec.  pd. ;  11*3  ft.,  6  “  180  deg.,  6.33  sec.  pd.; 

H  ■  4  ft.,  0  “  300  deg.,  9.00  sec.  pd. 


out  for  512  sample  points  of  0.25  second  interval.  A  portion  of  the 
simulated  time  record  of  the  water  surface  displacement  generated 
over  each  of  the  five  sampling  points  of  the  DPG  is  shown  in  Figure 
IV-5.  The  raw  spectral  energy  distribution,  or  energy  spectra, 
presented  in  Figure  IV-6,  illustrates  the  presence  of  three  distinct 
frequencies  well.  Although  the  simulated  lowest  and  highest 
frequency  waves  were  of  the  same  wave  height,  the  higher  frequency 
wave  demonstrates  a  smaller  energy  spike  due  to  the  spectral  leakage 
between  the  lower  frequency  bands.  The  non-smoothed  7-coef f icient 
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FIGURE  IV-7:  Non-weighted  7  coefficient  Fourier  series 

directional  energy  distribution.  Simulation  Case  1. 
Illustrates  presence  of  three  distinct  wave  trains  with 
different  directions. 
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Fourier  series  directional  distribution  was  plotted  over  two  degree 
increments  £or  the  three  bands  of  highest  energy  and  is  presented  in 
Figure  IV~7.  The  analysis  program  demonstrates  its  ability  to  define 
each  wave  train's  direction  accurately  by  presenting  the  maximum 
energy  at  the  correct  simulated  wave  direction.  The  half-power 
directional  width  of  each  is  approximately  60  degrees.  These  curves, 
of  course,  do  not  contain  all  of  the  simulated^ave  energy  since  only 
three  major  spectral  bands  are  plotted. 

Two  wave  trains  of  different  frequency  but  with  identical 
height  and  direction  were  next  simulated  over  512  points  of  0.25 
second  interval.  The  frequency  of  each  was  selected  to  fall  exactly 
between  two  spectral  frequencies.  The  higher  frequency  bands  were 
0.1184  rad/sec  apart  and  the  lower  were  0.464  rad/sec  apart.  The  raw 
energy  spectra,  shown  in  Figure  IV~8,  depicts  the  two  wave  trains' 
energy  distribution  over  frequency.  Slight  differences  in  leakage 
between  the  two  simulated  wave  frequencies  can  be  detected.  The 
non~smoothed  directional  distributions  of  the  two  waves,  shown  in 
Figure  IV-9,  agree  well  with  each  other  as  expected.  The  very  slight 
differences  are  the  result  of  differences  in  spectral  leakage  for 
each  wave. 

The  analysis  procedure,  however,  is  unable  to  differentiate 
between  two  waves  of  identical  frequency  and  different  direction. 
Two  waves  of  the  same  height  and  non-spectral  frequency  were 
simulated  with  a  130  degree  difference  in  direction.  The 
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FIGURE  IV~8:  Raw  energy  spectra;  Simulation  Case  2. 


ENERGY-SPECTRA  of  DISPLACEMENT  TERM 


FIGURE  IV~9:  Non-weighted  7-coefficient  Fourier  eerie*  dir¬ 
ectional  energy  distribution;  Simulation  Case  2. 
Illustrates  good  directional  spectra  agreement  with 


i  slight  differences  due  to  leakage. 


f' 
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non-smoothed  seven-coefficient  Fourier  series  directional 
distribution,  presented  in  Figure  IV-10,  demonstrates  a  wide  peak  of 
80  degree  half-power  bandwidth.  The  distribution  suggests 
omnidirectional  wave  power,  but  with  a  maximum  at  the  mean  of  the  two 
simulated  direction,  80  and  195  degrees. 

The  identical  case  was  re~simulated  with  the  80-degree  wave 
retaining  its  three  foot  height  and  the  195-degree  wave  assuming  a 
four  foot  height.  The  directional  distribution,  shown  in  Figure 
IV-11,  is  asymmetrical  with  peak  energy  at  the  weighted  average  value 
of  the  two  wave  directions. 

Inclusion  of  higher  order  terms  in  the  partial  Fourier  series 
representation  of  direction  should  increase  the  accuracy  of  the 
series  representation.  Figure  IV-12  illustrates  the  heightened 
directional  resolution  of  a  simulated  wave  train  at  a  principal 
frequency  using  the  f^.-st  seven  directional  Fourier  coefficients, 
(N*3)  ,  compared  t<?  the  first  five,  (N“2)  .  The  non-smoothed 

half-power  bandwidth  for  seven  coefficients  is  66  compared  to  76.3 
for  five. 


In  the  simple  simulation 
non-sensitive  to  wave  direction 
one  wave  train.  Four  waves  of 
different  direction,  were  each 
The  directional  distribution  for 


analysis,  the  software  appears  to  be 
—  at  least  for  the  presence  of  only 
equal  energy  and  frequency,  but 
simulated  individually  over  the  DFG. 
each  was  overlaid  upon  one  another 
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and  excellent  agreement  readily  seen. 


FIGURE  IV— 12:  Non-weighted  directional  spectra  using  the 
first  seven  directional  Fourier  coefficients  compared  to 
the  first  five.  Single  simulated  wave  of  spectral 
frequency. 
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CHAPTER  V 

LABORATORY  DEVELOPMENT,  CALIBRATION,  AND  TESTING 
OF  THE  INSTRUMENT 


A.  Introduction 

Some  modifications  to  the  instrument  were  made  during  the 
fluid  back-filling  of  the  sensors.  These  changes,  as  well  as  the 
static  bench  calibration  of  the  DPG  instrument  are  documented  and 
discussed  in  this  chapter.  The  preliminary  wave-tank  tests  carried 
out  at  the  Coastal  Engineering  Research  Center  (CERC) ,  Fort  Belvoir, 
Virginia,  are  also  described. 


B.  Fluid  Back-Filling  the  Sensors 
System  Modifications 

Certainly  the  most  difficult  problem  encountered  during  the 
assembly  of  the  DPG  was  the  fluid  back-filling  of  the  sensing 
diaphragms.  The  frustration  and  numerous  failures  during  the  very 
intense  six-week  period  of  back-filling  attempts  led  to  two  changes 
each  of  the  type  of  fluid  and  tubing  used. 
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The  original  design  called  for  1/16  inch  inner  diameter 
stainless  steel  tubing  between  the  water-tight  cylinder  and  isolation 
diaphragms.  The  fuselage,  or  center,  tubings  were  14  inches  in 
length,  while  the  arm  tubings  were  44  inches  in  length.  Methanol  was 
first  tried  as  the  back-filling  fluid,  but  then  discarded  because  of 
the  poisonous  hazard  it  presented.  Standard  automobile  transmission 
fluid  was  next  used  as  the  back-filling  agent.  Although  much  safer, 
it  easily  contaminates  laboratory  tools  and  work  areas. 

Hundreds  of  attempts  were  made  to  back-fill  the  sensors  and 
eliminate  all  air  in  the  system.  An  air  bubble  of  appreciable  size 
could  compress  under  a  dynamic  load  and  thereby  prevent  the 
transducer  from  sensing  the  full  pressure  exerted  upon  the  isolation 
diaphragm.  The  back-filling  fluid  itself  is  assumed  incompressible 
so  that  with  no  air  in  the  system  the  fluid  transmits  the  pressure 
sensed  by  the  isolation  diaphragm  immediately  to  the  transducer.  The 
final  technique  used  to  back-fill  the  DPG  is  outlined  in  Appendix  D. 

Once  satisfied  that  air  was  removed  from  the  system,  more 
problems  arose  ■  when  the  transducers,  sensing  through  the  1/16  inch 
I.D.  tubing  filled  with  transmission  fluid,  failed  to  respond  to  a 
load  faster  than  five  to  twelve  seconds.  It  was  decided  that  the 
frictional  head  loss  caused  by  the  high  viscosity  of  the  fluid  and 
long,  small  diameter  tubing  might  be  responsible  for  such  large 
response  times.  It  was  also  considered  possible  that  small  bubbles 
might  still  be  lodged  in  the  system  or  that  impurities  in  the  fluid 


may  have  blocked  the  transducer  filter  screens  and  diaphragm. 


The  frictional  pressure  gradient  along  the  inside  length  1  of 
a  circular  pipe  is  (Peerless , 1967)  : 


AP 


£  128  jiQ 

4 

ttD 


(5.1) 


where  aPj  “  frictional  pressure  drop 

p  ■  coefficient  of  dynamic  viscosity 
Q  •  volumetric  flow  rate 
□  “  inner  diameter  of  tubing. 


Figure  V-l  presents  a  schematic  representation  of  the  differential 
transducer  and  sensors.  The  fluid  flow  along  each  arm,  Q,  and  Q2 , 
that  results  from  difference  in  pressure  between  the  system's  sensors 
may  be  expressed  from  Eq.  5.1  as: 


£>1  =  B1(P1  -  P  >  (5.2a) 

^2  ~  B2^Pc2  ~  P2^  (5.2b) 


where 
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(5.3) 
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Re-writ ing; 

pi  -  %  '  «i/Bi 


And,  from  continuity; 


Ql  =  Q2  =  Q 


Adding  Equations  5.4a  and  5.4b; 


P 


1 


The  transducer  diaphragm  displaces  a  volume  V-c  proportional 
change  in  pressure  across  the  transducer  chamber: 


AV  =  aAP  = 
c  c 


) 


Re-writing  Equation  5.7  and  expressing  the  change  in  volume  in 
of  the  fluid  flow  Q; 


P  =  —  AV 
c2  a  c 


’  W  I  8  « 


(5.4a) 


(5.4b) 


(5.5) 


(5.6) 

to  the 


(5.7) 

terms 


(5.8) 
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Differentiating  both  sides', 


Q  =  a  -  (APc) 


(5.9) 


Substituting  (Equations  5.8  and  5.9  into  5.6,  one  obtains: 


aP  -  P1  -  +  iPc  +  +  5^>  h  APc 


(5.10) 


where  aP  is  the  difference  in  pressure  between  the  two  sensors.  If 
the  change,  aP,  occurs  instantaneously,  it  can  be  shown  that  the 
solution  to  Eq.  5.10  is 


AP  =  Ap(l  -  e'Kot) 
c 


(5.11) 


substitution  into  Eq.  5.10  yields: 


AP  =  AP  +  a(~~  +  K  AP  e  KDt 
c  B1  B2  D 


(5.12) 


where  K0_1  is  the  "folding  time"  of  the  response.  Solving, 


a(B1  +  b2) 
B1  B2 


128  a  v1  v2  ^  &2 

4  _  4 
tt  Dx  D2 


Vl 


U2*2 


(5.13) 


If  the  fluid  and  tubing  diameter  are  the  same  on  both  sides  of  the 
transducer,  Eq.  5.13  reduces  to: 


128  a  p  l2 


IT  D 


(5.14) 


The  form  of  Eq.  5.11  was  verified  qualitatively  by  examining 


the  strip  chart  output  of  a  differential  transducer  and  sensor  system 
loaded  in  an  approximately  instantaneous  manner. 
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It  is  clear  that  the  response  time  is  proportional  to  the 
fluid  viscosity  and  inversely  proportional  to  the  fourth  power  of  the 
tubing  diameter.  Changing  the  tubing  length  was  not  considered  since 
the  lengths  were  pre-determined  by  the  DPG  configuration  as  descirbed 
in  Chapter  II. 

It  was  first  thought  that  the  response  time  could  be 
sufficiently  decreased  by  reducing  the  fluid  viscosity.  A  number  of 
fluids  were  considered  but  all  discarded  because  of  their  degrading 
effects  on  the  system's  materials,  (i.e.,  acrylic,  neoprene,  Buna-n, 
nylon,  and  stainless  steel).  After  two  hundered  and  eighty  hours  of 
nearly  continuous  laboratory  work,  it  was  decided  to  use  gin  as  the 
back-filling  fluid.  It  is  inexpensive,  easily  obtainable,  safe,  and 
has  a  dynamic  viscosity  coefficient  one  and  one  half  orders  of 
magnitude  less  than  transmission  fluid.  The  transducers  and  sensors 
were  drained  of  transmission  fluid  and  flushed  briefly  with  hexane, 
then  methanol,  then  back-filled  with  gin.  A  50-50  combination  of  40 
proof  McCalls  and  80  proof  Gilbeys  gin  was  used.  The  water  content 
of  the  gin  required  that  the  fluid  rest  in  the  system  for  twenty-four 
hours  during  rising  temperatures  so  that  air  released  from  the  water 
could  be  collected  and  removed. 

Although  the  system  response  time  improved,  it  remained 
unacceptably  high.  It  is  crucial  for  each  of  the  five  DPG 
transducers  to  respond  quickly  and  simultaneously  to  ensure  proper 
instantaneous  measurement  of  the  sub-surface  dynamic  pressure. 
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C.  Bench  Calibration 

1.  The  Calibration  Technique 

An  extremely  simple  calibration  was  first  considered  to 
statically  determine  the  static  response  of  the  transducers  to  a 
known  load.  The  isolation  diaphragms  were  positioned  facing  upwards 
and  were  loaded  with  known  weights.  This  technique  was  ineffective 
because  of  the 
physical  nature 
of  the  dia¬ 
phragms  .  The 
bottoms  of  the 
weights  were 
rarely  in  full 
contact  with 
the  elastomer, 
as  is  depicted 
graphically  in 
Figure  V-2,  so 
that  the  ap¬ 
plied  pressure 
to  the  sensor 
could  not  be 

easily  calculated.  If  the  weight  was  too  heavy  and  large,  some  of 


FIGURE  V-2:  Cross-section  of  isolation  dia¬ 
phragm  with  an  object  load  applied. 


the  force  of  the  weight  would  be  exerted  upon  the  top  face  of  the 
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acrylic  housing  and  not  upon  the  fluid.  (intuitively,  one  would 
expect  this  to  happen  even  in  the  ideal  case  of  force  applied  by  the 
dynamic  water  column  —  if  the  exposed  elastomer  area  was  so  small 
that  it  could  not  deflect.) 

A  slightly  more  sophisticated  calibration  system  was  then 
designed.  The  system,  as  drawn  in  Figure  V-3,  consists  of  two 
pressure  chambers  which  secure  around  each  isolation  diaphragm  of  a 
differential  transducer,  (only  one  is  used  with  the  absolute 
transducer).  Each  chamber  is  equipped  with  its  own  pressure 
transducer  and  a  bicycle  tire  valve  stem.  The  chambers  are  connected 
by  tubing  with  a  valve  in  the  center  to  isolate,  bleed,  or  allow  the 
chambers  to  communicate  with  one  another.  The  chambers  are  built  of 
short  steel  cylinders  with  a  plate  welded  on  one  end.  Each  chamber 
is  secured  to  a  sensor  diaphragm  with  a  short  length  of  motorcycle 
tire  inner  tube  that  is  glued  to  the  chamber  and  hose-clamped  to  the 
isolation  diaphragm  housing.  The  system  is  pressured  with  a  bicycle 
tire  pump  and  the  pressure  sensed  in  each  chamber  and  by  the  DPG 
transducer  is  reported  on  a  multi-channel  strip  chart.  The  DPG 
transducers'  response  was  interfaced  with  the  strip  chart  using  a 
short  "bench  cable"  that  mated  with  the  instrumentation  connector 
located  within  the  DPG  transducer  stack. 
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time,  while  the  other  side  was  assumed  to  remain  at  atmospheric 
pressure.  The  calibration  technique,  then,  was  not  able  to  simulate 
differential  pressure  loading  while  both  sides  were  subjected  to  a 
pressure  head  as  would  be  felt  on  the  seafloor.  (Both  chambers 
leaked  somewhat,  so  that  it  was  not  feasible  to  pressure  each  side  to 
a  particular  head  and  then  begin  the  differential  loading.) 

The  calibration  technique  described  herein  entailed  static 
loading  whereas  the  operating  instrument  would  be  subjected  to 
dynamic  loading.  This  was  not  of  great  concern,  however,  since  the 
in~si tu  loading  is  of  relatively  low  frequency.  The  calibration 
technique  could  be  improved  by  using  an  oscillating  pump  to  pressure 
the  chambers  sinusoidally  over  time. 

2.  Testing  for  Air  Bubbles 

Theoretically  the  calibration  system  can  be  used  to  check  for 
the  presence  of  air  bubbles  in  the  transducer  sensor  lines.  Firstly, 
air  in  a  line  will  compress  under  a  load  and  bias  the  pressure  value 
that  the  transducer  will  report.  Secondly,  it  is  hypothesized  that 
the  loading  and  unloading  of  the  air  bubble  will  delay  the  response 
time  of  the  system.  If  this  is  true,  then  it  may  be  possible  to 
detect  an  air  bubble  by  loading  each  side  of  a  differential 
transducer  equally,  then  unloading  them  simultaneously.  The  side 
with  the  greater  air  volume  will  lag  in  response.  (It  is  assumed 
that  the  inherent  lag  in  arm-side  response  due  to  the  longer  arm 


tubing  length  is  negligible  for  large  tubing  diameters  (  >  1/8  inch 


inner  diameter).  The  length  of  small  diameter  nylon  tubing  is  the 
same  for  each  side  of  a  transducer.)  This  effect  was  never  tested 
exhaustively  in  the  laboratory,  but  one  possible  such  occurence  is 
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FIGURE  V-4:  Strip  chart  recording  illustrating  possible 

existence  of  air  in  arm  side  of  dP4.  A)  arm  sensor 
loaded;  B)  arm  and  center  sensors  communicate  to 
establish  zero  differential  pressure;  C)  both  sensors 
bled  simultaneously,  arm  side  shows  slight  lag.  Chart 
speed:  5  mm/sec. 


shown  in  Figure  V-4.  The  chamber  on  the  arm  sensor  of  dP4  was  loaded 

(A)  and  then  allowed  to  communicate  with  the  center  sensor  chamber 

(B)  .  The  differential  transducer  correctly  reported  zero  difference 
in  pressure.  When  the  entire  system  was  bled  (C) ,  the  transducer 
briefly  indicated  a  greater  pressure  on  the  arm  side.  The  arm  was 
re-back-filled  to  be  safe. 

3.  Static  Calibration  Results 

The  calibration  results  for  each  of  the  four  differential 
pressure  and  the  absolute  pressure  transducers  as  found  during  the 
pre-FRF  calibration  are  shown  in  Figures  V-5  (a)  through  (e) . 
Differential  transducers  dP2,  dP3,  dP4,  and  the  absolute  transducer 
each  demonstrated  acceptable  one-to-one  reporting  of  the  loads 
applied.  DPI,  however,  is  significantly  biased.  The  problem  was 
attributed  to  the  transducer  itself,  since  a  similar  error  appeared 
during  the  first  calibration.  The  zero-differential  value  of  dPl  was 
also  consistently  high  (5  to  7  volts  as  compared  to  2.5  volts,  as 
ordered);  and  unlike  the  other  transducers,  it  demonstrated 
sensitivity  to  excitation  voltage  as  is  discussed  in  Section  V.C.4. 
This  malfunction  did  not  appear  until  just  before  the  scheduled 
testing  was  to  begin  at  CERC  and  time  limitations  did  not  allow  for 
its  immediate  repair.  The  calibration  results  of  the  first 
investigation  for  the  absolute  transducer  and  dP2,  dP3,  and  dP4  are 
identical  to  those  of  the  second. 
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The  calibration  curves  were  generated  from  the  strip  chart 
output  reproduced  in  Appendix  B.  The  calibration  chamber  was 
pressured.  allowed  to  leak  slightly,  and  the  chamber  and  DPG 
transducer  readings  taken  from  the  strip  chart  at  the  point  before 
the  system  was  loaded  again.  The  response  time  for  the  DPG 
transducers  was  determined  by  measuring  the  passage  of  time  between 
loading  events  as  reported  by  the  chamber  transducer  and  the  DPG 
transducer  on  the  strip  chart.  Each  transducer  appears  to  respond  in 
an  acceptably  similar  short  time  period. 

A.  Power  Supply  Sensitivity 

Each  of  the  transducers  demonstrated  satisfactory 
insensitivity  to  power  supply  voltage  fluctuations  except  dPl. 
Whereas  the  absolute  and  dP2,  dP3,  and  dPA  transducers  were 
unaffected  by  voltage  changes  from  20  to  26  volts,  the  output  of  dPl 
decreased  in  an  unexplained  manner  at  excitation  voltages  below  22 
volts . 

5.  Transducer  Drift 

The  absolute  and  dPl,  dP2,  and  dP3  transducers  proved  steady 
in  output  with  steady  loading.  However,  dPA  intermittently  and 
erratically  drifted  during  the  second  set  of  bench'  calibration  tests. 
It  is  assumed  that  such  drift  will  be  simply  interpreted  as  a  long 
period  wave  by  the  FFT  analysis  of  the  data.  But  the  inability  of 


the  transducer 


to  establish 


a  zero-differential  tare  value 


invalidates  any  calculations  of  instrument  tilt  or  ambent  temperature 
using  dPA  signals  (Section  IV. C. A). 

6.  Temperature  Effect 

All  transducers  are  temperature  sensitive  to  some  degree. 
The  thermal  sensitivity  specifications  of  the  DPG's  Setra  absolute 
and  differential  pressure  transducers  are  less  than  +0.022  FS/°F  and 
+0.032  FS/°F  thermal  zero  shift  with  thermal  coefficients  of 

sensitivity  of  less  than  +0.0152  FS/°F  and  +0.022  FS/°F  respectively. 
The  transducer  and  sensor  systems  in  the  DPG  are  slightly  more  or 
less  sensitive  than  the  factory  specifications  because  of  their  fluid 
back-fill.  As  the  DPG  environment  changes  temperature  the 
back-filling  gin  expands  or  contracts.  Moreover,  a  temperature 

change  induces  a  slight  differential  pressure  across  the  differential 
transducers  because  of  the  greater  volume  of  fluid  in  the  (longer) 
arm  sensors.  Since  the  volumetric  expansion  of  a  fluid  is  linearly 
proportional  to  the  change  in  temperature  and  initial  volume  of  fluid, 
the  expanded  volume  in  the  arm  sensors  will  be  four  times  that  of  the 

center  sensors,  (since  the  arm  tubing  is  four  times  the  length  of  the 

center  tubing.)  But  since  temperature  changes  under  twenty  feet  of 
ocean  water  are  typically  small  and  of  low  frequency,  it  is  assumed 
that  any  thermal  transducer  drift  during  a  sampling  interval  will  be 
interpreted  as  a  very  long  period  wave  or  as  a  mean  which  changes 


from  record  to  record.  Such  thermal  effects  accordingly  affect  the 
mean  of  a  pressure  record.  This  could  be  mistakingly  interpreted  as 
instrument  tilt  if  two  differential  records  taken  at  different 
ambient  temperatures  are  compared  (Section  III.C.4).  It  could  also 
affect  the  tide  calculation  by  altering  the  back-fill  pressure  term, 
Pbf  •  of  the  absolute  sensor  in  Equation  4.35.  Future  development  of 
the  DPG  should  more  thoroughly  investigate  the  effects  of  temperature 
on  the  transducers  and  sensors. 

7.  Absolute  Sensor  Back-Fill  Pressure 

The  absolute  pressure  sensor  back-fill  pressure,  Pbf  ,  is 
obtained  in  the  laboratory  by  either  of  two  ways.  (1)  The  instrument 
is  positioned  such  that  there  is  no  fluid  head  above  the  transducer 
and  the  absolute  pressure  reported  above  the  atmospheric  is  recorded. 
(2)  The  instrument  is  oriented  upright  and  the  absolute  pressure  that 
is  reported  beyond  the  atmospheric  and  sensor  fluid  head  is  recorded. 

D.  CERC  Wave  Tank  Tests 

1.  Introduction 

The  first  DPG  prototype  was  laboratory  tested  at  the  large 
outdoor  wave  tank  at  the  Coastal  Engineering  Research  Center,  Fort 
Belvoir,  Va.  The  instrument  and  cradle  proved  themselves  easily 
maneuverable  —  both  above  and  below  water.  Although  three  of  the 
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four  differential  pressure  channels  failed  and/or  were  distrurbed  by 
the  wave  tank,  the  analysis  of  measurements  from  the  absolute  and  one 
differential  pressure  channel  suggest  satisfactory  performance  of  the 
DPG. 

2.  Installation 

The  cradle  was  secured  to  a  small  flatbed  trailer  and  the 
instrument  mounted  upon  it.  The  trailer  was  easily  pulled  and 
delivered  to  the  CERC  wave  tank  by  a  light  vehicle.  The  outdoor  wave 
tank  is  a  long,  15  foot  wide,  20  foot  deep  facility  with  a  large 
piston-type  wavemaker.  The  cradle  was  lowered  on  its  side  over  the 


FIGURE  V-6:  Orientation  of  the  DPG  during  CERC  wave  tank  tests. 
Configuration!:  0-55°;  Configuration  II:  0*35°. 


tank  wall  by  ropes  and  dropped  slowly  into  the  water.  It  was 
positioned  by  two  divers  under  10.2  feet  of  water  at  approximately 


mid-length  of  the  tank.  The  cradle,  in  the  center  of  the  tank,  was 
then  secured  by  ropes  to  large  U-bolts  on  return-flow  pipes  that 
terminated  at  this  location,  (Figure  V-6) .  It  was  not  clear  whether 
the  outflow  of  the  pipes  during  testing  would  influence  the  DPG. 

The  instrument,  carried  upon  the  shoulder  of  one  deck-hand, 
was  then  positioned  on  a  trestle  over  the  tank  and  lowered  into  the 
water  by  another  rope  slung  through  an  arm  on  the  trestle.  Care  was 
taken  to  keep  the  instrument  arms  away  from  the  trestle  and  tank 
walls  using  two  lateral  control  lines  that  were  held  on  each  side  of 
the  tank.  The  instrument  was  received  by  a  diver,  guided  and 
positioned  into  the  cradle,  and  secured  with  cable  ties.  The 
orientation  of  the  instrument  was  checked  by  measuring  the  distance 
from  the  wall  to  each  of  the  cradle  legs  with  a  tape  secured  along  a 
wooden  "T."  The  top  of  the  "T"  was  held  against  the  wall  to  ensure 
that  the  measurement  was  taken  perpendicular  to  the  wall.  A  100  foot 
length  of  unarmored  cable,  earlier  mated  to  the  instrument  above 
water,  was  interfaced  with  a  power  supply  and  strip  chart  recorder  in 
a  trailer  alongside  the  tank.  Despite  the  poor  underwater 
visibility,  (less  than  six  inches),  the  entire  installation  took  less 
than  two  hours  and  required  only  two  divers  and  one  deck-hand. 

To  remove  the  equipment,  the  procedure  was  reversed.  The 
cradle  was  lifted  from  the  water  by  securing  a  rope  between  it  and  a 
van,  then  the  van  was  driven  away  and  the  cradle  hoisted  up  along  the 


Diver  entry  and  exit  was  made  using  a  long  flat  aluminum 
ladder  tied  to  the  top  of  one  of  the  tank  walls. 

3.  Tests  and  Results 

The  facility  was,  at  the  time,  in  the  final  week  of  testing 
rubble  mound  breakwater  structures.  So  the  DPG  tests  conveniently 
"piggy-backed"  with  the  regularly  scheduled  tests.  Sets  of  eleven 
waves  of  1  meter  nominal  height  (3.28  feet)  and  5  second  period  were 
run  every  A  to  5  minutes.  Wave  height  and  period  were  monitored  by 
resistance  probes  and  recorded  on  a  strip  chart  located  inside  the 
trailer. 

Initially,  the  three  differential  channels  monitored  appeared 
ill-behaved  (Figure  V~7) •  The  fourth,  dP2,  had  already  been 
confirmed  inoperative  due  to  a  broken  arm  tubing  connector.  It  was 
then  discovered  that  the  end-cap  on  the  dP3  arm  had  accidentally  been 
left  on.  The  instrument  was  assembled  for  the  CERC  tests  with  the 
arm  isolation  diaphragms  flush  with  the  end  of  the  PVC  arms.  With 
the  arm  plugged,  the  diaphragm  could  communicate  with  the  ambient 
water  only  through  a  thin  crescent  part-way  around  the  diaphragm 
housing.  It  is  assumed  that  such  a  small  gap  dampened  the  dynamic 
pressure  field  of  the  passing  waves. 


FIGURE  V~7:  Sample  strip  chart  recording;  CERC  test.  Configi 
ara  end-cap  accidentally  left  on.  Chart  speed 
Calibration:  10  psia/V  and  0.2  psid/V. 
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Once  the  end-cap  was  removed,  the  third  channel  demonstrated 
sinusoidal  output  with  the  expected  TT/2  phase  difference  to  the 
absolute  signal  as  predicted  by  Equations  2.1  and  2.2.  The 
unexpected  characteristics  of  the  first  and  fourth  differential 
channels,  however,  remained  unexplained.  Both  channels  indicated  the 
passage  of  a  trough,  as  expected,  but  did  not  report  the  actual  crest 
passage  —  instead  suggesting  the  passage  of  two  crests  rr/4  before 
and  after  the  actual  passage.  This  phenomenon  was  noted  in  both 
orientations  of  the  DPG  tested.  It  is  thought  that  the  unexpected 
characteristics  of  these  differential  channels  could  have  been  due  to 
the  poor  positioning  of  the  sensors  at  the  unprotected  ends  of  the 
PVC  arms.  It  is  also  possible  that  the  arm  tubings  were  broken  at 
the  time  of  the  tests  or  that  the  sensitive  differential  gauges  were 
adversely  influenced  by  the  return  flow  pipes  to  which  the  cradle  was 
attached. 


The  instrument  was  originally  positioned  such  that  the  dPl 
and  dP3  arms  were  aligned  approximately  55  degrees  from  the 
centerline  of  the  tank,  or  wave  ray.  The  cradle  and  instrument  were 
next  rotated  about  20  degrees  such  that  arms  1  and  3  were 
approximately  35  degrees  from  the  wave  ray.  Figures  V~8  and  V-9  are 
strip  chart  recordings  from  each  of  the  configurations,  respectively. 
The  rotation  of  the  DPG  was  restricted  between  these  two  positions 
because  of  the  pipes  that  surrounded  either  side  of  the  cradle. 
These  pipes,  however,  offerred  the  only  structure  to  which  the  cradle 


FIGURE  V~8:  Sample  strip  chart  recording;  CERC  tes 


FIGURE  V~9:  Sample  strip  chart  recording;  CERC  test.  Configuration  II 
speed:  2  mm/sec.  Calibration:  10  psia/V  and  0.2  psid/V. 
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could  be  secured. 


Figure  V-10  compares  the  absolute  pressure  transducer  signal 
with  that  of  the  tank's  resistance  wave  probe  located  about  sixty 
feet  upstream  of  the  DPG  site.  The  resistance  probe  describes  the 
waves'  non-linear  profile  featuring  peaked  crests  and  long  shallow 
troughs.  Wave  period  is  5.0  seconds  and  the  height  is  between  3.0 
feet  (0.925  m)  and  3.36  feet  (1.025  m)  from  crest  to  trough.  The 
absolute  transducer  also  reports  a  period  of  5.0  seconds  and  a  wave 
height  of  about  3.2  feet  as  calculated  using  stream  function  theory 
(Dean, 1974).  Stream  function  theory  is  used  because  of  the  great 
non-linear  appearance  of  the  wave  form.  The  calculated  wave  height 
is  only  about  2.8  feet  if  linear  theory  is  assumed.  From  Eq.  2.1, 


amplitude  fr,st. 

trough 


pdyn  [crest 
[  trough 

1  Kq  [crest  ] 
[trough] 


(5.15) 


Here,  Y  “  62.23  lb/cu.ft.  (13  C) ,  and  Kp  is  approximately  0.874  and 
1.068  for  the  wave  crest  and  trough  respectively  from  stream  function 
theory  (Dean, 1974).  This  was  calculated  for  a  wave  of  five  second 
period  in  water  depth  10.2  feet  and  absolute  sensor  height  above  the 
bottom  of  4.25  feet.  The  intimation  of  static  water  level  is  also 


good.  The  no-wave  tare  value  was  steadily  reported  as  18.23  psia, 
which  corresponds  to  a  water  height  above  the  sensor,  from  Eq.  4.35 


of  5.91  feet. 
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z  =  P  -  y,  ,  -  P  -  P.  ,)/y 
s  mean  bf  bf  atm  bf 


(4.35) 


where  Y 


bf 

Zbf 

3 

’  atm 

Pbf 

Y 


55.6  lb/ cu. ft . 
18.5  inches 
14.74  psi 
0.34  psi 
62.23  lb/cu.ft. 


A  value  of  5.95  feet  is  expected  from  the  difference  of  the  tank 
water  depth,  10.2  feet,  and  the  sensor's  height  above  bottom,  4.25 
feet . 


An  estimate  of  the  wave  direction  was  attempted  from  the 
signal  of  dP3  using  wave  parameters  as  predicted  by  stream  function 
theory.  The  average  magnitude  of  the  channel's  signal  was  estimated 
and  the  maximum  differential  pressure  expected  along  the  wave  ray 
calculated  from  stream  function  theory.  The  angle©  between  the  wave 
ray  and  the  dP3  arm  is  then: 


0  =  cos'1 


C  d  P/d  x]  measured 

CdP/dx)max 


(5.16) 


The  wavenumber  k  was  calculated  using  stream  function  theory  for  a  5 
second  wave  in  the  water  depths  stipulated  earlier.  The  pressure 
response  function  Kp  was  similarly  estimated  for  a  depth  equal  to  the 
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mean  of  the  center  and  arm  sensor  vertical  positions  above  the 
bottom.  The  gage  length  dx  was  3.5  feet.  For  configuration  I,  the 
estimated  angle,  51.3°,  compares  well  with  the  55°  angle  expected. 
For  configuration  II,  the  estimated  angle  is  35.4°  compared  to  an 
expected  value  of  35°. 

The  experiment  at  CERC,  although  plagued  with  arm  sensor 
failures  and  highly  non-linear  waves,  demonstrated  the  capability  of 
the  instrument  to  measure  wave  direction.  The  tests  also  served  the 
purpose  of  defining  several  important  modifications  required  to 
improve  the  performance  of  the  DPG  system. 


E,  The  Second  DPG  Prototype 

The  DPG  was  returned  to  the  University  of  Delaware  laboratory 
and  partially  disassembled  to  install  the  new  hardware  modifications. 
The  fracture  of  the  arm  lines  of  dPl,  dP2,  and  dP4  were  confirmed. 
All  of  the  stainless  steel  tubing  exterior  to  the  water-tight 
instrument  cylinder  was  replaced  with  flexible  tubing  and 
back-filled,  and  the  system  was  calibrated  as  described  in  Section  C 
of  this  chapter.  New  PVC  arms  were  machined  to  house  the  flexible 
tubing  and  isolation  diaphragms  and  the  instrument  was  reassembled  as 
described  in  Appendix  D.  Care  was  taken  to  ensure  that  the  sensors' 
elastomer  diaphragms  were  well  back  from  the  ends  of  the  PVC  arms  and 
exposed  to  the  holes  in  the  sides  of  each  arm.  Some  dimensions  of 
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the  instrument,  reassembled  for  the  field  evaluation,  are  specified 
in  Appendix  C. 
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CHAPTER  VI 

FIELD  INSTALLATION  AND  EVALUATION 

A.  Installation  at  the  FRF 

The  newly  re-assembled  instrument  and  cradle  were  towed  on  a 
trailer  to  the  Coastal  Engineering  Research  Center  Field  Research 
Facility  at  Duck,  North  Carolina,  in  early  May  of  1982.  An 
installation  site  was  chosen  approximately  735  feet  directly  south  of 
the  seaward  end  of  the  1800  feet  long  research  pier.  This  site  fell 
safely  between  the  FRF 's  regular  bathymmetric  profile  lines,  allowed 
a  reasonalbe  safety  margin  of  length  for  cable  deployment  (18/0, 
exhibited  a  relatively  stable  20  foot  depth  contour  over  time,  and 
was  considered  to  be  sufficiently  far  from  the  large  scour  hole  at 
the  end  of  the  pier  to  avoid  any  related  refraction  effects.  The 
geographic  location  and  typical  nearshore  bathymetry  of  the  site  are 
illustrated  in  Figures  VI-1  and  VI-2. 

It  was  decided  to  install  the  DPG  using  the  facility's 
Coastal  Research  Amphibious  Buggy  (CRAB)  marine  vehicle.  The  CRAB  is 
a  three-wheeled  tripod  vehicle  capable  of  driving  across  the  beach 


FIGURE  VI-1:  Site  of  the  DPG  field  evaluation 
Birkemeier,  et.al.,  1981)- 


FIGURE  VI-2:  Nearshore  bathymetry  of  installation 


FIGURE  VI-3:  Preparing  for  field  inatallation:  the 

DPG  suspended  underneath  the  operator's  platform  of 
the  CRAB. 
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face  and  through  the  surf  zone.  The  instrument  was  secured  to  the 
cradle  and  hoisted  underneath  the  top  (operator's)  platform,  (Figure 
VI-3) .  The  cable  and  service  box  were  then  brought  to  the  end  of  the 
pier. 

On  May  14,  1982,  the  wave  activity  settled  sufficiently  to 
attempt  the  installation.  The  CRAB  was  driven  offshore  with  the 
instrument  and  positioned  about  200  feet  south  of  the  end  of  the 
pier.  Three  divers  and  a  boat  operator  entered  the  surf  in  a  Zodiac 
and  motored  underneath  the  pier.  The  service  box,  with  its  80  feet 
of  unarmored  cable  secured  inside,  was  lowered  into  the  Zodiac  and 
then  motored  towards  the  CRAB  —  pulling  the  cable  behind  it.  The 
service  box  was  carried  by  a  diver  below  the  rear  of  the  CRAB  and 
brought  underneath  the  center  of  the  vehicle.  In  what  became  the 
only  difficult  phase  of  the  operation,  the  service  box  was  pushed  up 
into  the  cradle  and  the  cable  mated  to  the  instrument.  Large,  long 
period  swell  made  this  maneuver  both  frustrating  and  often  dangerous 
—  as  the  sharp  corners  of  the  cradle  feet  posed  a  hazard  to  the 
divers.  The  cable  was  secured  to  the  bottom  of  the  CRAB  using  a 
prussick  knot  and  the  vehicle  was  driven  towards  the  installation 
site  —  pulling  out  the  cable.  Floats  were  placed  every  80  feet 
along  the  cable  from  the  end  of  the  pier. 
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Once  at  the  site,  originally  marked  by  a  buoy,  the  instrument 
was  lowered  into  the  water.  Divers  maneuvered  t"he  instrument  between 
the  two  rear  legs  of  the  CRAB  to  the  seafloor,  detached  the  lines 
from  the  CRAB  to  the  cradle  and  cable,  and  the  CRAB  drove  away  a 
short  distance.  A  four-foot  screw  anchor  was  turned  into  the  sandy 
bottom  approximately  three  and  a  half  feet  from  the  end  of  each  arm 
and  chain  was  attached  from  each  arm  to  its  respective  screw  anchor. 
The  chains  were  secured  taut  by  a  chain  binder  and  each  rested  at 
about  a  45  degree  angle  to  the  seafloor.  A  fifth  screw  anchor  and 
buoy  was  next  installed  twenty  feet  away  from  the  DPG  and  a  heavy 
line  secured  between  this  anchor  and  the  cradle.  Satisfied  with  the 
instrument  installation,  the  buoys  were  released  from  the  cable  and 
work  commenced  at  the  pier. 

The  cable  had  to  be  secured  under  the  pier  and  prevented  from 
chafing  on  the  cement  piling  jackets.  The  cable  descended  from  the 
pier  deck  inside  a  stainless  steel  case  along  the  northeastern-most 
piling.  A  steel  cable  was  looped  around  the  base  of  this  piling  and 
attached  to  a  wire  jacket  that  cinched  the  cable  under  tension.  To 
prevent  chafing  about  the  south  pilings  among  which  the  cable  passed, 
a  long  pipe  was  jetted  between  the  two  southeastern-most  pilings  and 
a  similar  steel  line  and  wire  jacket  attached  between  the  pipe  and 
cable.  The  cable  began  to  bury  itself  into  the  sand  within  hours  of 
the  deployment. 


The  cable  leads  were  interfaced  with  FRF  power  and 
data-logging  equipment  and  regular  data  recording  commenced  that 
evening . 

B.  Instrument  Orientation 

The  orientation  of  the  instrument  with  respect  to  magnetic 
north  was  determined  using  a  submersible  digital  compass  and 
conventional  diving  compasses.  The  digital  compass  was  mounted  at 
the  end  of  a  five-foot  length  of  aluminum  angle.  One  diver  recorded 
the  value  indicated  by  the  compass  while  another  diver  held  two  feet 
of  the  free  end  of  the  angle  onto  the  top  of  each  PVC  arm  of  the 
instrument.  In  this  way,  the  digital  compass  was  at  least  six  feet 
from  the  steel  cradle.  Compass  headings  were  recorded  for  each  arm 
twice  —  once  by  each  of  the  two  divers  —  and  then  averaged.  The 
orientation  of  each  arm  was  also  measured  by  securing  a  diver's 
compass  around  the  end  of  an  arm  using  the  compass  wrist  strap. 
Measurements  were  taken  three  tiroes  for  each  arm  using  two  different 
compasses  and  then  averaged.  All  of  the  values  recorded  for  each 
technique  along  each  arm  agreed  reasonably  well  within  a  technique. 
The  wrist  compass,  strapped  around  an  arm,  indicated  the  heading 
perpendicular  to  the  arm.  Accordingly,  ninety  degrees  was  subtracted 
from  each  wrist  compass  reading  to  give  the  heading  of  the  arm. 


As  an  experiment,  the  divers  compasses  were  placed  on  an  arm 


above  the  end  of  the  steel  channel  that  holds  the  arm.  The  compass 
was  then  slid  towards  the  end  of  the  arm.  The  heading  changed  as  the 
compass  was  slid  across  the  end  of  the  steel  channel  and  then 
remained  relatively  stable.  It  was  therefore  thought  at  the  time 
that  the  readings  taken  at  the  ends  of  the  PVC  arms  were  sufficiently 
far  from  the  steel  cradle  to  a/oid  any  bias  of  the  divers'  compasses. 
However,  as  can  be  seen  from  Table  6.1,  this  was  not  true.  Since 
each  arm  is  perpendicular  to  its  neighboring  arm,  each  of  the  compass 
readings  are  restrained  to  be  90  degrees  apart.  Although  the  digital 
compass  measurements  approximate  90  degree  separation,  those  of  the 
wrist  compasses  clearly  do  not. 


TABLE  6-1:  ORIENTATION  HEADINGS  OF  THE  DBG  ARMS 


ARM 

WRIST  COMPASS 
apparent  actual 

DIGITAL  COMPASS 
apparent  actual 

RESOLVED 

average 

■ 

228.33 

241.75 

238.80 

241.89 

241.82 

mm 

15503 

151-75 

153-90 

151-89 

151.82 

■ 

?8.00 

61.75 

61.20 

61.89 

61.82 

H 

32503 

331-75 

333-65 

331-89 

331.82 

*  Headings  are  with  respect  to  magnetic  north,  6/15/82. 


This  directional  anomaly  was  resolved  using  a  technique 
suggested  by  Dr.  Robert  Dean  (personal  conversation).  The  apparent 
(biased)  direction,/?,,  ,  indicated  by  the  compass  for  each  arm  n,  can 


Compass  needle 


FIGURE  VI-4:  Typical  influence  of  the  steel  cradle  upon  the 
magnetic  compasses  used  for  orientation  measurements. 


and 
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be  expressed  in  terms  of  the  actual  orientation  of  that  arm , /3n  » 
the  error  associated  with  that  arm,  en  ,  that  is  introduced  by  the 
metal  mass  of  the  cradle: 


8  '  =  8„  +  e 


(6.1) 


Since  the  compasses  were  placed  upon  each  arm  at  about  the  same 
distance  from  the  center  of  the  instrument,  and  if  one  assumes  that 
the  effect  of  the  metal  mass  upon  the  compass  readings  is  the  same 
for  each  arm  (i.e.,  radially  inward),  then  the  errors  associated  with 
each  of  two  col  linear  arms  should  be  equal  and  opposite  in  sense. 
(See  Figure  VI-4).  Accordingly,  the  sum  of  the  errors  associated 
with  the  reading  for  each  arm  should  equal  zero; 


4 

l  en  =  °  (6.2) 

n=l 

If  one  expresses  the  actual  orientation  of  each  arm  in  terms  of  the 
orientation  of  the  arm  with  the  smallest  value  of  apparent  direction. 
Equation  6.1  may  be  written  as: 


+  e. 


(6.3a) 


V  =  e3  +  90°  +  e4 


(6.3b) 


DPG  arm  was  taken  as  the  average  of  the  corrected  values  found  from 
the  digital  and  wrist  compass  measurements.  The  values  of  each 
correspond  to  the  heading  towards  which  the  end  of  each  arm  points 
with  respect  to  magnetic  north.  The  values  were  corrected  to  true 
north  using  the  1982  value  of  the  variation  for  Cape  Hatteras,  N.C. 

This  investigation  indicates  the  importance  of  careful 
redundant  checks  of  instrument  orientation  when  working  near  a  steel 
structure.  The  analysis  technique  outlined  here  requires  at  least 
one  pair  of  compass  measurements  made  on  opposite  sides  and 
equidistant  from  the  appoximate  center  of  the  metal  structure.  The 
investigation  suggests  that  simple  wrist  compasses  may  be  adequate 
for  determining  instrument  orientation  if  special  care  is  taken  during 
the  measurement  and  interpretation  of  the  results. 

C.  First  Month's  Inspection 

The  instrument  was  re-inspected  on  June  15,  1982  —  one  month 
after  installation.  The  cradle  and  instrument  had  settled  into  the 
seafloor  only  four  or  five  inches  with  minimal  disturbance  to  the  bed 
around  the  gauge.  The  cable  had  buried  itself  deeply  immediately 
outside  of  the  service  box.  The  system's  anti-fouling  agents  had 
thus  far  worked  excellently. 
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D.  Results  of  DPG  Data  Analysis 

1.  Introduction 

Multiple  sets  of  data  were  analyzed  to  determine  the 
integrity  of  the  gravity  wave  information  obtained  from  the  DPG. 
Although  output  was  available  from  the  facility's  HF  radar  system, 
there  were  no  other  functioning  in~s i tu  directional  wave  monitors  at 
the  FRF  during  the  weeks  after  the  instrument  installation.  Hence, 
DPG  directional  information  was  compared  with  visual  and  radar 
estimates.  Wave  height  and  period  could,  however,  be  compared  with 
CERC  Baylor  gauge  data  from  near  the  end  of  the  pier.  The  data  were 
also  inspected  for  irregularities  in  the  time  domain  that  could 
signal  system  malfunction. 

2.  Time  Domain  Signal  Analysis 

The  five  transducers'  records  corresponding  to  a  single  set 
of  wave  data  were  input  to  the  DPG  software  package  and  scanned  for 
signals  beyond  three  standard  deviations  from  the  mean  of  the  record. 
The  means  were  then  re-calculated  and  subtracted  from  the  records 
which  were  then  converted  to  pressure  signals.  (See  Sections 


IV. D. 1-3.) 


In  each  data  set  analyzed,  the  number  of  points  beyond  the 
three  standard  deviation  limit  was  satisfactorily  small  —  less  than 
one  percent  of  the  record,  as  expected.  The  mean  of  channels  dP2, 
dP3,  and  the  absolute  agreed  well  between  data  sets,  (possibly 
indicating  minimal  instrument  tilt  over  time  due  to  settling  —  see 
Section  IV. D. 4).  The  means  of  dPl  and  dP4  often  varied  between  sets, 
however.  It  was  expected  that  the  mean  of  dP4  might  vary  given  the 
signal  drift  observed  during  bench  calibration.  Similarly,  it  was 
not  surprising  to  encounter  difficulty  with  dPl  after  considering  its 
irregularities  on  the  bench. 

The  mean  values  of  the  absolute  records  predicted  tide  levels 
reasonably  well  as  is  shown  in  Table  6~2.  DPG  tidal  calculation  is 
compared  with  that  determined  by  a  stilling  well  located  near  the 
middle  of  the  pier.  Differences  might  be  attributed  to  temperature 
effects  upon  the  absolute  sensor,  (see  Section  V.C.6).  Atmospheric 
pressure  fluctuations  were  accounted  for. 

The  mean  values  of  dP2  and  dP3  differ  significantly  from 
these  two  channels'  tare  values  recorded  on  land.  In-situ  mean 
values  are  2.31  V  and  2.56  V  for  dP2  and  dP3  respectively,  and  0.60  V 
and  0.71  V  on  land.  Both  subaerial  and  subaqueous  checks  were  made 


with  the  same  1100 

foot  long  cable. 

The  subaerial  tare 

values 

were 

checked  on 

a  hot 

afternoon  and  the 

ocean  was  some  8°C 

cooler . 

But 

it  was  not 

1  ikely 

that  this  rise 

in  mean  voltage 

was  due 

to 

temperature  changes  since  one  expects  a  voltage  drop  with  a  decrease 


TABLE  6-2:  COMPABISON  OF  TIDE  LEVEL 


1.982 

m 

J. 

'IDE 

AMBIENT 

CONDITIONS 

DPG 

FRF 

T  (°C) 

P  (mt) 

DATE 

TIME 

O) 

(m) 

5/18 

0700 

-0.26 

-0.26 

16.8 

1018.0 

1300 

+0.12 

+0.18 

5/19 

0700 

-0.15 

-O.15 

I6.3 

1020.0 

1300 

-0.07 

-0.02 

5/20 

0700 

-0.05 

-0.05 

13.0 

1018.3 

1300 

-O.36 

-0.30 

5/21 

0700 

+0.13 

+0.12 

13.3 

1015.5 

1300 

-0.53 

-0.54 

5/22 

0700 

+0.46 

+0.43 

17.5 

1013.2 

1300 

-O.53 

-0.57 

5/23 

0700 

+O.58 

+0.57 

16.8 

1019.5 

1300 

-0.64 

-0.66 

5/24 

0700 

+O.56 

+0.56 

I6.5 

1018.0 

1300 

-O.56 

-0.62 

NOTES: 

FRF  data  is  from  stilling  well  at  pier. 
Tide  elevations  are  referenced  to  NGVD. 
DPG  Pbf(at>s)  =  1.32  psi. 

TIME  is  Eastern  Standard. 

Temperature  is  for  sea-water. 


in  temperature  in  the  present  DPG  configuration,  (Section  V.C.6). 
The  difference  could  indicate  slanted  installation  of  the  instrument. 
The  changes  in  mean  voltages  suggest  that  arms  2  and  3  could  be 
tilted  9.7°  and  11.8°  above  the  horizontal,  respectively;  (from 
Equation  A. 34).  This  could  be  confirmed  by  comparison  with  dPl  and 
dP4;  however,  neither  of  these  channels  show  a  clear  trend  in  the 
change  of  mean  value  between  land  and  sea. 

A  portion  of  the  time  record  of  one  data  set  is  shown  in 
Figure  VI-5.  The  mean  value  has  been  removed  from  each  record  and 
the  signs  of  dPl  and  dP2  have  been  reversed.  The  dPl  signals  have 
been  multiplied  by  the  channel  1  calibration  factor  of  1.65.  The 
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absolute  record  corresponds  to  pounds  per  square  inch  and  the  signals 
of  the  differential  records  to  pounds  per  square  inch  differential 
per  inch  of  gage  length.  It  is  immediately  obvious  that  the  signal 
of  dPl  is  extremely  erratic  with  high  frequency  energy  that  would 
alias  into  the  energy  spectrum. At  other  times,  it  agrees  in  form  with 
dP3  satisfactorily  although  it  is  shifted  downwards  because  of  the 
large  mean  value  subtracted  from  the  record.  This  large  mean  is 
apparently  the  result  of  the  large  positive-valued  high  frequency 
oscillations .  This  disturbance  is  intermittent  and  does  not  appear 
in  all  of  the  records  from  dPl.  These  high  frequency  oscillations 
have  never  been  observed  on  the  strip  chart  output  at  the  PRF,  and  so 
it  is  possible  that  the  problem  lies  within  the  digitization  of  the 
signal . 


The  records  of  dP2  and  dP4  agree  well  as  expected  for  two 
collinear  gauges.  There  are  slight  phase  and  amplitude  differences 
as  would  be  expected  since  the  two  gauges  are  separated  in  space  a 
finite  amount. 

There  is  reasonable  phase  agreement  between  the  dP2,  dP3,  and 
dPA  signals  with  those  of  the  absolute  gauge  —  as  best  as  can  be 
seen.  Again,  one  expects  a  ninety  degree  phase  difference  between 
the  absolute  and  differential  signals. 


The  "graininess"  of  the  absolute  gauge  signal  is  due  to 


resolution  limitations  of  the  transducer  and  digitizer. 

It  is  reasonable  to  conclude  from  inspection  of  the  raw  data, 
then,  that  the  dPl  records  should  not  be  included  in  the  spectral 
analyses,  and  that  less  faith  should  be  placed  in  dP4  because  of  its 
tendency  to  drift  (as  is  exhibited  by  its  changing  mean).  If  these 
two  channels  do  not  provide  valid  data,  it  is  not  possible  to 
generate  the  water  surface  curvature  terms  as  described  in  Chapter  IV 
with  any  confidence. 

3.  Record  Length 

In  beginning  the  conversion  of  the  data  to  the  frequency 
domain,  it  was  discovered  that  the  DEC-10  computer  system  could  not 
handle  the  17  minute  long  f our-samples~per-second  data  records,  (4096 
points) .  So  the  data  size  was  halved  by  averaging  two  adjacent 
points  in  the  time  series.  The  Fast  Fourier  Transform,  then, 
operated  upon  a  17  minute  record  with  an  effective  sampling  rate  of 
two  Hertz.  (The  corresponding  Nyquist  frequency  was  still  safely 
above  the  two  rad/sec  cut-off  frequency  used  in  the  analysis 
program.)  The  number  of  data  points  in  the  time  series  was  halved 
after  the  initial  scan  for  bad  points  and  calculation  of  the  record 


mean. 


It  was  discovered  that  this  modification  was  very  important. 
The  results  of  the  analysis  of  DPG  data  varied  considerably  between  a 
17  minute  long  record  with  an  effective  sampling  interval  of  0.5 
seconds  compared  to  the  analysis  of  an  8.5  minute  long  record  with 
0.25  second  sampling.  Further,  there  are  considerable  differences 
between  the  analysis  results  of  the  first  8.5  minutes  of  a  record 
compared  to  the  last  8.5  minutes.  This  variability  may  be  due  to 
sampling  and  the  confidence  limits  on  the  spectrum.  Figures  VI-6 
through  VI-8  are  copies  of  the  analysis  output  for  one  data  set 
analyzed  with  a  17  minute  record,  the  first  8.5  minutes  of  the 
record,  and  the  last  8.5  minutes,  respectively.  The  17  minute  record 
with  0.5  second  sampling  interval  block-averages  8  adjacent  frequency 
bands  compared  to  4  bands  for  the  0.25  second  sampled  records.  The 
spectra  developed  from  the  17  minute  long  record,  then,  is  more 
stable  than  for  the  8.5  minute  long  record,  yet  retains  the  same 
frequency  resolution.  Figure  VI~9  and  Figure  VI-10  illustrate  the 
energy  and  non-smoothed  peak-energy  directional  spectra  for  each 
case.  The  analysis  results  most  closely  match  other  observations 
made  at  the  FRF  when  the  full  17  minute  long  record  with  effective 
2-Hertz  sampling  is  utilized.  Table  6.3  lists  the  results  of  each 
type  of  analysis  and  compares  the  results  with  independent  FRF 
observations  over  a  period  of  four  days. 
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DEV 

YR 

MO 

DAY 

TIME 

1  1  1 

82 

6 

1  1 

7  0 

121 

82 

6 

1  1 

7  0 

141 

82 

6 

1 1 

7  0 

4096  POINTS  ANALYZED  AT  .50  SECOND  SAMPLING 

8  BANDS  BLOCK  AVERAGED 

GAGE  NUMBER  2 

8  LOW  truncations;  6  HIGH  truncations 

Mean*  2.299  S.Dev*  0.09652  Lower  Limit*  2.010  Upper  Limit*  2.589 

GAGE  NUMBER  3 

0  LOW  truncations;  8  HIGH  truncations 

Mean*  2.039  S.Dev*  0.19138  Lower  Limit*  1.465  Upper  Limit*  2.613 

GAGE  NUMBER  5 

O  LOW  truncations;  1  HIGH  truncations 

Mean*  2.364  S.Dev-  0.03900  Lower  Limit*  2.247  Upper  Limit-  2.481 

TIDE  LEVEL  IS  0.824  FEET. 

POSSIBLE  INSTRUMENT  TILT  (deg) 

DP2  DP3 

-0.061  -3.365 


Significant  Wave  Height  (ABS  gage)  *  4.32  feet. 

Effective  period  *  7.817  seconds. 


Maximum  Energy  located  at  bands; 

11.070  exceeds  average  energy  by 
10. 189  exceeds  average  energy  by 
12.118  exceeds  average  energy  by 
13.386  exceeds  average  energy  by 
8 . 790  exceeds  average  energy  by 
8 . 225  exceeds  average  energy  by 


715.00  X 
266. 13  X 
194 . 16  X 
159.03  X 
115.97  X 
108.46  X 


Compass 

heading  from 

which  .aves  propogate 

Band 

1 1 . 070 

secs : 

Greatest 

Energy 

at 

73 

Band 

lO. 189 

secs : 

Greatest 

Energy 

at 

71 

Band 

12.118 

secs : 

Greatest 

Energy 

at 

67 

Band 

13.385 

secs : 

Greatest 

Energy 

at 

69 

Band 

8.790 

secs : 

Greatest 

Energy 

at 

69 

Band 

8.225 

secs : 

Greatest 

Energy 

at 

75 

deg 

deg 

deg 

deg 

deg 

deg 


FIGURE  VI-6:  DPG  analysis  result*  using  the  full  record 
length  of  17  ninutes  with  effective  2  hertz  sampling. 
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DEV 

YR 

MO 

DAY 

TIME 

1  1  1 

82 

6 

1  1 

7  0 

12  1 

82 

6 

1  1 

7  0 

141 

82 

6 

1  1 

7  0 

2048  POINTS  ANALYZED  AT  .25  SECOND  SAMPLING 

4  BANDS  BLOCK  AVERAGED 


GAGE  NUMBER  2 

7  LOW  truncations;  2  HIGH  truncations 

Mean-  2.299  S.Dev*  0.10032  Lower  Limit*  1.998  Upper  Limit*  2.600 


GAGE  NUMBER  3 

O  LOW  truncations:  4  HIGH  truncations 

Mean*  2.037  S.Oev*  0.18856  Lower  Limit*  1.471  Upper  Limit*  2.603 


GAGE  NUMBER  5 

0  LOW  truncations;  4  HIGH  truncations 

Mean*  2.360  S.Dev*  0.03693  Lower  Limit*  2.250  Upper  Limit*  2.471 


TIOE  LEVEL  IS  0.743  FEET. 


POSSIBLE  INSTRUMENT  TILT  (deg) 
DP2  DP3 

-O  065  -3.377 


Significant  Wave  Height  (ABS  gage)  *  4.24  feet. 


Effective  period  =  7.206  seconds. 


Maximum  Energy  located  at  bands: 

11.011  exceeds  average  energy  by 
12.047  exceeds  average  energy  by 
8.192  exceeds  average  energy  by 
10. 139  exceeds  average  energy  by 
9.394  exceeds  average  energy  by 
6.872  exceeds  average  energy  by 


435. 16  % 
343.94  ’/. 
233.40  % 
138.87  % 
136.97  % 
127.69  % 


Compass  heading  from  which  waves  propogate 


Band  11.011 
Band  12.047 
Band  8.192 
Band  10. 138 
Band  9 . 394 
Band  6.872 


secs;  Greatest 
secs;  Greatest 
secs :  Greatest 
secs:  Greatest 
secs:  Greatest 
secs :  Greatest 


Energy 

at 

77 

Energy 

at 

77 

Energy 

at 

71 

Energy 

at 

81 

Energy 

at 

61 

Energy 

at 

61 

deg 

deg 

deg 

deg 

deg 

deg 


FIGURE  VI-7:  DPG  analysis  results  using  the  first  8.5  min¬ 
utes  of  the  record,  (4  hertz  sampling). 
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DEV 

YR 

MO 

DAY 

TIME 

1  1  1 

82 

6 

1  1 

7  0 

121 

82 

6 

1  1 

7  0 

14  1 

82 

6 

1  1 

7  0 

2048  POINTS  ANALYZED  AT  .25  SECOND  SAMPLING 

4  BANDS  BLOCK  AVERAGED 


GAGE  NUMBER  2 

O  LOW  truncations;  4  HIGH  truncations 

Mean=  2.300  S.Dev*  0.09259  Lower  Limit*  2.022  Upper  Limit*  2.578 


GAGE  NUMBER  3 

O  LOW  truncations;  4  HIGH  truncations 

Mean*  2.041  S.Dev*  0.19418  Lower  Limit*  1.458  Upper  Lfmit*  2.623 


GAGE  NUMBER  5 

0  LOW  truncations;  0  HIGH  truncations 

Mean*  2.367  S.Dev*  0.04066  Lower  Limit*  2.246  Upper  Limit*  2.489 


TIDE  LEVEL  IS  O  904  FEET. 


POSSIBLE  INSTRUMENT  TILT  (deg) 
DP2  DP3 

-0.058  -3.351 


Significant  Wave  Height  (ABS  gage)  =  4.43  feet. 

Effective  period  =  8.245  seconds. 


Maximum  Energy  located  at  bands: 

11.011  exceeds  average  energy  by 
10.139  exceeds  average  energy  by 
13.299  exceeds  average  energy  by 
8.752  exceeds  average  energy  by 
7.699  exceeds  average  energy  by 
12.047  exceeds  average  energy  by 


931.01  % 
301.30  % 
291.62  % 
199.04  % 
185.02  % 
118.44  % 
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Energy 
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at 
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FIGURE  VI-8:  DPG  analysis  results  using  the  last  8.5  min¬ 
utes  of  the  record,  (4  hertz  sampling). 


FIGURE  VI~9:  Block-averaged  energy  spectra  of  water  surface 
displacement  for  a  17  minute  long  record  and  8-5  minute 
long  records. 


ENERGY -SPECTRA  of  DISPLACEMENT  TERM 


ENERGY  (In2  sec  ) 


FIGURE  VI-10:  Non-weighted  directional  spectra  of  peak 
energy  band  (11.0  seconds)  for  a  17  minute  long  record 
and  8.5  minute  long  records. 
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TABLE  6- 3s  COMPARISON  OF  DPG  RESULTS  WITH  FRF  OBSERVATIONS 
8  June  -  11  June,  1982 


DATE 

TIME 

1 

first 
half  of 
record 

-  DPG  — 

second 
half  of 
record 

“1 

FULL 

RECORD 

FRF 

OBSERVATIONS 

6/8 

0700 

dirxn 

59° 

59° 

59° 

58°  (50°)b 

Hsig 

1.48 

1.30 

1.39 

1 .41  c 

period 

8.75 

8.19 

8.23 

8. 00c 

1300 

dirxn 

59° 

59° 

59° 

- 

Hsig 

1 .31 

1.26 

1.28 

l.33c 

period 

8.75 

6.39 

8.79 

9 -57c 

6/9 

0700 

dirxn 

69° 

65° 

67° 

68°  ?  (50°)b 

Hsig 

1-37 

I.32 

I.34 

1.33c 

period 

9-39 

10.14 

9.44 

9.66  c 

6/10 

0700 

dirxn 

630 

69° 

65° 

66°  (55°)b 

Hsig 

1.31 

1.16 

1.23 

1.56c 

period 

10.14 

11.01 

10.19 

10.56c 

1300 

dirxn 

69° 

91° 

85° 

- 

Hsig 

1.69 

1.74 

1.71 

1.82  c 

period 

11.01 

12.05 

11.07 

10.34  c 

6/11 

0700 

dirxn 

77° 

73  ° 

73° 

69°  (60°)„ 

Hsig 

1.29 

1.35 

1.32 

1.53c 

period 

11.01 

11.01 

11.07 

10.89c 

NOTES: 

"dirxn"  listed  is  the  principal  direction  (peak  energy),  TN. 
"Hsig"  is  the  significant  wave  height  in  meters. 

"period"  corresponds  to  the  frequency  band  of  greatest  energy. 
a  CERC  Radar  (+  2°  ) 

b  Visual  estimate  from  the  end  of  the  pier 
cCERC  Baylor  Gauge  near  the  end  of  the  pier 
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5,  Directional  Spectra  Smoothing 

The  analysis  results  of  six  data  sets  from  the  DPG  over  a 
four  day  period  are  presented  in  Appendix  A.  The  analysis  was 
carried  out  using  the  absolute  pressure  gauge  and  differential  gauges 
dP2  and  dP3. 

The  smoothed  directional  spectra  of  the  six  peak  energy  bands 
are  illustrated  for  each  data  set.  In  the  development  of  these 
results,  the  directional  spectra  is  expressed  as  the  partial  Fourier 

series  in  terms  of  the  wave  directions  as  in  Equation  4.  14: 

N  N 

s(a,0)  =  w  a  (o)  +  £  w  a  cos  n6  +  I  Wnbn  sin  n6 

°  n=l  n  n  n=l 

where  N-2.  This  is  smoothed  by  the  weighting  function 


-J  COS4  ~(0'  -  0) 


(6.6) 


such  that  w0  -  1 ,  w^  ”  2/3,  and  w2  -  1/6.  Such  smoothing  eliminates 
negative  energy  side-lobes  in  the  directional  spectra  by  imposing  an 
always-positive  function  of  the  same  angular  frequency  upon  the 
directional  distribution  (Longuet-Higgins ,  Cartwright  and  Smith, 
1963).  The  function 


4 

B  cos  b0 


(6.7) 
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is  first  considered,  and  then  expressed  in  complex  form: 


=  Re 


ibfl  -ib0 
e  +  e 


(6.8) 


Expanding, 


=  R 


B  4ib9  .  2ib0  ^  .  -2ib0  4ib6 
—(e  +  4e  +  4e  +  e  +6) 

lb 


*] 


Equation  6.7  can  be  written  as: 

g 

=  —(cos  4b0  +  4cos  2b0  +  3) 

O 

Letting  b”l/2,  and  factoring  3,  the  function  becomes: 


=  B (1  +  cos  0  +  cos  20)  (6.9) 

o  3  3 


The  weighted  angular  distribution  must  satisfy  the  requirement 


S(o) 


•2n 

s(o,0)  d0 

'0. 


(6.10) 


to  retain  the  original  energy  calculated  in  the  spectra.  Hence,  B 
must  equal  8/3.  The  coefficients  w0,  w, ,  and  w2  are  taken  as  1,  2/3, 
and  1/6  in  accordance  with  the  definition  of  the  Fourier  coefficients 


described  in  Equations  4-18  through  4-24.  If  the  first  seven 
directional  Fourier  coeff icients ,  (i.e.,  N=3) ,  were  developed,  then 
the  directional  spectra  might  be  smoothed  by  a  weighting  function 

fcos6i(9'-91  (6.11) 

such  that  w0  =  1,  w1  *  66/77,  w2  *  15/28,  and  w3  *  5/21. 

6.  Redundancy  of  dP2  and  dP4  Measurements 

If  differential  pressure  gauges  dP2  and  dP4  work  correctly, 
then  the  results  of  directional  analysis  of  wave  data  using  either  of 
these  two  collinear  gauges  should  be  nearly  identical.  To  check 
this,  a  number  of  wave  records  were  analyzed  using  the  signals  from 
the  absolute  pressure  channel,  dP3,  and  dP2,  and  then  the  signals 
from  the  absolute,  dP3,  and  dP4  channels.  Two  of  the  results  are 
shown  in  Figures  VI-12  and  VI-13.  Figure  VI-12  illustrates  the 
energy  spectrum  and  the  non-smoothed  directional  spectra  of  the  first 
two  peak  energy  bands  for  one  data  set  using  both  combinations  of 
gauges.  The  energy  spectrum,  taken  from  the  signals  of  the  absolute 
gauge,  is  the  same  for  either  gauge  combination.  There  is  good 
directional  agreement  in  this  case  of  well-defined  swell.  Figure 
VI-13  illustrates  the  energy  spectrum  and  non-smoothed  directional 
spectra  of  the  first  two  peak  energy  bands  for  another  data  set.  The 


FIGURE  VI-12:  Non-weighted  directional  apectra  of  the  firat  two 
peak  energy  banda  uaing  ABS,  dP2,  dP3  gauge  combination  and 
ABS,  dPA ,  dP3  gauge  combination.  Caae  1:  Well-defined  awell. 
Block-averaged  water  aurface  diaplacement  energy  apectrum  ahown 
in  inaet. 
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FIGURE  VI-13:  Non-weighted  directional  spectra  for  the  first  two 
peak  energy  bands  using  ABS,  dP2,  dP3  gauge  combination  and 
ABS,  dPA ,  dP3  gauge  combination.  Case  2:  Omnidirectional  sea. 
Block-averaged  water  surface  displacement  spectrum  shown  in 
inset. 
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agreement  between  different  combinations  of  gauges  is  not  as  good  for 
this  set  that  represents  a  "ragged"  sea  state.  This  is  most  probably 
because  each  of  the  two  gauges,  separated  a  finite  distance, 
independently  measure  more  slightly  different  wave  activity  than  in 
the  first  data  set  shown.  For  such  a  confused  sea  state  of 
relatively  short  period  waves,  it  was  not  possible  to  verify  which 
channel,  dP2  or  dP4,  gave  the  more  accurate  estimate  of  wave 
direction  (if  indeed  there  was  any  one  principal  direction). 
Accurate  radar  and  visual  observations  are  practically  impossible 
during  such  conditions.  As  expected,  the  omnidirectional  energy 
spectra  generated  from  the  DPG  data  for  this  confused  sea  state 
agrees  in  form  with  that  of  the  computer-simulated  case  of  one  wave 
frequency  with  two  directions  as  shown  in  Figures  IV-10  and  IV-11. 

7.  Energy  Spectra  from  Differential  Pressure  Gauges 

Figure  VI-14  illustrates  one  record's  wave  energy  spectra 
generated  from  the  absolute  pressure  gauge  and  the  water  surface 
slope  energy  spectrum  from  each  of  the  two  perpendicular  differential 
pressure  gauges  dP2  and  dP3.  The  record  was  taken  during  wave 
activity  from  the  northwest  —  the  direction  in  which  DPG  arm  3  is 
oriented  towards.  Accordingly,  one  observes  that  there  is 
considerably  more  energy  contained  in  the  slope  spectrum  from  dP3 
than  from  dP2.  This  is  because  there  are  relatively  small  slope  (and 
pressure)  differences  along  the  wave  crests  and  troughs  which  are 


FIGURE  VI-14:  Block-averaged  power  apectra  of  water  surface  disp- 


placenent  (developed  from  the  absolute  pressure  gauges)  and 
6  water  surface  slope  (developed  from  differential  pressure 

gauges  dP2  and  dP3) . 
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aligned  with  arm  2  whereas  there  are  relatively  large  slope  (and 
pressure)  differences  along  the  wave  ray  —  aligned  with  arm  3.  (See 
Figure  VI-15.) 


As  is  illustrated  in  Figure  — . ---  —  11  - 

VI-14,  one  expects  that  the  wave 

frequencies  of  greatest  energy  should  s'ssSsSy'^^  s'  s 

be  the  same  for  both  the  water  sssssS^\.S//^s'  s's^ 

surface  displacement  and  slope  terms.  s'  s'  s'  . Jsy^sssSssSs' 

It  is  instructive  to  note,  however,  s'  s's^SSssss'  s' 

that  greater  energy  appears  in  the  s' /sss/%  s'  s 

higher  frequency  bands  of  the  slope  SsSsSSSS s'  s'  S^ s, 

spectra  than  in  those  of  the  wave  FIGURE  VI-15;  Wave  crest  ap¬ 

proximately  aligned  with 
spectrum.  This  too  is  an  expected  jjpj  arB 

result  since  the  energy  spectrum  of 

of  slope  (developed  from  a  differential  pressure  signal)  differs  from 
that  of  water  surface  displacement  (generated  from  the  absolute 
pressure  signal)  by  a  factor  of  wavenumber-squared.  The  energy 
spectrum  of  the  water  surface  displacement  is  developed  through  the 
frequency-by-frequency  division  of  the  dynamic  pressure  energy 
spectrum  by  the  specific  gravity  of  seawater  (squared)  and  the  square 
of  the  pressure  response  function,  Kp.  The  water  surface  slope 
energy  spectrum  is  similarly  developed  through  the  division  of  the 
differential  pressure  energy  spectrum  by  the  squares  of  the  specific 


gravity  of  seawater  and  the  pressure  response  function.  The  energy 


spectrum  of  the  water  surface  displacement  is  then  as  Equation  4.46: 


S 

nn 


(a) 


(4.46) 


and  the  energy  spectra  of  the  water  surface  slopes  are  as  Equations 
4.47  and  4.48: 
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Since  wavenumber  increases  with  wave  frequency,  energy  leaked  into 
the  higher  frequency  differential  pressure  bands  is  enhanced  in  the 
generation  of  the  water  surface  slope  spectrum  relative  to  the 
displacement  spectrum.  It  is  for  this  reason  that  the  frequency 
bands  of  greatest  energy  are  selected  for  directional  analysis  from 
the  water  surface  displacement  spectrum.  (Spectral  leakage  is  the 
"spilling"  of  energy  from  the  frequencies  that  are  actually  present 
in  the  record  to  other  frequencies.  It  occurs  naturally  in  the 
analysis  of  a  finite  length  of  record  due  to  the  inability  of  the  FFT 
to  exactly  resolve  the  spectral  content  of  a  time  record  on  the 
frequency  axis.) 


8-  Directional  Estimates  without  the  Absolute  Pressure  Gauge 


Estimates  made  of  the  wave  direction  without  the  absolute 
pressure  gauge  are  generally  unreliable.  Typical  results  calculated 
using  Equation  4.52  and  the  dP2  and  dP3  signals  for  the  data  sets 
illustrated  in  Appendix  A  are  listed  at  the  end  of  the  appendix. 
Equations  4.52  and  4.54  are  valid  if  there  exists  only  one  wave 
direction  per  frequency.  When  waves  of  one  frequency  arrive  at 
slightly  different  directions  and  spectral  frequency  bands  are 
contaminated  by  leakage,  the  validity  of  Equations  4.52  and  4.54 
fails.  As  a  general  statement,  it  can  be  said  that  these  equations 
do  not  usually  provide  reliable  estimates  of  the  principal  wave 


direction  in  typical  ocean  wave  analysis.  The  estimates  might  be 


9.  Arithmetic  Development  of  the  Curvature  Terms 


For  investigative  purposes,  it  was  attempted  to  create  the 
water  surface  curvature  terms  as  described  earlier.  Typically,  this 
was  only  done  along  the  y-axis  using  the  dP2  and  dP4  signals  since 


FIGURE  VI-16:  Typical  portion  of  time  series  of  dP2  and  dP4 
"slope"  signals  with  the  associated  dP4-dP2  "curvature" 
signals.  Signal  of  dP4  is  dashed  line. 


the  digitized  x-axis  signal  of  dPl  was  erratic.  The  terms  were 
arithmetically  created  in  the  time  domain  through  the  difference  of 
the  dP2  and  dP4  signals  after  adjacent  values  in  each  channel's 


record  were  averaged  and  the  number  of  points  in  the  17  minute  long 


record  halved.  A  portion  of  the  resulting  y-axis  curvature  term  in 
the  tiae  doaain  is  shown  in  Figure  VI-16.  The  tiae  series  was 
processed  by  the  FFT  and  the  energy  spectra  developed.  Figure  VI-17 
illustrates  the  energy  spectrum  of  the  curvature  term  compared  to 
that  of  the  displacement  term.  The  wave  frequencies  corresponding  to 
peak  energy  levels  of  the  data  set  correspond  well  between  the 
spectra.  The  curvature  spectrum  contains  greater  energy  in  the  high 
frequency  bands  than  the  displacement  spectrum.  This  is  expected 
since  the  curvature  spectrum  differs  from  the  displacement  spectrum 
by  a  factor  of  the  fourth  power  of  the  wavenumber. 

Despite  the  reasonable  appearance  of  the  curvature  spectra, 
the  directional  Fourier  coefficient  b3 ,  developed  using  the  water 
surface  curvature  along  the  y-axis,  was  at  least  an  order  of 
magnitude  larger  than  the  other  directional  coefficients  for  each 
case  tested.  The  dominance  of  this  term  would  create  three  high 
energy  peaks  in  the  dirctional  spectra  and  invalidate  directional 
estimates . 

The  curvature  terms  were  also  developed  in  the  time  domain 
before  adjacent  values  in  the  record  were  averaged.  The  results  were 
similar. 

The  inability  to  generate  higher  order  directional  Fourier 
coefficients  through  the  arithmetic  creation  of  the  water  surface 
curvature  terms  is  not  surprising.  The  DPG  was  developed  upon  the 
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FIGURE  VI-17:  Block-averaged  power  spectra  of  water  sur¬ 
face  displacement  and  water  surface  curvature  terms. 
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premise  that  there  is  inherent  error  in  the  creation  of  small-valued 
water  surface  slope  terms  through  the  subtraction  of  two  large 
pressure  values.  Similarly,  the  attempt  was  made  in  this  case  to 
create  the  very  small-valued  water  surface  curvature  terms  through 
the  subtraction  of  two  small  differential  pressure  values. 


CHAPTER  VII 


RECOMMENDATIONS 


During  the  design,  development,  testing,  installation,  and 
evaluation  of  the  DPG  directional  wave  monitor  system,  some  problem 
areas  and  phenomena  of  interest  were  cited  for  further  investigation 
and/or  development.  Foremost,  thorough  comparisons  should  be  made  of 
the  DPG  analysis  results  with  those  of  other  in~s i tu  directional  wave 
monitors  at  the  Field  Research  Facility  —  particularly  those 
instruments  that  typically  generate  water  surface  slope  terms  through 
the  subtraction  of  signals  from  adjacent  pressure  sensors.  Also,  the 
digitization  of  the  dPl  signals  from  the  DPG  should  be  checked  and/or 
differential  transducers  1  and  4  should  be  repaired  or  replaced. 

There  are  certain  design  characteristics  of  the  first  DPG 
prototype  that  might  be  investigated  and  improved.  One  might 
re-design  the  fuselage  to  align  the  differential  pressure  sensors 
exactly  upon  the  x~  and  y-  axes.  In  the  first  prototype,  as 


installed  at  the  FRF,  the  dP2,  dP3,  and  dP4  sensors  are  misaligned  by 
one  to  two  degrees  as  discussed  in  Section  III.C.l.  The  instrument 
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might  be  re-designed  such  that  the  arm  and  center  sensors  are  also  at 
the  same  elevation  above  the  seafloor.  However,  for  the  first  DPG 
prototype,  the  nine  inch  difference  in  height  between  the  sensors 
leads  to  very  small  error  in  the  calculation  of  an  average  pressure 
response  function  used  in  the  development  of  the  water  surface  slope 
terms.  (See  Section  III.C.2.)  In  water  depths  as  shallow  as  10  feet 
(3.1  m) ,  the  error  is  still  small  —  less  than  three  percent  for 
waves  of  three  second  period. 

The  instrument  might  be  re~designed  such  that  the  lengths  of 
each  tubing  used  in  the  sensor  systems  are  identical.  This  would 
help  ensure  that  the  response  time  of  each  channel  is  identical  and 
that  there  is  no  difference  in  response  time  for  the  loading  of  an 
arm  sensor  compared  to  a  center  sensor.  This  would  eliminate 
differential  pressure  drift  induced  by  changes  in  ambient 
temperature,  thus  disabling  the  instrument  to  measure  ambient  water 
temperature. 

The  arrangement  of  the  transducers  inside  the  water-tight 
cylinder  might  be  re-designed  to  allow  servicing  of  individual 
pressure  transducers  without  the  need  to  bring  the  entire 
instrument  to  the  surface.  Such  a  design  change  could,  however,  add 
additional  problems  to  the  dependability  of  the  transducers  since 
this  would  most  probably  require  the  addition  of  underwater-pluggable 
connectors  for  each  transducer.  The  present  instrument  utilizes  one 
such  connector  that  is  mated  only  above  water  in  routine  maintenance. 


For  DPG  operation  in  water  depths  of  20  feet,  the  absolute 
pressure  transducer  should  be  changed  from  50  psia  rated  capacity  to 
30  or  35  psia  to  improve  resolution. 

Different  designs  of  the  isolation  sensor  diaphragms  could  be 
considered  in  order  to  create  a  sensor  that  is  easier  to  back-fill 
and  less  difficult  to  machine.  When  calibrating  the  sensors,  an 
oscillating  pressure  pump  used  in  conjunction  with  the  calibration 
system  described  in  Section  V.C.l  might  better  estimate  the  system's 
dynamic  loading  characteristics.  Temperature  effects  upon  the 
transducers  and  their  respective  sensors  should  be  investigated  more 
thoroughly  during  instrument  calibration. 

The  effects  of  vortex  shedding  about  the  sensor  diaphragms 
might  be  investigated.  The  sensors  could  be  re-configured  within  the 
PVC  arms  and  fuselage  to  minimize  vortex  shedding  effects  if  they  are 
found  to  be  of  sufficient  magnitude  for  some  application  of  the  DPG. 
Tests  in  the  large  wave  tank  at  the  Coastal  Engineering  Research 
Center  suggest  that  the  configuration  of  the  sensors  within  the  arms 
is  important  for  accurate  directional  estimates,  (Section  V.D.3). 

Although  a  tripod  was  at  one  time  considered  for  the 
instrument  cradle,  the  four-legged  cradle  used  for  the  DPG  at  the 
Field  Research  Facility  proved  stable  and  easy  to  maneuver  in  the 
installation  and  first  weeks  of  operation  in  the  field.  Future  DPG 
systems  may  retain  this  four-legged  cradle  or  experiment  with  a 
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tripod.  Any  sharp  corners  on  the  cradle  should  be  rounded  and 
smoothed  for  safety. 

The  DPG  might  be  better  protected  against  high  voltage 
transients,  (i.e.,  lightning),  by  bonding  the  steel  armor  of  the 
cable  to  the  water-tight  instrumentation  cylinder  at  the  seaward  end 
and  to  the  ground  system  of  the  FRF  at  the  landward  end  of  the  cable, 
(J.  Anderson  Plumer,  Lightning  Technologies,  Inc.,  personal 
conversation) . 

Various  configurations  of  the  DPG  arms  used  for  differential 
pressure  measurements  might  be  considered  for  investigative  purposes. 
Intuitively,  one  might  expect  that  the  most  effective  directional 
measurements  made  by  the  DPG  would  be  for  the  case  of  a  maximum 
signal  reported  by  each  differential  pressure  transducer.  This 
requires  that  the  differential  transducers  be  oriented  45  degrees  to 
the  wave  ray  or  crest.  In  this  way,  the  peak  signals  from  each 
transducer  would  be  of  the  same  amplitude,  and  for  sufficiently  large 
wave  activity,  above  the  noise  level  of  each  transducer  channel.  If 
the  incident  wave  direction  at  a  site  was  normally  distributed  about 
the  shore  normal,  then  one  might  expect  the  most  effective 
directional  measurements  to  be  made  by  a  DPG  with  its  arms  oriented 
45  degrees  to  the  shore  normal.  The  DPG  prototype  was  placed  at  the 
FRF  (with  little  consideration  for  orientation)  such  that  the  x~axis 
(dPl  and  dP3)  arms  and  the  y~axis  (dP2  and  dP4)  arms  are  17  and  73 
degrees  from  the  shore  normal,  respectively.  Various  configurations 
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of  the  arms  with  respect  to  one  another  might  be  considered.  The 
angle  between  arms  might  be  varied  to  "tune"  the  DPG  to  the  incident 
wave  field  —  much  the  same  as  a  directional  antenna.  Changes  in  arm 
orientation  could  also  permit  the  estimation  of  higher  order 
directional  Fourier  coefficients,  (see  Appendix  E) . 

In  future  versions  of  the  DPG  system,  wave  data  could  be 
transmitted  to  ship  or  shore  by  telemetry  —  eliminating  difficult 
cable  installations.  Data  might  also  be  integrated  in~situ  and 
reported  in  concise  segments  by  telemetry  or  on  tape.  Present 
coastal  field  measurements  are  often  made  over  long  periods  of  time, 
where  the  results  of  such  measurements  represent  coastal  processes 
that  are  essentially  integrated  over  time.  In  such  cases,  specific 
hour-by-hour  wave  data  may  not  be  necessary  or  even  desirable. 
Concise,  integrated  wave  information  is  just  as  useful.  Software 
capability  on  site  with  the  DPG  could  be  developed  to  integrate  large 
amounts  of  wave  data  and  report  one  or  two  numbers  (say,  measures  of 
P,s  or  Sxy)  to  characterize  the  wave  climate  over  time,  (Robert  Dean, 
personal  conversation). 

A  technique  has  also  been  considered  to  directly  measure 
water  surface  curvature  using  differential  pressure  gauges.  This 
measurement  technique  could  be  developed  and  evaluated  in  the  hope  of 
accurately  developing  higher  orders  of  directional  Fourier 
coefficients  from  directly  measured  data. 
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Various  types  of  directional  data  analysis  schemes  should  be 
investigated.  Among  them  are  the  Maximum  Likelihood  Method,  (Clarke 
and  Gedling, 1981 ;  Jefferys,  et.al.,1981;  Goda,1981;  LaCoss , 1971) . 
and  the  estimate  of  higher  order  directional  Fourier  coefficients  as 
alluded  to  in  Section  II.B.l  and  discussed  in  Appendix  E.  Proper 
conditioning  of  DPG  data  might  also  allow  directional  estimates  made 
without  the  signal  of  the  absolute  pressure  gauge. 

In  all  future  developments  of  the  DPG,  as  with  any 
well-designed  ocean  instrument,  careful  consideration  must  be  given 
to  the  problems  of  corrosion,  biological  fouling,  scour  and 
instrument  settlement  into  the  seabed,  installation  and  maintenance, 
electronics  failure  and  noise,  and  interference  from  local  marine 
activity.  The  limitations  of  directional  wave  spectra  analyses 
should  always  be  recognized  when  dealing  with  the  data  obtained  from 


the  instrument. 
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CHAPTER  VIII 
CONCLUSIONS 


This  thesis  has  presented  a  description  of  the  design, 
development,  and  evaluation  of  a  new  type  of  directional  wave 
monitor.  The  work  undertaken  in  the  preparation  of  this  thesis 
indicates  that  differential  pressure  transducers  used  with  an 
absolute  pressure  transducer,  as  employed  in  the  DPG  directional  wave 
monitor,  appear  capable  of  generating  accurate  estimates  of  wave 
direction  using  the  first  five  directional  Fourier  coefficients  as 
determined  from  direct  measurements  of  water  surface  displacement  and 
slope.  This  conclusion  is  based  upon  the  comparison  of  analyzed  data 
from  the  DPG  with  estimates  from  the  CERC  radar  located  at  the  CERC 
Field  Research  Facility,  Duck,  N.C.  The  relatively  small  dimensions 
and  light  weight  of  the  DPG  system,  (two-fifths  the  size  and  less 
than  one-seventh  the  weight  of  conventional  pressure  sensor  arrays), 
allowed  a  smooth,  relatively  simple  field  installation.  The  DPG,  as 
designed,  should  similarly  allow  for  easy  retrieval  and  re-deployment 
of  the  instrument  for  maintenance  purposes. 
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It  has  been  demonstrated  that  the  differential  pressure 
transducer  can  be  considered  to  be  an  inherently  more  effective 
instrument  for  determining  wave  direction.  The  design  is  well  suited 
to  the  measurement  of  water  surface  slope  and  the  response  of  surface 
slope  with  depth  is  a  maximum  over  frequencies  that  are  more 
representative  of  typical  ocean  gravity  waves  than  the  response  of 
conventional  absolute  pressure  transducers. 

In  the  development  of  the  DPG  system,  it  was  discovered  that 
careful  attention  must  be  given  to  the  diameter  of  the  tubing  and  to 
the  nature  of  the  back-filling  fluid  used  in  the  transducers' 
sensors.  This  is  to  ensure  a  small  response  time  for  each 
transducer.  The  DPG  prototype,  as  installed  at  the  CERC  Field 
Research  Facility,  utilizes  gin  and  1/8  inch  I.D.  flexible  tubing  in 
its  sensor  systems.  Ambient  temperature  changes  induce  differential 
pressure  changes  as  measured  by  the  transducers  since  the  tubing  used 
in  the  arm  sensors  is  four  times  the  length  of  that  of  the  center 
sensors . 

A  simple  bench-calibration  device  was  developed  to  monitor 
the  response  of  the  transducers  and  sensor  systems  to  applied  loads. 
It  is  important  that  the  load  be  distributed  uniformly  across  the 
sensor  diaphragm. 


Experiments  conducted  with  the  first  DPG  model  in  the  large 
wave  tank  at  the  Coastal  Engineering  Research  Center  were  plagued  by 
highly  non-linear  waves  and  the  failure  of  three  differential 
pressure  channels.  However,  results  from  the  absolute  pressure  gauge 
and  one  differential  pressure  gauge  indicated  the  ability  of  the  DPG 
to  estimate  wave  height,  period,  and  direction.  The  tests  were  also 
valuable  in  defining  several  important  modifications  required  to 
improve  the  performance  of  the  DPG  system. 

Measurements  of  the  orientation  of  the  instrument  were 
discussed  and  the  importance  of  careful  consideration  of  compass 
readings  in  the  vicinity  of  steel  structure  was  stressed.  Although 
it  initially  appeared  that  diving  compasses  were  placed  upon  the 
instrument  sufficiently  far  from  the  steel  cradle  to  avoid  its 
magnetic  effects,  redundant  measurements  of  the  instrument 
orientation  indicated  that  the  compass  readings  were  biased.  The 
error  was  resolved  using  a  simple  technique  outlined  in  the  thesis. 
It  was  suggested  that  measurements  from  simple  wrist  compasses,  when 
taken  and  interpreted  carefully,  may  be  adequate  in  determining 
instrument  orientation. 

The  DPG  utilizes  four  differential  pressure  gauges  placed  to 
the  bow,  port,  stern,  and  starboard  of  an  absolute  pressure  gauge 
such  that  there  are  two  adjacent  differential  gauges  along  each  of 
two  perpendicular  axes.  The  DPG  system  develops  the  first  five 
directional  Fourier  coefficients  using  signals  from  the  absolute 


pressure  gauge  and  one  differential  pressure  gauge  along  each  axis. 
The  other  differential  transducers  on  each  axis  were  installed  for 
redundancy  and  for  the  generation  of  higher  order  directional  Fourier 
coefficients.  The  absolute  pressure  transducer  and  one  differential 
pressure  transducer  along  each  axis  (dP2  and  dP3)  appear  to  work 
satisfactorily.  The  redundant  differential  pressure  gauges  dP4  and 
dPl  intermittently  exhibit  drift  and  high  frequency  oscillation, 
respectively.  The  failure  of  either  gauge  invalidates  reliable 
generation  of  higher  order  directional  Fourier  coefficients  by 
conventional  techniques  as  discussed  in  this  paper.  Accurate, 
reliable  generation  of  the  higher  order  (N*3)  directional  Fourier 
coefficients  through  the  subtraction  of  collinear  differential 
pressure  signals  —  similar  to  the  generation  of  lower  order 
directional  coefficients  through  the  subtraction  of  adjacent  absolute 
pressure  signals  —  remains  questionable. 

For  well-defined  swell,  directional  spectra  generated  from 
the  signals  of  dP2,  dP3,  and  the  absolute  gauge  correspond  well  with 
the  directional  spectra  generated  from  the  signals  of  dP4,  dP3,  and 
the  absolute  gauge.  There  is  less  agreement  for  poorly-defined  sea. 
Directional  estimates  made  without  the  absolute  pressure  gauge 
signals  are  generally  unreliable  using  techniques  discussed  in  this 


Estimates  of  wave  direction  developed  from  DPG  data  agree 
best  with  those  of  the  CERC  radar  for  a  17  minute  long  DPG  record 
with  2  Hertz  effective  sampling  and  16  degrees  of  freedom.  Analysis 
results  for  an  8.5  minute  record  with  4  Hertz  sampling  and  8  degrees 
of  freedom  differ  from  the  results  of  a  17  minute  record  depending 
upon  which  half  of  the  record  is  analyzed. 

Possible  modifications  to  the  instrument  design  and 
recommendations  for  future  DPG  development  have  been  outlined.  The 
configuration  of  the  DPG  arms  and  the  manner  in  which  the  data  is 
processed,  transmitted,  and  analyzed  might  be  altered  to  obtain  more 
accurate  or  effective  directional  wave  data. 

The  DPG  directional  wave  monitor  prototype  has  thus  far 
provided  reliable,  accurate  directional  wave  data  through  the  use  of 
a  smaller  and  lighter  instrument  than  conventional  pressure  sensor 
arrays.  It  generates  directional  wave  spectra  using  data  that  is 
directly  measured  instead  of  arithmetically  created  and  could 
potentially  develop  higher  order  directional  Fourier  coefficients 
than  conventional  arrays.  Overall,  the  DPG  appears  to  be  a  promising 
advancement  in  the  technology  of  directional  wave  measurement. 
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Mean*  2.381  S.Dev*  0.03698  Lower  Limit*  2.270  Upper  Limit*  2.492 

TIDE  LEVEL  IS  1.245  FEET. 

POSSIBLE  INSTRUMENT  TILT  (deg) 

DP2  DP3 

-0.048  -2.962 


Significant  Wave  Height  (ABS  gage)  *  4.05  feet. 


Effective  period  *  8.073  seconds. 


Maximum  Energy  located  at  bands: 

10.189  exceeds  average  energy  by  704.75  % 

9.438  exceeds  average  energy  by  448.02  % 

11.070  exceeds  average  energy  by  385.94  % 

8  790  exceeds  average  energy  by  196.96  % 

8.225  exceeds  average  energy  by  182.72  % 

7. 728  exceeds  average  energy  by  117.25  % 

Compass  heading  from  which  waves  propogate 


Band 

10. 189 

secs : 

Greatest 

Energy 

at 

65. 

deg 

Band 

9.438 

secs : 

Greatest 

Energy 

at 

61  . 

deg 

Band 

1 1 .070 

secs : 

Greatest 

Energy 

at 

67  . 

deg 

Band 

8.790 

secs : 

Greatest 

Energy 

at 

69  . 

deg 

Band 

8 . 225 

secs  : 

Greatest 

Energy 

at 

63. 

deg 

Band 

7  728 

secs : 

Greatest 

Energy 

at 

69  . 

deg 
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-  3 .14  secs  j 

-  11.07  secs 

-  8.73  secs 

-  8.23  secs 

-  7.73  secs 

i 
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UNIVERSITY  OF  DELAWARE 
DEPARTMENT  OF  CIVIL  ENGINEERING 
DPG  DIRECTIONAL  WAVE  MONITOR  /  FRF.  DUCK.  N.C. 


DEV 

YR 

MO 

DAY 

TIME 

1  1 1 

82 

6 

10 

13  0 

121 

82 

6 

10 

13  0 

141 

82 

6 

10 

13  0 

4096  POINTS  ANALY2ED  AT  .50  SECOND  SAMPLING 

8  BANDS  BLOCK  AVERAGED 

GAGE  NUMBER  2 

6  LOW  truncations;  55  HIGH  truncations 
Mean*  2.272  S.Dev*  0.21523  Lower  Limit*  1.626  Upper  Limit*  2.917 

GAGE  NUMBER  3 

3  LOW  truncations;  41  HIGH  truncations 
Mean*  2.116  S.Dev*  0.26262  Lower  Limit*  1.329  Upper  Limit*  2.904 

GAGE  NUMBER  5 

0  LOW  truncations;  2  HIGH  truncations 

Mean*  2.368  S.Dev*  0.05197  Lower  Limit*  2.212  Upper  Limit*  2.524 

TIDE  LEVEL  IS  0.893  FEET. 

POSSIBLE  INSTRUMENT  TILT  (beg) 

DP2  DP3 

-0.285  -2.871 


Significant  Wave  Height  (ABS  gage)  *  5.60  feet. 


Effective  period  *  8.412  seconds. 


Maximum  Energy  located  at  bands: 

11.070  exceeds  average  energy  by  1016.50  % 
9.438  exceeds  average  energy  by  377.97  % 
10.189  exceeds  average  energy  by  349.36  % 
12.118  exceeds  average  energy  by  248.42  % 
8.790  exceeds  average  energy  by  184.15  7. 
8 . 225  exceeds  average  energy  by  93 . 46  7. 

Compass  heading  from  which  waves  propogate 

Band  11.070  secs;  Greatest  Energy  at  85.  deg 
Band  9.438  secs:  Greatest  Energy  at  87.  deg 
Band  10.189  secs;  Greatest  Energy  at  83.  deg 
Band  12.118  secs:  Greatest  Energy  at  83.  deg 
Band  8.790  secs:  Greatest  Energy  at  91.  deg 
Band  8.225  secs:  Greatest  Energy  at  85.  deg 
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UNIVERSITY  OF  DELAWARE 
DEPARTMENT  OF  CIVIL  ENGINEERING 
DPG  DIRECTIONAL  WAVE  MONITOR  /  FRF.  DUCK.  N.C. 


DEV 

YR 

MO 

o 

J> 

< 

TIME 

111 

82 

6 

1 1 

7  0 

121 

82 

6 

1 1 

7  0 

141 

82 

6 

1 1 

7  0 

4096  POINTS  ANALYZED  AT  .50  SECOND  SAMPLING 

8  BANDS  BLOCK  AVERAGED 

GAGE  NUMBER  2 

8  LOW  truncations;  6  HIGH  truncations 

Mean*  2.299  S.Dev*  0.09652  Lower  Limit*  2.010  Upper  Limit*  2.589 

GAGE  NUMBER  3 

0  LOW  truncations;  8  HIGH  truncations 
Mean*  2.039  S.Dev*  0.19138  Lower  Limit*  1.465  Upper  Limit*  2.613 

GAGE  NUMBER  5 

0  LOW  truncations;  1  HIGH  truncations 

Mean*  2.364  S.Dev*  0.03900  Lower  Limit*  2.247  Upper  Limit*  2.481 

TIDE  LEVEL  IS  0.824  FEET. 

POSSIBLE  INSTRUMENT  TILT  (deg) 

DP2  OPS 

-0.061  -3.365 


Significant  Wave  Height  (ABS  gage)  «  4.32  feet. 


Effective  period  »  7.817  seconds. 


Maximum  Energy  located  at  bands: 

1 1 . 070  exceeds  average  energy  by  715.00% 
10.189  exceeds  average  energy  by  266.13  % 
12.118  exceeds  average  energy  by  194.16  % 
13.386  exceeds  average  energy  by  159.03  % 
8 . 790  exceeds  average  energy  by  115.97  % 
8.225  exceeds  average  energy  by  108.46  % 

Compass  heading  from  which  waves  propogate 

Band  11.070  secs:  Greatest  Energy  at  73.  deg 
Band  10.189  secs:  Greatest  Energy  at  71.  deg 
Band  12.118  secs:  Greatest  Energy  at  67.  deg 
Band  13.385  secs:  Greatest  Energy  at  69.  deg 
Band  8.790  secs:  Greatest  Energy  at  69.  deg 
Band  8.225  secs:  Greatest  Energy  at  75.  deg 


6.25  7.50  8.75  10.00 


S/ll/82  0700  EST 


-  11.0?  secs  j 

- . —  10.19  secs  ; 

• -  12.12  sec*  > 

-  12.12  sec*  j 

- -  8.79  secs  1 

-  8.29  secs  j 


! 

i 

i 

i 

J 


\ 


DIRECTION  (deg? 
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TABLE  A-l :  COMPARISON  OF  OTHER  DIRECTIONAL  ESTIMATES 
June  8-11,  1982 


DATE 

TIME 

PERIOD 

(secs) 

SPECTRAL 

PEAK 

(deg) 

ESTIMATE 

dP2*dP3 

(deg) 

MEAN  DIRECTION 

tan1 

(deg)  (deg) 

6/8 

0?00 

8.22 

59 

144 

60.1 

65.3 

8.79 

61 

148 

64.1 

65.7 

9.44 

61 

194 

60.0 

68.3 

1320 

8.79 

59 

238 

58.8 

62.1 

9.44 

61 

104 

61.0 

65.O 

8.23 

57 

188 

57-9 

62.1 

6/9 

0700 

9-44 

67 

109 

66.7 

79.9 

10.19 

63 

144 

63.2 

73.8 

8.79 

63 

113 

61.5 

74.8 

6/10 

0700 

10.19 

65 

119 

65.I 

74.5 

9.44 

61 

194 

61.1 

65.9 

11.07 

67 

276 

68.6 

7  8.4 

1300 

11.07 

85 

85 

84.6 

117.7 

9.44 

87 

85 

86.6 

120.5 

10.19 

83 

200 

82.4 

110.9 

6/11 

0700 

11.0? 

73 

76 

74.8 

94.2 

10.19 

71 

204 

68.2 

88.0 

12.12 

67 

227 

73.0 

76.8 

APPENDIX  B 


TYPICAL  STRIP  CHART  RECORDINGS  FROM 
BENCH  CALIBRATION 

Results  from  Calibration  Test 
Performed  1  Week  before  Field  Installation 


ill  1 1  ill  min  ii  It  i  ii  tiiiliiiiliiiiliiiiliiH  iiiiIiii»IiiiiIiiiiIiiii|iihIi»iiIuiiIiiiiIiiiiJiii|Ii»iiIiiiiIh,*I 


Chamber  Pressure  Xdcr. 

Cal.  Factor  (top}:  .05  psi/mm  Calib.  F 
Cal.  Factor  (bottom) s  0.1  psi/mm 

Chart  Speed:  5  mm/sec 


Calib.  Factor:  .02  psid/mm 
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APPENDIX  D 


ASSEMBLY  AND  MAINTENANCE  OF  THE  DPG 


INSTRUMENT  ASSEMBLY  AND  MAINTENANCE 


Initial  Assembly 

1.  The  transducer  stack  is  created  by  placing  the  transducers 
within  the  plexiglass  wafers  and  attaching  the  acrylic  stand-off 
rods  between  the  wafers. 

2.  The  stack  is  attached  to  the  bottom  of  the  water-tight 
cylinder  (can)  end-cap. 

3.  Nylon  tubing  is  connected  between  the  the  bottom  of  the 
end-cap  and  the  transducer  pressure  ports. 

4.  The  wiring  is  then  completed  terminating  to  a  16-pin  connector 
positioned  just  below  the  can  end-cap. 

5.  The  exterior  stainless  steel  tubing  and  acrylic  diaphragm 
housings  are  attached  to  the  top  of  the  end-cap. 

6.  The  center  housings  are  larger  in  diameter  and  should  connect 
to  the  slightly  "elevated"  Swage-Locks  on  the  transducers. 


Fluid  Back-Filling 


6.  The  center  differential  lines  are  filled  first.  The 
transducer  stack  is  laid  on  its  side  with  the  bottom  slightly 
elevated.  Two  center  bleed  ports  face  up.  The  line  to  be  filled 
is  flexed  up  so  that  the  housing  is  above  the  stack.  Fluid  is 
poured  in  to  the  housing  until  a  steady  flow  streams  from  the 
bleed  port.  The  port  is  then  closed.  (Figure  D~l.) 

7.  Enough  fluid  must  be  poured  into  the  housing  so  that  it  forms 
an  inverted  meniscus  about  the  o-ring.  The  elastomer  diaphragm  is 
slid  horizontally  over  the  fluid  until  its  holes  line  up  with  the 
threaded  holes  in  the  housing.  The  copper-nickel  ring  is  then 
lowered  carefully  onto  the  elastomer  and  secured  to  the  housing 
with  Monel  screws.  Each  housing,  elastomer,  and  ring  is  labelled 
as  a  set.  The  labels  should  fall  in  line  so  that  the  holes  in 
each  of  the  three  match. 

8.  The  process  is  repeated  for  each  of  the  center  differential 
1 ines . 

9.  The  arm  sensors  are  filled  next. 

10.  A  large  plate  is  attached  to  the  bottom  of  the  stack  using 
the  threaded  (  6~32)  holes  in  the  stand-offs. 

11.  The  stack  is  inverted  (end-cap  down)  on  top  of  a  stand  or 


FIGURE  D-2:  Back-filling  the  arm  sensor*.  (Brian 
Jacobs  shown  with  instrument.) 


platform. 


12.  The  arm  line  to  be  filled  is  lifted  above  its  transducer  and 
the  filling  process  is  repeated.  (Figure  D~2.) 

13.  After  each  arm  line  is  filled,  the  stack  is  positioned 
upright  to  fill  the  absolute  transducer. 

14.  The  absolute  transducer  has  no  bleed  port,  so  the  nylon 
tubing  is  removed  from  the  transducer  and  the  port  filled  to 
overflowing.  Fluid  is  then  poured  into  the  absolute  line's  sensor 
housing  until  it  streams  out  of  the  bottom  of  the  nylon  line. 
While  it  is  still  flowing,  the  line  is  inserted  into  the  pressure 
port  and  secured.  The  housing  is  then  sealed  with  the  elastomer 
diaphragm  ana  ring  as  usual. 

15.  The  exterior  tubing  should  be  flexed  and  the  assembly 
re-positioned  and  inspected  several  times  to  ensure  that  no  air 
has  been  trapped  in  the  lines.  It  should  then  be  placed  in 
sunlight  or  a  warm  temperature  —  with  the  isolation  diaphragms 
above  the  transducers  —  to  allow  any  small  bubbles  in  the  fluid 
to  appear. 

16.  Air  bubbles  that  appear  can  be  eliminated  by  removing  the 
copper-nickel  ring  and  elastomer,  adding  more  fluid  to  the  housing 
and  then  re-sealing. 


17.  Once  the  system  is  air-free,  the  instrument  is  ready  for 


calibration. 


C.  Assembly  of  the  Housing 

18.  After  calibration,  the  can  is  slid  over  the  stack  —  guiding 
the  cable  from  the  bottom  of  the  can  into  a  large  cut-out  on  the 
outer  edge  of  the  wafers. 

19.  When  the  can's  cable  reaches  the  transducer  stack's 
electrical  connector,  the  two  are  mated  and  the  can  is  secured  to 
the  end-cap. 

20.  The  can  is  then  slid  into  the  bottom  of  the  fuselage,  sensors 
first.  The  arm  sensors  are  guided  out  the  four  large  fuselage 
holes  and  the  can  pushed  upwards  until  it  hits  the  stops  inside 
the  fuselage. 

21.  The  bottom  stops  are  then  screwed  into  the  fuselage  walls  to 
secure  the  can. 

22.  This  assembly  is  then  stood  up  on  blocks  (to  protect  the 
titanium  bolts  on  the  underside  of  the  fuselage). 

23.  The  PVC  spider  is  pushed  in  through  the  top  of  the  fuselage 
until  it  lines  up  with  the  arm-holes.  The  PVC  arms  are  pushed 
through  the  fuselage  walls  and  into  the  spider.  They  are  secured 
by  bolts  that  thread  through  the  top  of  the  spider. 
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24.  The  sensors  are  secured  to  the  inside  of  the  PVC  arms  by  two 
brass  screws  that  pass  through  the  arms  and  thread  into  the 
acrylic  housings. 

25.  The  center  plate  is  next  placed  into  the  top  of  the  fuselage 
and  the  center  sensors  positioned  within  it. 

26. '  The  plexiglass  end-cap  is  then  attached  to  the  top  of  the 
fuselage. 

D.  Field  Preparation  and  Installation 

27.  Copper  scour  pads  are  placed  into  the  ends  of  the  arms  and 
the  arm  end-caps  screwed  on. 

28.  The  cable  is  retreived  and  mated  to  the  instrument  on  the 
underside  of  the  can  and  a  plexiglass  plate  attached  to  the  bottom 
of  the  fuselage. 

29.  With  the  cradle  already  on  the  seafloor,  the  instrument  is 
brought  underwater  and  positioned  onto  it.  The  plate  breaching 
the  gap  about  the  center  of  the  cradle  must  be  removed  to  allow 
passage  of  the  cable,  then  replaced  after  installation. 

30.  Free  cable  may  then  be  stored  in  the  service  box;  the  top  of 
which  is  removed  by  unthreading  a  titanium  nut  on  each  top  corner 
of  the  box. 


Bergman  (1969)  has  developed  a  technique  to  estimate  higher 
order  directional  Fourier  coefficients  which  does  not  require  the 
subtraction  of  signals  from  adjacent  pressure  gauges.  If  the 
distance  and  angle  between  a  pair  of  pressure  gauges  in  an  array  is  D 
and  ,  respectively,  then  the  cross-spectral  density,  normalized  by 
the  unidirectional  spectral  density  S(6)  can  be  expressed  as: 


CO 

c'+iq"  =  J  (kD)  +  2  Y  (a  cos  ng  +  b  sin  ng)ir(i)nJ  (kD) 
o  L‘ ,  n  n  n 

n=l 


(F.l) 


for  some  frequency  c3  •  Jn  (argument)  and  Jr1(argument)  are  Bessel 
functions.  In  terms  of  real  and  imaginary  parts,  Eq.  F.l  can  be 


>-• 


-  202  - 


c'  =  JQ(kD) 

-  2ir  J2(kD)(a2  cos  26  +  b2  sin  26) 

+  2 it  J^CkD)  (a4  cos  46  +  b4  sin  46) 

-  2n  J/.(kD)  (a,  cos  66  +  b,  sin  66) 

6  o  o 

+  ... 

q'  =  2ir  (kD)  (a1  cos6  +  b1  sin  6) 

-  2ir  (kD)  (a3  cos  36  +  b3  sin  36) 

+  2ir  J,.  (kD)  (a^  cos  56  +  b5  sin  56) 


The  Fourier  series  can  be  accurately  truncated  to  N  if  J  (kD)  is 
negl igibte  for  n>N  or  if  an  and  bn  are  negligible  for  n>N.  There  is  a 
set,  then,  of  2N  equations;  N  equations  from  each  of  the  real  and 
imaginary  parts  of  Equation  F.l,  and  each  with  the  unknowns  an  and 
bn •  The  maximum  possible  value  of  N  is  determined  by  the  number  of 
different  non-ambiguous  cross-correlations  possible  from  the  gauges 
in  the  array.  If  there  are  no  ambiguous  cross-correlations  between 
the  gauges,  the  number  of  possible  combinations  for  M  gauges  is 


M (M— 1) /2 .  This  is  the  number  of  coefficient  pairs  N  that  one  can 
develop.  It  is  suggested  that  one  does  not  develop  all  N  pairs,  but 
stops  at  least  two  pairs  sooner  so  that  there  are  more  equations  than 
unknowns.  The  coefficients  may  then  be  fitted  by  least  squares 
analysis . 

As  an  example,  four  pressure  gauges  positioned  at  the  axes  of 
an  irregular  three-faced  pyramid  could  potentially  develop  six 
coefficient  pairs  (or  the  first  thirteen  directional  Fourier 
coefficients).  It  could  more  effectively  develop  four  coefficient 
pairs  (or  the  first  nine  coefficients)  using  least  squares  analysis 
to  fit  the  coefficients.  Signals  from  four  pressure  gauges 
positioned  at  the  corners  of  a  square,  such  as  in  the  Scripps  Sxy 
gauge  (Seymour , 1978) ,  develop  only  four  non- ambiguous 
cross-correlations  and  thusly  could  generate  the  first  four 
coefficient  pairs  (or  the  first  nine  coefficients).  There  are  two 
other  possible  (ambiguous)  correlations  that  can  be  utilized  to  fit 
the  first  nine  coefficients  by  least  squares  analysis  —  just  as  the 
pyramid-shaped  array  does  by  excluding  the  development  of  the  last 
two  coefficient  pairs. 

Since  the  DPG  utilizes  differential  pressure  gauges,  the 
relations  for  the  cross-spectral  densities  differ  from  Equation  F.l. 
The  cross-  spectral  density  for  two  colinear  differential  pressure 
signals  can  be  shown  to  be: 
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00050  C 

OOIOO  c 
00150  C 
00200  C 
00250  C 
00300  C 
00350  C 
00400  C 
00450  C 
00500  C 
00550  C 
00600  C 
00650  C 
00700  C 
00750  C 
00800  C 
00850  C 
00900  C 
00950 
01000 
01050 
Ol  100 
01  150 
01200 
01250 
01300 
01350 
01400 
01450 
01500 
01550 
01600 
01650 
01700 
01750 
01800 
01850 
01900 
01950 
02000 
02050 
02100 
02150 
02200 
02250 
02300  C 
02350  C 
02400  C 
02450  C 
02500  C 
02550  C 
02600  C 
02650  C 
02700  C 
02750  C 
02800  C 
02850  C 
02900  C 
02950 
03000 


This  program  is  designed  for  use  with  the  University 
of  Delaware  DPG  Directional  Wave  Monitor  as  Installed 
at  the  FRF  on  14  May.  1982.  The  instrument  and  software 
were  designed  by  Kevin  R.  Bodge  under  the  direction  of 
Dr.  Robert  G.  Dean. 

The  program,  as  utilized  by  the  FRF.  uses  only 
differential  pressure  channels  dP2  and  dP3  with  the 
absolute  gauge.  It  generates  the  directional  spectra 
using  the  first  5  directional  Fourier  coefficients. 

The  program  as  listed  utilizes  all  four  differ¬ 
ential  gauges  with  the  absolute  gauge  and  Is  capable 
of  generating  the  directional  spectra  using 
the  first  seven  directional  Fourier  coefficients. 

Fngl Ish  units  are  used  throughout. 


1 

1 

1 


1 

1 


DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

INTEGER 


dpi (4 100) , dp2( 4 100 ) , dp3( 4100) ,dp4( 4 100) . 

abs(4100) . dpxx(4100) ,dpyy(4100) 

d 1 1 ( 2050 ) . d2 1 ( 2050 ) , d3 1 ( 2050 ) . d4 1 ( 2050 ) . 

ab 1 ( 2050 ) . dxx 1 ( 2050 ) . dy y 1 ( 2050 ) 

d  1  (  2050 )  ,  d2  (  2050  )  .  d3  (  2050 )  .  d4  <  2050 ) . 

ab( 2050 ) . dxx ( 2050 ) , dyy ( 2050 ) 

prsp( 2050) . wvnr( 2050) . f req( 2050) . maxf rq( 20) 

thta 1(400 ) , thta2( 400 ) , thta3(400) ,pd( 20) 

aO( 400) , a  1 ( 400 ) , a2( 400 ) . a3( 400) . 

b 1 ( 400 ) . b2 ( 400 ) . b3 ( 400 ) 

dl r 1 ( 400) . d I r2( 400) ,d 1 r3( 400 ) , dl r4( 400) . 
ang 1 ( 400) , ang2( 400) ,ang3( 400) . ang4(400) 
blk 


OPEN (uni t*  18 , device* 'dsk ' . f 1 1e= 'DP  1 . dat ' ) 
OPEN( un 1 t= 19, device* 'dsk ' , f 1 1 e* ' DP2 . dat ' ) 
0PEN(un1 t= 20, device* 'dsk ' .file* 'DP3 . dat ' ) 
0PEN(un1t=21, device* ' dsk ' , f 1 le= 'DP4 . dat ' ) 
OPEN( unit  =  22. dev  ice* 'dsk' . f 1 le= ' DP5 . dat ' ) 
0PEN(un1 t = 1 7 , device* 'dsk ' , f 1 1 e* ' PWR . DAT ' ) 
0PEN(un1t=14, device* ' dsk ' .file*' RAW . dat ' ) 
OPEN( un 1 1  =  3. dev  Ice* 'dsk' . f 1 le* 'DIR . dat ' ) 
OPE N( uni t*  4 .device* 'dsk' ,f11e='FRF  dat' ) 


INPUT  ANALYSIS  DETAILS: 

ATM  *  Atmospheric  Pressure  (mb) 

GAMMA  *  Seawater  Specific  Gravity  (Ib/cu.ft.) 

N  ■  I  of  sampled  points  to  analyze  (power  of  2 ) 

K  *  such  that  2* *K  =  N 

dT  *  sampling  Interval  (secs) 

BLK  *  0  of  freq.  bands  to  b 1 ock -average 
NUFF  *  M  of  peak  energy  bands  to  identify 
Icut  *  w  of  non-b 1 ock -averaged  cut-off  freq. 


ATM  *  1013.0 

ATM  *  ATM»0. 0145038 


-  209  - 


03050 
03100 
03  150 
03200 
03250 
03300 
03350 
03400 
03450 
03500 
03550 
03600 
03650 
03700 
03750 
03800 
03850 
03900 
03950 
04000 
04050 
04100 
04150 
04200 
04250 


05300 

05350 

05400 

05450 

05500 

05550 

05600 

05650 

05700 

05750 

05800 

05850 

05900 

05950 

06000 


GAMMA  =  64  18 
GAMMA  =  GAMMA/ 1728. 
N  =  4096 
K=  1  1 

N2  =  N/2 
DT  =  0.5 
MOST  -  N+4 
BLK  =  8 
NUFF  =  6 
Icut  =  327 


OPG  PHYSICAL  DETAILS: 

ORIENT  «  compass  heading  of  arm  3  w.r.t.  true  north 


SGAGfc 
OX  Ini 
DX1n2 
DX1n3 
DX  1n4 
ARM  13 
ARM24 


mean  depth  of  sensors  below  mean  water  (ft) 
length  of  arm  1  dlff.  measurement  (In) 


length  of  arm 
length  of  arm 
length  of  arm 


measurement  (In) 
measurement  (In) 
measurement  (In) 


distance  between  DPI  &  DPS  measurement  (in) 
distance  between  DP2  8  0P4  measurement  (In) 


Pbf  =  fluid  back-fill  pressure  on  absolute  (psl) 


04300 

ORIENT  *  53. 

04350 

SGAGE  =  16.49 

04400 

DXInl  =  47  25 

04450 

DX1n2  *  53.75 

04500 

DX1n3  =  47.75 

04550 

DX1n4  *  53.875 

04600 

ARM  1 3  =  53.25 

04650 

ARM24  *  50.94 

04700 

Pbf  *  1.23 

04750 

04800 

04850 

pi  =  3. 14159265 

04900 

twopl  *  2.0  *  pi 

04950 

05000 

6008 

Format ( 14 ) 

05050 

WRITE (4,819) 

05100 

819 

Format(//,40X, 'UNIVERSITY  OF 

05150 

1 

36X, 'DEPARTMENT  OF  CIVIL 

05200 

2 

29X , 'DPG  DIRECTIONAL  WAVE 

05250 

3 

'FRF.  DUCK,  N.C. ' ,//) 

Load  data  arrays  with  zeroes  for  cleanliness. 
DO  115  1*1, most 
DP  1  ( 1 )  *  0.00 
DP2( 1 )  «  0.00 
DP3( 1  )  «  0.00 
DP4(  1  )  -  0.00 
15  ABS( 1 )  «  0.00 


3055 


CALL  INPUT  (dpi. dp2 . dp3 , dp4 , abs ) 

WRI TE ( 4 , 3055 )  N.dT 

Format (/, 35x , 15, '  POINTS  ANALYZED  AT  ',f4.2,'  SECOND', 


(8050 

18100 

18150 

18200 

18250 

18300 

18350 

13400 

18450 

18500 

18550 

18600 

18650 

18700 

18750 

18800 

18850 

18900 

18950 

19000 

19050 

19100 

19150 

19200 

19250 

19300 

19350 

19400 

19450 

19500 

19550 

19600 

19650 

19700 

19750 

19800 

19850 

19900 

19950 

20000 

20050 

20100 

20150 

20200 

20250 

20300 

20350 

20400 

20450 

20500 

20550 

20600 

20650 

20700 

20750 

20800 

20850 

20900 

20950 

21000 


READ! 19.900)  dev . yr , mo , day , hr . m 1 n 
WRITE(4,901 )  dev, yr, mo, day. hr, min 
n=  1 

00  202  J  =  1 , 4 10 

RE AD( 19 , 905 )  (A< 1 ) . 1= 1 . 10) .LAST 
00  102  1-1,10 
0PV2(n)  *  A( 1 )/200 . 0 

102  n=n+ 1 

202  Continue 

READ( 20, 900)  dev . yr . mo , day , hr , m 1 n 
WRITE (4, 901 )  dev.yr , mo . day , hr , min 
n=  1 

DO  203  j-1.410 

RE AD( 20, 905 )  (  A  (  1 ) .  1  =  1 . 10 ) . LAST 
DO  103  1*1.10 
DPV3(n)  =  A( 1 )/200 . 0 

103  n=n+1 

203  Continue 

READ! 2 1 , 900 )  dev , yr , mo , day , hr , ml n 
WRITE (4 ,901 )  dev.yr . mo . day . hr , min 
n=  1 

DO  204  j  *  1 , 4  10 

RE AD( 2 1 , 905 )  ( A(  1  )  .  1  =  1 . 10) , LAST 
DO  104  1=1,10 
DPV4(n)  *  A( 1 )/200 . 0 

104  n=n+1 

204  Continue 

READ( 22 ,900)  dev , yr , mo , day , hr , m 1 n 
WRITEI4. 901)  dev , yr , mo , day , hr , m 1 n 
n»  1 

DO  205  J« 1,410 

RE AD( 22 , 905 )  ( A( 1 ) , 1  -  1 . 10) , LAST 
00  105  1=1,10 
ABSV(n)  =  A ( 1 )/200 . 0 

105  n=n+1 

205  Continue 


890  Format (//, 15X. '  DEV  YR  MO  DAY  TIME') 

900  Format! 614) 

901  Format ( 15x , 6 1 4 ) 

905  Format! 1015, 1 1 ) 

Return 

End 


SUBROUTINE  PERIOD  ( ab , dl , d2 , d3 , d4 , n, lout ,dT , prsp, wvnr , 
1  freq) 

C  Computes  the  centroid  of  the  power  spectra 
C  and  reports  the  corresponding  frequency 

C  as  the  "effective  period*. 

DIMENSION  ab(2050) ,d1( 2050) ,d2( 2050) ,d3( 2050) ,d4( 2050) 
1  prsp( 2050) , wvnr ( 2050) , f req( 2050) 


-  215  - 


21050 
21100 
21150 
2  1200 
21250 
21300 
21350 
21400 
21450 
21500 
21550 
21600 
21650 
21700 
21750 
21800 
21850 
21900 
21950 
22000 
22050 
22100 
22150 
22200 
22250 
22300 
22350 
22400 
22450 
22500 
22550 
22600 
22650 
22700 
22750 
22800 
22850 
22900 
22950 
23000 
23050 
23100 
23150 
23200 
23250 
23300 
23350 
23400 
23450 
23500 
23550 
23600 
23650 
23700 
23750 
23800 
23850 
23900 
23950 
24000 


twopt  -  6.283185 
N2  =  N/2 

DATA  sumabk , sumlk , sum2k , sum3k , sum4k . sabsk . slk . s2k . 
s3k,s4k/0. .0. .0. .0. .0. .0. .0. .0. .0. .0./ 

DO  420  J*2. Icut 

SUMABK  *  sumabk  +  ab(J) 

SUM IK  -  sumlk  +  dl(J) 

SUM2K  *  sum2k  +  d2(j) 

SUM3K  >  sum3k  +  d3(J) 

SUM4K  »  SUM4K  +  d4  (  j  ) 

SABSK  «  SABSK  +  AB( J )  *  freq(J) 

S IK  »  slk  ♦  d 1 ( J )  *  freq( j ) 

32K  =  s2k  +  d2(j)  *  freq(J) 

S3K  =  s3K  +  d3( J )  *  freq(J) 

S4K  «  s4k  +  d4( j )  *  freq(j) 

Cent inue 


Tabsk  >  twopl  *  sumabk/sabsk 
Tlk  *  twopl  *  sumlk/stk 
T2k  •  twopl  *  sum2k/s2k 
T3k  «  twopl  *  sum3k/s3k 
T4k  »  twopl  *  sum4k/s4k 

C  Report  the  effective  period  from  the  absolute  gage. 
WR I TE ( 4 , 929 )  Tabsk 

929  Format( '  Effective  period  *'.f7.3.'  seconds.') 

Return 

End 


SUBROUTINE  SURF  ( Icut , dl ,d2 ,d3 . d4 , dpi . dp2 , dp3. dp4 . 

1  d11,d21.d31,d41,d1r1,d1r2.d1r3.d1r4.blk, 

2  orient) 

DIMENSION  d 1 ( 2050 ).d2(2050),d3( 2050 ) . 

1  d4 ( 2050 ) 

DIMENSION  dp1(4100),dp2(4100), dp3( 4 10O) , dp4( 4 100) , 
1  d 1 1 ( 2050 ) . d2 1 ( 2050 ) , d3 1 ( 2050 ) . d4 1 ( 2050 ) 

DIMENSION  dir  1 ( 400) td1r2(400) , d1r3( 400) ,d1r4(400) 
INTEGER  blk 


Generates  directional  estimate  without  absolute  gage. 
0IR1  uses  gages  DPI  and  DP2; 

DIR2  uses  gages  0P3  and  DP4; 

DIR3  uses  gages  DP2  and  0P3; 

0IR4  uses  gages  DPI  and  0P4 . 

00  435  <J«2 , tcut 

1 f (dpi ( j ) . ne.0 .0) tan1«atan2(d1 1  ( J )  .dpi  (  J ) ) 

If (dp1( j ) .eq.O. )tanl  ■  0.0 

If (dp2( J ) . na. O. O )tan2*a tan2 (d21 (J ) ,dp2( J ) ) 

If (dp2( J ) .eq.O. )tan2  ■  0.0 

1f(dp3( j) .ne.0  0)tan3  *atan2(d3t ( J ) ,dp3( J ) ) 

If (dp3(J ) .eq.O. )tan3  ■  0.0 


212 


12050 

DY V ( 1 )  =  2.0*(dpyy< i )*«2  +  dyy 1 ( 1 ) *  *  2 ) 

12100 

67 

Cont 1 nue 

12150 

12200 

C  Store  the  raw  power  spectra.  If  desired 

12250 

GO  TO  4003 

12300 

WRI TE ( 17,4007) 

12350 

WRITE( 17,3098 ) 

12400 

00  4003  1=2. 1  cut 

12450 

WRITE( 1 7 , 9 1 4 ) f  req( 1 ) . ab( 1 ) , dl ( 1 ) , d2 ( 1 ) . d3( 1 ) , 

12500 

1 

d4( 1 ) ,  dxx( 1 ) ,dyy( 1 ) 

12550 

4003 

Cont 1 nue 

12600 

4007 

Format( '  RAW  POWER  SPECTRA') 

12650 

12700 

12750 

C  Calculate  significant  wave  height  from  absolute  gage. 

12800 

E  =0 . 0 

12850 

DO  405  j  =2 , 1  cut 

12900 

405 

E  =  ab( j )  +  E 

12950 

Hslg  =  (4 . *sqrt(E ) )/l2 .0 

13000 

WRITE ( 4 . 92 1 )  Hslg 

13050 

921 

Format(/,'  Significant  Wave  Height  (ABS  gage)  =',f6.: 

13100 

1 

'  feet. ' .//) 

13150 

13200 

13250 

13300 

C  Determine  effective  period  from  all  gauges. 

13350 

13400 

CALL  PERIOD  ( ab . d 1 . d2 . d3 . d4 . n , i cut , dT . pr sp . wvnr , f req 

13450 

13500 

13550 

C  Estimate  direction  without  absolute  gage. 

13600 

13650 

GO  TO  6002 

13700 

CALL  SURF  ( icut.dl ,d2.d3,d4.dp1 , dp2 , dp3 , dp4 . 

13750 

1 

dl 1 ,d21 ,d31 ,d41 .dir  1 ,d1r2 ,d1r3.d1r4 ,blk , 

13800 

2 

orient ) 

13850 

13900 

13950 

CALL  SURF2(ab,d1 , d2 . d3 , d4 , angl , ang2 , ang3 , ang4 , icut . 

14000 

1 

wvnr , bl k , or  lent ) 

14050 

14100 

6002 

Cont Inue 

14150 

14200 

14250 

C  Generate  the  directional  Fourier  coefficients 

14300 

14350 

CALL  FCOEFS(ab,abs,abi ,d1 .dpi ,d1 1 , d2 . dp2 . d2 1 ,d3. 

14400 

1 

dp3 , d3 1 . d4 , dp4 . d4 1 ,dxx , dpxx ,dxx 1 , dyy , dpyy . dyy 1 , 

14450 

2 

wvnr ,freq,thta1,thta2,thta3. 

14500 

3 

aO. a  1 , a2 . a3. bl . b2 . b3. blk ,  icut) 

14550 

14600 

14650 

C  Store  the  b lock-averged  power  spectra 

14700 

WRITE <  17.3097 )  blk 

14750 

WRITE( 17.3098) 

14800 

3097 

Format(/,'  BLOCK  AVERAGED  BY'. 13. 

14850 

1 

'  BANOS') 

14900 

3098 

Format('  freq  .  ab  .  dl  .  d2  d3  .  d4'. 

14950 

1 

'  .  dxx  .  dyy ' . / ) 

15000 

DO  77  1=2. icut 

1 


15050 

15100 

15150 

15200 

15250 

15300 

15350 

15400 

15450 

15500 

15550 

15600 

15650 

15700 

15750 

15800 

15850 

15900 

15950 

16000 

16050 

16100 

16150 

16200 

16250 

16300 

16350 

16400 

16450 

16500 

16550 

16600 

16650 

16700 

16750 

16800 

16850 

16900 


WRITE< 17,914)  freq(l).  ab( 1 ) . d 1 ( i ) , d2( 1 ) . d3( I ) , 
d4( 1 ) ,dxx( 1 ) ,dyy( 1 ) 

Cont i nue 

Format ( f 11.3,f11.5.4F14.9,2F2018) 


C  Identify  bands  of  highest  energy 

Cal  1  MAXS I G( ab , 1  cut , maxf rq , n , dT , nuf f . f req) 

C  Generate  directional  spectra  for  highest  energ”  bands 

can  DIRXSP  ( AO, A 1.A2.A3.B1 , B2 , 83 ,  f req . n , dT , nuf f , 

1  maxf rq, or  lent ) 

C  GO  TO  9999 

WR ITE(4,910) 

DO  666  nn=1,nuff 
LKM  =  maxfrq(nn) 

PD(nn)  =  twop l/f req( 1  km) 

I F ( thta 1 ( lkm) . 1 t .0. ) thta1( 1  km)  =  360 . +thta  1  (  1km) 

IF( thta2( 1km) . 1 1 .0. )thta2( 1  km)  =  360 . + thta2 ( 1  km) 

IF( thta3( 1km) . 1 1  O.  )thta3( 1  km )  =  360  +thta3( 1km) 

WRITE (4 ,913)  PO(nn) , thta1( 1km) . thta2( 1km) , 

1  thta3( lkm) .dir  1 ( 1km) .dl r2( 1km) , d1r3( 1  km) . dl r4(  1km) 

2  angl ( lkm) ,ang2( 1km) , ang3< lkm) , ang4( lkm) 

666  Continue 

910  Format ( / , 5x . '  MEAN  DIRECTION  OF  PROROGATION' . 6x ESTIM ' 

1  'ATE  W/OUT  ABSOLUTE  GAGE : ' . 1 2x , ' WI TH  1  CHANNEL', 

2  '  AND  ABSOLUTE',/,'  Per  1 od ' . 5x , ' n» 1 ' , 7x , ' n=2 ' , 

3  7x, 'n»3' ,6x, ' 1-2' ,6x. '3-4' ,6x, '2-3' ,6x. ' 1-4' , 

4  1 1 x , ' dp  1 ' , 6x , ' dp2 ' , 6x , ' dp3 ' , 6x , ' dp4 ' ) 

913  Format( f6.2,3f10.2,4f9. 1,5x.4f9. 1) 


16950 

17000 

9999 

END 

17050  C 

17  100 

17150 

SUBROUTINE  INPUT(dpv1 , dpv2 , dpv3 , dpv4 . absv ) 

17200 

17250 

DIMENSION  dpv 1(4100), dp V2(4100) ,dpv3(4100) 

17300 

1 

dpv4(4 100) , absv( 4 100) , M( 4 100) 

17350 

INTEGER  dev . yr , day , hr , A( 10) 

17400 

17450 

WRITE<4 .890) 

17500 

17550 

REAO( 18,900)  dev , yr , mo , day , hr , m 1 n 

17600 

WRI TE( 4 , 901 )  dev , yr , mo , day , hr , ml n 

17650 

n=  1 

17700 

00  200  J  =  1 ,4  10 

17750 

REAO( 1 8 , 905 )  (A( 1 ), 1=1. 10). LAST 

17800 

DO  lOO  1*1,10 

17850 

DPV 1 ( n)  *  A ( 1 1/265.0  *  1.65 

17900 

100 

n=n+  1 

17950 

200 

Cont 1 nue 

10000 

-  210  - 


06050 

06100 

06150 

06200 

06250 

06300 

06350 

06400 

06450 

06500 

06550 

06600 

06650 

06700 

06750 

06800 

06850 

06900 

06950 

07000 

07050 

07100 

07150 

07200 

07250 

07300 

07350 

07400 

07450 

07500 

07550 

07600 

07650 

07700 

07750 

07800 

07850 

07900 

07950 

08000 

08050 

08100 

08150 

08200 

08250 

08300 

08350 

08400 

08450 

08500 

08550 

08600 

08650 

08700 

08750 

08800 

08850 

08900 

08950 

09000 


'  SAMPLING') 

WRITE ( 4 , 3056 )  blk 

Format ( 57x , 12 , '  BANDS  BLOCK  AVERAGED') 


C  Scan  each  record,  calculate  std. deviation,  truncate 
C  unreasonable  points,  report  them,  calculate  and  subtract 
C  mean  from  record,  convert  to  psl. 

NGAGE  »  1 

Cal  1  SCAN(dp1 , n , ngage , avgl .dxlnl ) 

NGAGE  *  2 

Call  SCAN( dp2 , n , ngage , a vg2 , dx 1 n2 ) 

NGAGE  -  3 

Cal  1  SCAN(dp3,n, ngage, avg3,dx1n3) 

NGAGE  =4 

Ca 1 1  SC AN( dp 4 , n , ngage , a vg4 , dx 1 n4 ) 

NGAGE  =5 

Cal  1  SCAN( abs , n , ngage , avgabs , dx 1 n4 ) 

N  *  N/2 
N2  *  N2/2 

C  Reverse  signs  on  DPI  &  DP2  for  sign  convention. 

DO  1057  JJ-1.N 

DPI(jj)  =  -DPI(jj) 

DP2 ( J  J )  =  -0P2 ( J  j ) 

1057  Continue 


Record  portion  of  pressure  signals  In  scratch  file 
If  desired. 

GO  TO  1056 
DO  1056  1-1,600 

WRITE( 14.4096)1 .  dpt( 1 ) ,dp2( 1 ) ,dp3( 1 ) ,dp4( 1 ) ,abs( 1 ) 
56  Cont 1 nue 
36  Format( 16. 5f 14.7) 


Calculate  tide  and  water  depth.  Neglect 

temperature  changes  on  back-filling  fluid. 
DEPTH  -  ((AVGABS* 10. -Pbf-atm)/gamma  )/12. 
TIDE  -  (DEPTH  -  SGA3E ) 

17  Formate/.'  TIDE  LEVEL  IS',F7.3,'  FEET.') 
WRITE ( 4,717)  TIDE 

Compare  record  mean  and  no-wave  condition  tare 
to  approximate  Instrument  tilt. 

TARE  1  -7.5 
TARE2  -  2.31 
TARE3  -2.56 
TARE4  -  3.825 

tl  «  (avgl  -  tare1)/5. 

t2  -  ( avg2  -  tare2)/5. 

t3  »  (avg3  -  tare3 )/5. 

t4  »  (avg4  -  tare4)/5. 

gl  «  gamma  *  dxlnl 

g2  •  gamma  *  dxin2 


09050 
09100 
09150 
09200 
09250 
09300 
09350 
09400 
09450 
09500 
09550 
09600 
09650 
09700 
09750 
09800 
09850 
09900 
09950 
10000 
10050 
10100 
10150 
10200 
10250 
10300 
10350 
10400 
10450 
10500 
10550 
10600 
10650 
10700 
10750 
10800 
10850 
10900 
10950 
1  1000 
11050 
1  1  100 
11150 
1  1200 
1  1250 
11300 
1  1350 
1  1400 
1 1450 
1  1500 
11550 
1  1600 
11650 
11700 
11750 
11800 
1 1850 
1 1900 
11950 
12000 


g3  =  gamma  *  dxin3 
g4  -  gamma  *  dx1n4 

TILT  1  -  atan2( 1 1 ,g1 )  *  360./twop1 
TILT2  *  atan2(t2.g2)  *  360./twop1 
TILT3  =  atan2(t3,g3)  *  360/twopl 
TILT4  *  atan2(t4,g4)  *  360/twopl 

WRITE  (4. 119)  t11t1.t11t2.tnt3.tnt4 
119  Format!/,'  POSSIBLE  INSTRUMENT  TILT  (deg)',/, 

1  5x . 'OP  1 ' ,6x . 'DP2 ' . 6x, 'DP3' .6x. 'DP4' ,/,4f9.3 ,/) 

C  Create  acceleration  terms 
DO  56  1*1  .N 

DPXX(I)  »  (0P3( 1 )  -  DP  1 ( I ) )/arm13 
DPVY(I)  *  ( DP4 ( 1 )  -  DP2( 1 )  )/arm24 
56  Continue 

DO  7034  1*1 .600 

WRITE!  14.7033)  dpxx< 1 ) . dpyy (  1 ) 

7034  Continue 
7033  Format(2f 16.9) 

CALL  F F T ( dp  1 . d 1 1 . k , 0 ) 

CALL  FFT(dp2 . d2 1 , k ,0) 

CALL  FFT(dp3.d31 ,k .0) 

CALL  F F T ( dp4 , d4 1 .  k  ,  O ) 

CALL  FFT( dpxx . dxx 1 , k . O ) 

CALL  FFT(dpyy ,dyy 1 , k , 0) 

CALL  FFT( abs . ab 1 , k .0) 


C  Generate  Pressure  Response  Function  and  Wavenumber 
C  using  the  calculated  water  depth.  Sensors  are  average 

C  height  of  3.79  feet  above  bottom. 

bo-. torn  *  depth  +  3.79 
DO  762  1*2. Icut 

sigma  -  twop1/(N*dT)  *  (1-1) 
f  req( 1 )  *  sigma 
call  WVLEN  (bottom, sigma, www) 
wvnr(1)  =  www/bottom 

762  PRSP( 1 )*cosh(wvnr( 1 )*(bottom-sgage) ) 

1  /cosh(wvnr( 1 )*bottom) 


C  Divide  Fourier  coefficients  by  pressure  response 
C  function  and  seawater  specific  gravity. 

CALL  HYDRO( Icut , prsp, wvnr , abs , ab 1 .dpl.dll ,dp2,d21 , 
1  dp3.d31 ,dp4,d41 , dpxx, dxx 1 .dpyy.dyyl .gamma) 

C  Create  power  spectra 
DO  67  1-2, Icut 

ABO)  *  2 ,0*(abs(  1  )**2  +  ab1(  1  )**2) 

01(1)  -  2.0*(dpl( 1 )**2  +  dll( 1  )**2 ) 

02(1)  -  2  0*(dp2( 1 )**2  +  d2 1 ( 1 )**2 ) 

03(1)  *  2.0*(dp3(  1  )**2  d3 1(1  )**2 ) 

D4( 1 )  «  2.0*(dp4( 1 )**2  +  d41(1)**2) 

DXX(  1  )  *  2 ,0*(dpxx( 1 ) **2  +  dxx1( 1 )**2) 


24050 

24100 


1 f (dp4 (J).ne.0.0)tan4=atan2(d41(J), dp4 ( j ) ) 
If (dp4( j ) .eq.O. )tan4  =  0.0 


24150 

tanu 

’  tan2  -  tanl 

24200 

tanv 

=  tan4  -  tan3 

24250 

tanw 

=  tan2  -  tan3 

24300 

tanz 

’  tan4  -  tanl 

24350 

D12  = 

sqr t (dl ( J ) )*sqrt(d2( j ) ) 

♦ 

24400 

034  ’ 

sqrt (d3( J ) ) *sqr t ( d4( j  )  ) 

* 

24450 

023  ’ 

sqrt(d3( j ))*sqrt(d2( j  )) 

* 

24500 

014  ’ 

sqrt(d1 ( J ) ) *sqrt(d4( j ) ) 

* 

24550 

znum 

’  2*012 

24600 

denom 

=  dl ( J  )  -  d2( j  ) 

24650 

24700 

24750 

24800 

24850 

24900 

24950 

25000 

25050 

25100 

25150 

25200 

25250 

25300 

25350 

25400 

25450 

25500 

25550 

25600 

25650 

25700 

25750 

25800 

25850 

25900 

25950 

26000 

26050 

26100 

26150 

26200 

26250 

26300 

26350 

26400 

26450 

26500 

26550 

26600 

26650 

26700 

26750 

26800 

26850 

26900 

26950 

27000 


( cos( tanu)+s 1n( tanu ) ) 
( cos( tanv )+s 1 n( tanv ) ) 
(cos( tanw)+s1n( tanw) ) 
(coa( tanz )+s 1 n( tanz ) ) 


DIRI(J)  ’ORIENT-  ( 180. /3. 14 1 59 ) *0 . 5*a tan2 ( znum , denom ) 
znum  =  2*D34 
denom  =  d3(J)  -  d4(j) 

D I R2 ( j  )  ’ORIENT-  ( 180 . /3 . 14 159 ) *0 . 5*atan2( znum , denom ) 
znum  =  2*D23 
denom  *  d3(J)  -  d2(j) 

DIR3( J  )  ’ORIENT-  (  180. /3 .  14159)*0. 5*atan2(znum, denom) 
znum  =  2*014 
denom  *  d 1 ( J )  -  d4(j) 

DIR4(j)  ’ORIENT-  (  180 . /3 .  14  1 59 ) *0 . 5*atan2 ( znum . denom ) 


5 * a t an2 ( znum , denom ) 


«  dl ( J ) 
“ORIENT- 


-  d4( J  ) 

(  180. /3.  14 159  )  *0  . 


5 *atan2( znum, denom) 


IF(d1r1(J). )t.O. )d1r1(] )=d1r1( J )+360.0 
IF(d1r2(J).lt.O. )d1r2( J )=d1r2( j )+360.0 
IF(d1r3(j).lt.O.)d1r3(J)=d1r3(j )+360.0 
IF(d1r4(J).1t.O.)d1r4(J)=d1r4(j )+360.0 

Continue 

CALL  BLOCK( icut.dirl.maxblk.bik) 

CALL  BLOCK( lcut.dlr2.maxblk.blk) 

CALL  BLOCK( 1cut.d1r3.maxblk.blk) 

CALL  BLOCK ( 1 cut , d 1 r4 , maxbl k ,bl k ) 


Return 

End 


SUBROUTINE  MAXSIG  (R.icut  naxf rq , n . dT , nuf f . f req ) 

Locates  "nuff"  number  of  peak  energy  bands 

and  reports  their  location  and  relative  magnitude. 

Dimension  R(2050),  maxf rq( 20) . f req( 2050) , PD( 20) , 

1  exceed) 20 ) . A ( 2050 ) 

NN  =  Icut  -  1 
TOT  ’  0.0 

Identify  average  energy  In  spectra 
00  440  1*1, Icut 
A(l)  =  R( 1 ) 

440  TOT  *  A ( 1 )  +  tot 
AVRG  •  TOT/nn 

Identify  "nuff*  number  of  peak  bands.  Calculate  each 
peak  band's  magnitude  beyond  the  average.  Set  band  =  0 
after  Identification  as  a  maximum  In  the  temporary  array. 


27050 

DO  445  m= 1 , nuf f 

27100 

maxfrq(m)  =  1 

27  150 

egymax  =  0.0 

27200 

00  450  j»1 ,nn 

27250 

' f! A( j  +  1  ).  1 t . EGYMAX )  GO  TO  450 

27300 

EGYMAX  =  A( j  +  1  ) 

27350 

maxfrq(m)  =  J+1 

27400 

450 

Cont Inue 

27450 

k  =  maxfrq(m) 

27500 

EXCEEO(m)  =  (A(k)-avrg)  /  avrg  +  100.0 

27550 

PD(m)  =  6 . 283 18/f req! k ) 

27600 

A ( k )  =  0.0 

27650 

445 

Cont  1  ruie 

27700 

27750 

WRITE (4, 945) 

27800 

945 

Format  (//,'  Maximum  Energy  located  at  bands:') 

27850 

WRITE ( 4 , 946  )  ( pd( m ) , exceed! m ) , m= 1 . nuf f  ) 

27900 

946 

Format! f 15. 3. '  exceeds  average  energy  by'. 

27950 

1 

f  8 . 2  , '  %' ) 

28000 

28050 

Return 

28100 

End 

28150 

28200 

C 

28250 

28300 

subroutine  FCOEFS  ( ab . abs . ab i . d 1 . dp  1 , d 1 1 , d2 . dp2 . 

28350 

1 

d2 1 , d3 . dp3 . d3 1 . d4 , dp4 . d4 1 , dxx . dpxx , dxx 1 , dyy . 

28400 

2 

dpyy .  dyy  1  .  wvnr  ,  f  req ,  thta  1  .  thta2  ,  thta3 , 

28450 

3 

aO.a1.a2.a3.b1.b2.b3.blk, icut) 

28500 

28550 

c 

Finds  the  first  nine  directional  fourler  coefs 

28600 

c 

for  a  few  selected  frequencies. 

28650 

c 

This  uses  only  one  possible  equation  for  each  coef, 

28700 

c 

but  calculates  each  coeflclent  with  all  possible 

28750 

c 

combinations  of  gages  for  that  equation. 

28800 

28850 

DIMENSION  ab( 2050) . abs! 4 100) . abl ( 2050) . 

28900 

1 

dl (2050). dpi (4 100) ,d1 1(2050), 

28950 

2 

d2 ( 2050 ) , dp2 (41 00 ) , d2 1 ( 2050 ) 

29000 

DIMENSION  d3( 2050 ),dp3(4100),d31(2050). 

29050 

1 

d4 ( 2050 ) , dp4 ( 4 1 OO ) . d4 1 ( 2050 ) . 

29100 

2 

dxx(2050) ,dpxx( 4 100) .dxx 1(2050) . 

29150 

3 

dyy(2050) ,dpyy(4100) ,dyy1(2050) 

28200 

DIMENSION  wvnr (2050) . freq( 2050) 

29250 

DIMENSION  thta1(400).thta2(400),thta3(400) 

29300 

DIMENSION  a0(400) , a  1 ( 400) , a2(400) . a3( 400) , 

29350 

1 

bl ( 400) ,b2( 400) . b3( 400) 

29400 

INTEGER  blk 

29450 

DIMENSION  aldl ( 400) , a1d3(400) ,a2d1d2(400) . 

29500 

1 

a2d1d4(400) ,a2d3d2( 400) ,a2d3d4( 400) ,a3d 1(400) , 

29550 

2 

a3d3( 400) ,b 1d2 ( 400) , b 1d4( 400) , 

29600 

3 

b2d2d1 (400) .  b2d2d3( 400) . b2d4d 1 ( 400) .  b2d4d3(400). 

29650 

4 

b3d2( 400) ,  b3d4 ( 400 ) 

29700 

29750 

pi  «  3.14159 

29800 

twopl  =  6.283185 

29850 

29900 

DO  510  1-2. Icut 

29950 

30000 

222 


42050 

DO  58  1 

*2, 

icut 

42100 

PK  » 

prsp!  1 ) 

* 

gamma 

42150 

abs! 1  ) 

=  abs! 1 )  /  PK 

42200 

abl ( i ) 

-  abl (1 )  /  PK 

42250 

a(1) 

M 

a!  1) 

/ 

PK 

42300 

b!  1  ) 

* 

b!  1) 

/ 

PK 

42350 

c(1  ) 

* 

c(i) 

/ 

PK 

42400 

d(1) 

a 

d(  1 ) 

/ 

PK 

42450 

e!  1) 

3 

e(  1 ) 

/ 

PK 

42500 

f!  1) 

= 

fd> 

/ 

PK 

42550 

g!  i ) 

3 

g!  i ) 

/ 

PK 

42600 

h!  1 ) 

e 

h!  1 ) 

/ 

PK 

42650 

p(  1 ) 

3 

p!  i ) 

/ 

PK 

42700 

q(  1 ) 

3 

q(D 

/ 

PK 

42750 

r(  1 ) 

= 

r!  1 ) 

/ 

PK 

42800 

si  1 ) 

= 

s  ( 1 ) 

/ 

PK 

42850  58 

Continue 

42900 

42950 

Return 

4  3000 

End 

43050 

43100  C 

43150 

43200 

43250 
43300 
43350  C 
43400  C 
43450  C 
43500  C 
43550 

43600  DIMENSION  aO< 400) . a  1 ( 400) . a2( 400) . b3( 400) . b 1 ( 400) , 

43650  1  b2( 400) , b3( 400) , f req( 2050) . maxf rq( 20) , 

43700  2  S(20, 186), deg! 186),pd(20) 

43750 

43800  C  DEGINC  Is  the  Increment  In  degrees  for  which  the 
43850  C  directional  spectra  Is  calculated  over. 

43900  DEGINC  «  2.0 

43950 

44000  C  Weighting  factors  to  smooth  directional  spectra  and 
44050  C  eliminate  negative  side  energy  lobes: 

44100  W1  «  66-/77. 

44150  W2  =*  15. /28. 

44200  W3  =  5./21 . 

44250 

44300  C  Weighting  factors  for  non- smoothed  directional  spectra: 

44350  C  W1  -  1. 

44400  C  W2  »  1 . 

44450  C  W3  *  0. 

44500 

44550 

44600  C  Print  out  the  directional  Fourier  coefficients, 

44650  C  If  desired.  The  header  follows: 

44700 
44750 

44800  438 

44850  1 

44900  2 

44950  3 

45000  4005 


WRITE! 4 , 438 ) 

Format!/, '  PERIOD' . I6x , ' AO' . 14x . ' A 1 ' , 14x, ' A2 ' , 14x . 

'A3' . 14X. '81' . 14x, '82' , 14x, 'B3' ./, '  _ ' . 16x, 

' _ ' , 14x, ' _ ' , 14x, ' _ ' . 14x, ' _ ' , 14x, ' _ ' , 14x, 

' _ ' , 14x, ' _ ' ./) 

Cont Inue 


Subroutine  DIRXSP  ( aO. a  1 ,a2 , a3 . 

1  bl . b2 ,b3 , f req, n,dT ,nuf f , maxf rq, or  lent ) 

We  weight  the  fourler  series  directional  coefs  with 
cos1ne**2s  factors.  Here,  s  =  2.  Then  we  create  a 
matrix  of  the  frequencies  of  interest  analyzed  at 
Incremental  degrees. 


45050 

45100 

00  601  NN» 1 , nuf f 

45150 

save  •  0.0 

45200 

I  ■  maxfrq(nn) 

45250 

PD(nn)  ■  6 . 2831/f req( 1 ) 

45300 

45350 

WRITE (4. 439)  pd( nn) . aO( 1 ) . al ( 1 ) . a2( 1 ) ,a3( 1 ) . 

45400 

1 

b 1 ( 1 ) , b2 ( 1 ),b3(  1 ) 

45450 

439 

Format(f6.2, '  secs' ,7f 16.4) 

45500 

45550 

DO  604  m«1, 181 

45600 

deg(m)  «  save 

45650 

as  •  a0( 1 )  +  a1(1)*w1  *  cosd(deg(m)) 

45700 

1 

♦  a2(1)*w2  *  cosd(2 ,0*deg(m) ) 

45750 

2 

+  a3( 1 )*w3  *  cosd(3.0*deg(m)) 

45800 

bs  ■  b1(1)*w1  *  s1nd(deg(m)) 

45850 

1 

+  b2(l)*w2  *  s1nd(2 .0*deg(m) ) 

45900 

2 

+  b3( 1 ) *w3  •  s1nd( 3 . 0*deg(m) ) 

45950 

46000 

S(nn,m)  ■  ai  ♦  bs 

46050 

46100 

save  •  save  ♦  deg Inc 

46150 

604 

Continue 

46200 

601 

Continue 

46250 

46300 

46350 

C  Determine  the  direction  of  greatest  energy  for 

46400 

C 

each  band  and  convert  to  true  north  heading  from 

46450 

C 

which  waves  propogate. 

46500 

46550 

WRITE (4, 6002) 

46600 

DO  613  NN-1.NUFF 

46650 

mBIG  «  O 

46700 

BIG  «  0.0 

46750 

00  612  *“1,181 

46800 

IF(S(nn,m) .gt .BIG)mBIG  ■  m 

46850 

612 

IF(S(NN,M).gt .BIG)BIG“S(NN,M) 

46900 

BIGOEG  •  mBIG  *  deg  Inc  -  deg  Inc 

46950 

BIGOEG  «  ORIENT  -  BIGOEG 

47000 

IF(BIG0EG. 1 t .0. )BIGOEG-BIGDEG+ 360.0 

47050 

WRITE ( 4 , 976 )  nn.BIGDEG 

47100 

613 

Cont Inue 

47150 

976 

Format ( '  Band', 12,':  Greatest  Energy  at' ,f7.2. 

47200 

1 

'  deg') 

47250 

47300 

C  Record  the  directional  spectra  for  each  band. 

47350 

C 

Directions  correspond  to  angle  from  which  waves 

47400 

C 

propagate  with  respect  to  true  north. 

47450 

47500 

WRITE (3, 969)  (pd(1),  1-l.nuff) 

47550 

00  615  m-1 , 181 

47600 

DEG(m)  •  ORIENT  -  DEG(m) 

47650 

I F ( DEG(m) . 1 t . O . )DEG(m)-DEG(m)+360.0 

47700 

WRITE (3, 970)  deg(m) , (S(nn,m) ,nn“1 ,nuff ) 

47750 

615 

Continue 

47800 

47850 

969 

Format  (//.'  DIRECTIONAL  -  FREQUENCY  SPECTRUM' 

47900 

1 

'  dlrxn' ,9x.6f 17.3, 'secs' ,/) 

47950 

970 

Format  (  f8.1.  '  deg' ,3x.6f 17.3) 

48000 

6002 

Format(/) 

220 
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*€ 


I 


36050 

36100 

CALL  BL0CK( 1  cut . ab , maxb 1  k  ,  b  1  k  ) 

36150 

CALL  BLOCK ( Icut , aO, maxbl k , bl k ) 

36200 

CALL  8L0CK( Icut . a  1 , maxbl k . b 1 k ) 

36250 

CALL  BLOCK (  1  cut ,  a2  , maxblk  .blk ) 

36300 

CALL  BLOCK( Icut . a3, maxblk , blk ) 

36350 

CALL  BLOCK( Icut. bl .maxblk, blk) 

36400 

CALL  BLOCK( Icut, b2 .maxblk. blk) 

36450 

CALL  BLOCK(  Icut ,b3, maxblk .blk ) 

36500 

CALL  BLOCK( Icut . thtal .maxblk.blk) 

36550 

CALL  BLOCK ( Icut . thta2 . maxblk.blk) 

36600 

CALL  BLOCK ( Icut , thta3 .maxblk .blk ) 

36650 

CALL  BLOCK( Icut , freq, maxbl k . blk ) 

36700 

CALL  BLOCK( Icut, dl .maxblk.blk) 

36750 

CALL  BLOCK( Icut ,d2 .maxblk .blk ) 

36800 

CALL  BLOCK( Icut, d3, maxblk, blk) 

36850 

CALL  BLOCK( Icut. d4. maxblk. blk) 

36900 

CALL  BLOCK( Icut .dxx .maxbl k ,bl k ) 

36950 

CALL  BLOCK( Icut. dyy .maxblk.blk) 

37000 

37050 

ICUT  =  MAXblk 

37100 

37150 

37200 

Return 

37250 

End 

37300 

37350 

C  __ 

37400 

37450 

SUBROUTINE  FFT ( ER . F I . K . I CO ) 

37500 

DIMENSION  FR( 4 100 ) , F I ( 4 100) 

37550 

N=2**K 

37600 

IF(ICO.EQ.O)  GO  TO  10 

37650 

DO  8  1  =  1  ,N 

37700 

8 

FI(I)=-FI(I) 

37750 

10 

Cont Inue 

37800 

MR=0 

37850 

NN=N-1 

37900 

DO  2  M= 1 . NN 

37950 

L=N 

38000 

1 

L  =  L/2 

38050 

IF(MR+L.GT.NN)  GO  TO  1 

38100 

MR=MOD( MR , L )+L 

38150 

IF(MR.LE.M)  GO  TO  2 

38200 

TR=FR(M+1) 

38250 

FR(M+ 1 )=FR(MR+1 ) 

38300 

FR( MR+ 1 ) =TR 

38350 

T I =F I ( M+ 1 ) 

38400 

FI(M+1)=FI(MR+1) 

38450 

FKMR+1  )«TI 

38500 

2 

Continue 

38550 

L=  1 

38600 

3 

IF(L.GE.N)  GO  TO  7 

38650 

ISTEP-2*L 

38700 

EL-L 

38750 

DO  4  M-1.L 

38800 

A »3 . 14 15926535* FLOAT ( 1 -M)/EL 

38850 

WR*COS( A ) 

38900 

WI'SIN(A) 

38950 

DO  4  I »M.N, ISTEP 

39000 

vJ-I-tL 

i 


i 
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39050 

I F ( ICO . EO . 1)  GO  TO  11 

39100 

TR=WR*FR(d)-WI*FI(d) 

39150 

T I -WR*  F 1(d) +WI *FR( J ) 

39200 

GO  TO  12 

39250 

1 1 

TR=WR*FR( J)+WI*FI(d) 

39300 

TI=WR*FI(d)-WI*FR(d) 

39350 

12 

FR( d ) =FR( I ) -TR 

39400 

FI(d)=FI(I)-TI 

39450 

FR( I  )  -FR( I )+TR 

39500 

4 

FI(I)=FI(I) +TI 

39550 

L  =  I  STEP 

39600 

GO  TO  3 

39650 

7 

Cont 1 nue 

39700 

AN=N 

39750 

I F ( ICO. EO. 1 )  GO  TO  6 

39800 

DO  5  I  =  1  , N 

39850 

FR( I ) =FR( I )/AN 

39900 

5 

FI( I )  =  —  FI(I )/AN 

39950 

6 

RETURN 

40000 

END 

40050 

40100 

40150 

C 

40200 

40250 

40300 

SUBROUTINE  WVLEN(OPT . SIG. XKH) 

40350 

c* 

THIS  SUBROUTINE  CALCULATES  LINEAR  WAVELENGTH  BY  NEWTONS 

40400 

c* 

METHOD 

40450 

P I =3 . 14159265 

40500 

TWOP 1=2. 0*PI 

40550 

XKH0=SIG**2*DPT/32.2 

40600 

I F ( XKHO-6 . 3 )  2,1.1 

40650 

1 

XKH=XKHO 

40700 

GO  TO  9 

40750 

2 

XKH=SQRT ( XKHO ) 

40800 

3 

SH=SINH( XKH ) 

40850 

CH=COSH(XKH) 

40900 

EPS -XKHO- XKH* SH/CH 

40950 

SL0PE=-XKH/CH**2-SH/CH 

41000 

DXKH=- EPS/SLOPE 

4  1050 

IF (A8S(DXKH/XKH) -0.0001 )  9,9,4 

41100 

4 

XKH= XKH+OXKH 

41150 

GO  TO  3 

4  1200 

9 

XLENTH-TWOPI*DPT/XKH 

4  1250 

RETURN 

4  1300 

end 

4  1350 

4  1400 

c 

41450 

4  1500 

subroutine  HYDRO  ( Icut ,prsp , wvnr , abs , abl , a , b, c . d, e , 

41550 

1 

f,g,h,p,q,r,s. gamma ) 

41600 

c 

Divides  through  by  pressure  response  function 

4  1650 

c 

and  gamma . 

4  1700 

c 

Operates  on  the  a(n)-'  and  b(n)'  terms  from 

4  1750 

c 

each  of  the  five  gages'  FFT  results. 

4  1800 

4  1850 

DIMENSION  prsp(2050) , wvnr (2050) 

4  1900 

DIMENSION  abs(2050),  abi(2050),  a(4 100) ,b( 2050) . 

41950 

1 

c(4 100) ,d( 2050) ,e(4 100) . f ( 2050) . g(4 100) ,h( 2050) , 

42000 

2 

!  p( 4 100) ,q( 2050) , r( 4 10O ) . s( 2050) 

r 


30050 

wvnr(1)  »  wvnr(  0/12.0 

30100 

tanab  ■  atan2(abl ( 1 ) ,abs( 1 ) ) 

30150 

I F ( OP  1(1). NE. 0.0) 

30200 

1 

TAN1  •  atan2(D1 1( I ) ,DP 1 ( I ) ) 

30250 

If (dp1( 1 ) .aq.O.O)  tanl  ■  0.0 

30300 

If (dp2( 1 ) .na.O.O) 

30350 

1 

tan2  •  atan2(d2 1 ( 1 > ,dp2( 1 ) ) 

30400 

1f(d21<1). aq.O.O)  tan2  *  000 

30450 

If  (dp3(  1  )  .na.O.O) 

30500 

1 

tan3  ■  atan2(d31 ( 1 ) ,dp3( 1 ) ) 

30550 

If (dp3( 1 ) -aq.O.O)  tan3  «  0.0 

30600 

if (dp4( 1)  na.O.O) 

30650 

1 

tan4  ■  atan2(d41 ( 1 ) ,dp4( 1) ) 

30700 

If (d41 ( 1 )  aq.O.OO)  tan4  ■  0.00 

30750 

If (dpxx( 1 ) -na.O.O) 

30800 

1 

tanxx  »  atan2(dxx1 ( 1) ,dpxx( 1) ) 

30850 

If (dpxx( 1 ) -aq.O.O)  tanxx  *  0.0 

30900 

If (dpyy( 1 ) .na.O.O) 

30950 

1 

tanyy  ■  atan2(dyy1 ( 1 ) ,dpyy( 1 ) ) 

31000 

1f(dyy1(1 ). aq.O.O)  tanyy  «  0.0 

31050 

31100 

A0( 1 )  -  ab( 1 )  /  (2.0  •  pi) 

31150 

31200 

tanM  ■  tanl  -  tanAB 

31250 

tanN  »  tan3  -  tanAB 

31300 

A1d1(1)  *  -1.0/(pt  *  wvnr(O)  *  aqrt(ab(1))  * 

31350 

1 

sqrt(d1(1))  *  sln(tanM) 

31400 

A 1d3( 1 )  -  -1.0/(p1  *  wvnr(O)  *  sqrt(ab(1))  * 

31450 

1 

sqrt(d3(1))  *  aln(tanN) 

31500 

31550 

A2d1d2( 1 )  -  - ( d2 ( 1 )  -  dl ( 1 ) )/(wvnr( 1 )**2*p1 ) 

31600 

A2dld4( 1 )  -  -<d4(l)  -  d1( 1 ) )/(*ivnr( 1 )**2*pl ) 

31650 

A2d3d2( 1 )  -  -(d2( 1 )  -  d3(1))/(wvnr(f)**2*p1) 

31700 

A2d3d4( 1 )  -  - ( d4 ( 1 )  -  d3( 1 ) )/(wvnr( 1 )**2*p1 ) 

31750 

31800 

C 

GO  TO  4001 

31850 

tanO  ■  tanxx  -  tanl 

31900 

tanP  »  tanxx  -  tan3 

31950 

tanO  »  tanl  -  tanAB 

32000 

tanR  -  tan3  -  tanAB 

32050 

A3d1(<)  ■  4 . 0/ ( p 1  •  wvnr(t)**3) 

32100 

1 

*  (-aqrt(d1( 1 ) )  *  aqrt(dxx(1))  *  aln(tanO) 

32150 

2 

+  0.75  *  wvnr(1)**2  *  aqrt(ab(0)  *  aqrt(dKl)) 

32200 

3 

*  aln(tanO)) 

32250 

A3d3(  1  )  »  4 . 0/ ( p i  *  wvnr( l)**3) 

32300 

1 

*  (-aqrt(d3(  1 ))  *  sqrt(dxx(1))  *  sin(tanP) 

32350 

2 

♦  0.75  *  wvnr( 0**2  *  aqrt(ab(1))  *  sqrt(d3(l)) 

32400 

3 

*  aln(tanR)) 

32450 

4001 

Cont Inua 

32500 

32550 

tans  *  tan2  -  tanAB 

32600 

tanT  •  tan4  -  tanAB 

32650 

B Id2(  1 )  *  - 1 . 0/ ( p  1  *  wvnr(i))  *  aqrt(ab(D) 

32700 

1 

*  sqrt(d2(D)  *  aln(tanS) 

32750 

Bld4(  1 )  •  -1.0/(p1  *  wvnr(1))  *  sqrt(ab(D) 

32800 

1 

*  sqrt(d4( 1 ) )  *  sln( tanT) 

32850 

32900 

tanU  «  tan2  -  tanl 

32950 

tanV  ■  tan2  -  tan3 

33000 

tanW  *  tan4  -  tanl 
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33050 

33100 

33150 

33200 

33250 

33300 

33350 

33400 

33450 

33500 

33550 

33600 

33650 

33700 

33750 

33800 

33850 

33900 

33950 

34000 

34050 

34100 

34150 

34200 

34250 

34300 

34350 

34400 

34450 

34500 

34550 

34600 

34650 

34700 

34750 

34800 

34850 

34900 

34950 

35000 

35050 

35100 

35150 

35200 

35250 

35300 

35350 

35400 

35450 

35500 

35550 

35600 

35650 

35700 

35750 

35800 

35850 

35900 

35950 

36000 


tanX  *  tan4  -  tar>3 
B2d2d1 ( 1 )  -  2 . 0/ ( p 1  *  wvnr( 1 )**2) 

1  *  sqrt(dKO)  *  sqrt(d2(D)  *(cos(  tanu)+s1n(  tanU)  ) 

82d2d3( 1 )  «  2 . 0/ ( p 1  *  wvnr( I )**2) 

1  *  sqrt(d3(1))  *  sqrt(d2(1))  *(coa( tanV)+s1n( tanV ) ) 

B2d4d1 ( 1 )  ■  2.0/(p1  *  wvnr(1)*»2) 

1  *  sqrt(d1(1))  *  sqrt(d4(1))  »(cos( tanW)+sfn( tanW) ) 

B2d4d3( 1 )  *  2 . 0/ ( p 1  *  wvnr  ( 1 )»*2 ) 

1  *  sqrt(d3(1))  *  sqrt(d4(1))  * (cos( tanX )+s 1n( tanX ) ) 

C  GO  TO  4002 

tany2  «  tanyy  -  tan2 
tanz2  «  tan2  -  tanab 
tany4  »  tanyy  -  tan4 
tanz4  *  tan4  -  tanab 
B3d2( 1 )  =  -4.0/(p1  •  wvnr( 1 )**3 ) 

1  *  ( -sqrt(d2( 1 ) )  •  sqrt(dyyd))  *s1n(tany2) 

2  +  0.75  *  wvnr(1)**2  *  sqrt(ab(1))  *  sqrt(d2(1)) 

3  *  s1n(tanz2)) 

B3d4( 1 )  ■  -4.0/(p1  •  wvnr( 1 )**3 ) 

1  *  ( -sqrt(d4( 1 ) )  •  sqrt(dyyd))  *  Sln(tany4) 

2  *  0.75  *  wvnr(1)**2  *  sqrt(ab(1))  *  aqrt(d4(l)) 

3  *  s1n(tanz4)) 

4002  Continue 


510  CONTINUE 


DO  511  1-2. 1 Cut 

C  IF  ALL  GAGES  ARE  WORKING.  USE  THIS  LOOP: 

GO  TO  7 

A1(1)  -  (AldK 1 )  +  A1d3( 1 ))/2.0 

A2( 1 )  ■  ( A2d1d2( 1 )  ♦  A2d1d4( 1 )  ♦  A2d3d2( 1 ) 

1  ♦  A2d3d4( 1 ) )  /  4.0 

A3( 1 )  «  ( A3d1 ( 1 )  +  A3d3( 1 ) )/2 .0 
B 1 ( 1 )  »  (B1d2( 1 )  +  B1d4(1))/2.0 
B2( 1 )  •  (B2d2d1 ( 1 )  +  B2d2d3( 1 )  ♦  B2d4d1(1) 

1  +  B2d4d3( 1 ) )/  4.0 

B3( 1 )  -  (B3d2( 1 )  +  B3U4 ( 1 ) )  /  2.0 
7  Continue 

C  IF  ONLY  TWO  DIFFERENTIAL  CHANNELS  ARE  WORKING.  USE: 

C  GO  TO  3077 

A  1 ( 1 )  •  A1d3( 1 ) 

A2( 1 )  -  A2d3d2( 1 ) 

A3( 1 )  •  A3d3( 1  ) 

B 1 ( 1 )  •  B1d2( 1 ) 

B2( 1 )  •  B2d2d3( 1 ) 

B3( 1 )  •  B3d2( 1 ) 

3077  Continue 


C  Calculate  MEAN  direction  for  n*1,2,3 

THTAI(I)  •  atan2( B 1 ( 1 ) , A  1 ( 1 ) )  *  180. /3. 14159 
THT A2( 1 )  -  atan2(B2( t ) , A2( 1 ) )  *  180. /3. 14159 
THTA3( 1 )  ■  atan2(B3( 1 ) , A3( 1 ) )  *  180. /3. 14159 

511  Continue 
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51050 

51100 

51150 

51200 

51250 

51300 

51350 

51400 

51450 

51500 

51550 

51600 

51650 

51700 

51750 

51800 

51850 

51900 

519SO 

52000 

52050 

52100 

52150 

52200 

52250 

52300 

52350 

52400 

52450 

52500 

52550 

52600 

52650 

S2700 

52750 

52800 

52850 

52900 

52950 

53000 

53050 

53100 

53150 

53200 

53250 

53300 

53350 

53400 

53450 

53500 

53550 

53600 

53650 

53700 

53750 

53800 

53850 

53900 

53950 

54000 


L0W(K1o+ 1 )  *  900. 

HIGH(Kh1+ 1 )  «  900. 

LOSTOP  ■  O 
HI STOP  «  0 

WRI TE (4 , 875 )  NGAGE 

WRITE (4, 890)  Klo.Khl . Avg.Dev , Lol Im.HIl 1m 
c  WR I T  E  (  4 , 900  ) 

GO  TO  400 
DO  400  1-1. N 

I F ( LOW( 1 ) . eq . 900 . )  LOSTOP-1 
IF(HIGH( 1 ) . eq . 900 . )  HISTOP-1 
I F( LOSTOP . eq . 0  .and.  HISTOP.eq.O) 

1  WR I T  E ( 4 , 90 1 )  LOW( 1 ) ,HIGH( 1 ) 

I F ( LOSTOP . eq . 1  .and.  HISTOP.eq.O) 

1  WRI TE ( 4 , 902 )  H1GH( 1 ) 

I F ( LOSTOP . eq . 0  .and.  HISTOP.eq.1) 

1  WRITE ( 4 , 903 )  LOW( 1 ) 

IF(LOSTOP.eq. 1  .and.  HISTOP.eq.1) 

1  GO  TO  410 

400  Continue 
410  CONTINUE 

C  Average  2  adjacent  points  to  halve  the  record  size. 
N2  «  N/2 
DO  480  1-1. N2 
d-I-2-1 

X(D  -  (X(j)+X(J+1))/2.0 
480  CONTINUE 

C  Re-ca leu late  mean. 

405  SUM  -  0.0 

DO  440  I-1.N2 
440  SUM  -  X{ 1 )  ♦  SUM 

AVG  -  SUM/N2 

C  Subtract  mean  from  record  and  convert  from  volts 
C  to  PSI.  Divide  dlff  records  by  distance  between 
C  sensing  points  to  give  pressure  differential 
C  per  length  of  arm. 


DO  450  1-1. N 

I F ( NGAGE . ne . 5 )  X(1)  -  ( (X( 1 >-AVG)/5 . 0)/DXINCH 
IF (NGAGE .eq.5)  X(1)  -  (X(l)-avg)  *  10.0 


450 

Continue 

875 

Format(/f 10*. 'GAUGE  NUMBER'. 12) 

890 

Format( 16. ’  LOW  truncations; '. 16, '  HIGH 

truncations' ,/ 

1 

4x. 'Mean-' .f6.3. '  S.Dev-' .fB.5. ' 

Lower  Limit-' , 

2 

f6.3, '  Upper  Limit-' ,f6. 3) 

900 

Format (4x, 'LOW' .7x. 'HIGH' ) 

901 

Format ( F8 . 2 , F 10 . 2 ) 

902 

Format (8x , F 10. 2 ) 

903 

Format (F8 . 2 ) 

Return 

End 
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54050 

C 

54100 

54150 

54200 

SUBROUTINE  BL0CK( n , X , maxb 1 k . bl k ) 

54250 

DIMENSION  X ( 4 10O) 

54300 

INTEGER  BLK.STOP 

54350 

54400 

STOP  »  (N-BLK) 

54450 

NM  »  1 

54500 

54550 

DO  100  10*2, stop.blk 

54600 

SUM-0.0 

54650 

Md= Id+BLK- 1 

54700 

Kd-Id 

54750 

DO  130  KKU=k j , mj 

54800 

130 

SUM-SUM  +  X(kkj) 

54850 

X(nm)  -  SUM/BLK 

54900 

NM-NM+1 

54950 

lOO 

Continue 

55000 

55050 

SUM-0 . 0 

55100 

DO  140  kkj-1 J ,n 

55150 

140 

SUM-SUM  ♦  X(kkj) 

55200 

X(nm)  «  SUM/(n-mj) 

55250 

55300 

MAXBLK-nm 

55350 

55400 

Return 

55450 

End 

55500 

55550 

55600 

SUBROUTINE  SURE2(ab.d1  .  d?. , d3 , d4 ,  angl ,  ang2  ,  ang3 ,  ang4 ,  Icut 

55650 

1 

wvnr.blk, orient) 

55700 

55750 

DIMENSION  ab(2050) ,d1( 2050 ) . d2 ( 2050 ) , d3 ( 2050 ) , d4 ( 2050 ) . 

55800 

1 

angl ( 400) ,ang2(4001 ,ang3(400) ,ang4(400) , 

55850 

2 

wvnr ( 2050 ) 

55900 

INTEGER  blk 

55950 

56000 

c 

Estimates  direction  using  only  absolute  gage  and  one 

56050 

c 

differential  channel  at  a  time. 

56100 

56150 

pi  =  3. 14159 

56200 

56250 

DO  300  J-2. Icut 

56300 

z  -  1./wvnr(J)  *  sqrt(d1( J )/ab( J ) ) 

56350 

I F( 2 . LE . 1 . ) ANG1 ( J ) -ORIENT-acos( 2 )* 180. /pi 

56400 

2  -  1./wvnr(J)  *  sqrt(d1( J )/ab( j ) ) 

56450 

IF  (  2 . LE . 1 . ) ANG2( J ) -OR I  ENT -as  1 n( 2) *  180. /pi 

56500 

2  -  1 ,/wvnr( J )  •  sqrt (d3( J )/ab( J ) ) 

56550 

IF(2.LE. 1 . )ANG3( j )-ORIENT-acos( Z )* 180 . /pi 

56600 

2  «  1./wvnr(J)  *  sqrt (d4( J )/ab( J ) ) 

56650 

IF( 2 . LE . 1 . )ANG4( j )«ORIENT-as 1n( Z )* 180. /pi 

56700 

56750 

IF (angl ( J ) . 1 t .0. )ang1( J )-ang1( J )+360.0 

56800 

IF(ang2( J ) . 1 t .O. )*ng2( j )«ang2( J ) +360.0 

56850 

IF(ang3( J ) . 1 t .0. )artg3( J )«ang3( J ) +360.0 

56900 

IF(ang4( J ) . 1 t .0. )ang4( j )*ang4( J ) +360.0 

56950 

57000 

300 

Continue 
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57050 

57100  CALL  BLOCK ( Icut , angl , maxbl k , bl k ) 

57150  CALL  BLOCK( (cut . ang2 .maxblk . bl k ) 

57200  CALL  BLOCK( Icut , ang3 . maxblk . bl k ) 

57250  CALL  BLOCK{ 1 cut . anfl4 , maxbl k , b ) k ) 

57300 

57350  Return 

57400  End 

57450 
57500  C 
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